Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,853 C→G 6.8% V246L (GTA→CTA)  fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,8530CG 6.8% 196.7 / 6.3 59V246L (GTA→CTA) fumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base G (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

ATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCC  >  NC_000913/1685707‑1685997
                                                                                                                                                  |                                                                                                                                                
atcatTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTacgc                                                                                                                                                <  5:90394/149‑1 (MQ=255)
    ttGGCGATTTTCATCAGTGCCGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGGGCACAGGTAATGACTGCCAGTTCATCTGCTacgcgccg                                                                                                                                            >  5:304827/1‑149 (MQ=255)
      ggCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTCACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTccgcgacgcg                                                                                                                                          >  6:340680/1‑149 (MQ=255)
      ggCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTacgcgacgc                                                                                                                                           <  8:383712/148‑1 (MQ=255)
       gCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTacgcgacgcg                                                                                                                                          >  7:142026/1‑148 (MQ=255)
       gCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTacgcgacgcg                                                                                                                                          <  6:285485/148‑1 (MQ=255)
        cGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCa                                                                                                                                        <  3:275714/149‑1 (MQ=255)
          aTTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATa                                                                                                                                      <  3:122509/149‑1 (MQ=255)
              tcatcaGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTcc                                                                                                                                  >  1:235869/1‑149 (MQ=255)
                 tcaGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGa                                                                                                                               <  2:93166/149‑1 (MQ=255)
                   aGTTACGCAGCGAACCCTTTCTACGCGCGGTGCGCCTTAACCAGTTCATCACTGGTCGCGAGCGCTTCAATTTTGTTCGGCGCGGTAACAAACGGTGCACTGGTAATGTCTGCCAGTTCATCTGCTACTCGAAGCTCATAATCCGGATg                                                                                                                             <  8:233134/149‑1 (MQ=255)
                    gTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATgg                                                                                                                            <  2:116397/149‑1 (MQ=255)
                    gTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATgg                                                                                                                            <  7:343435/149‑1 (MQ=255)
                       aCGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATTGGt                                                                                                                          <  2:155523/148‑1 (MQ=255)
                       aCGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTa                                                                                                                         <  8:142026/149‑1 (MQ=255)
                         gCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTCACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTAtt                                                                                                                       >  6:353297/1‑149 (MQ=255)
                         gCAGCCAACCCTTTCAACGCGCCGTGCGCCTCAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACACACGGTGCACAGGTAATGACGGCCAGTTCATCTGCTACGCGACGCGCATACTCCGCATGCGTAtt                                                                                                                       >  5:206181/1‑149 (MQ=255)
                           aGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACATGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTAttca                                                                                                                     <  8:40696/149‑3 (MQ=255)
                              catcccTTTCAAGGCGCCGTGCGCATGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTc                                                                                                                  <  6:137533/146‑1 (MQ=255)
                                  ccTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAg                                                                                                               >  1:190841/1‑148 (MQ=255)
                                   cTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTa                                                                                                             <  1:394010/149‑1 (MQ=255)
                                   cTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTa                                                                                                             <  8:57888/149‑1 (MQ=255)
                                   cTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTGCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTa                                                                                                             <  7:274063/149‑1 (MQ=255)
                                    tttCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTa                                                                                                                     <  6:9061/140‑1 (MQ=255)
                                    tttCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTa                                                                                                                     >  5:9061/1‑140 (MQ=255)
                                           gcgcCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGc                                                                                                     >  1:361640/1‑149 (MQ=255)
                                             gcCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTg                                                                                                   <  4:143823/149‑1 (MQ=255)
                                               cGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTa                                                                                                 <  2:165244/149‑1 (MQ=255)
                                                 tGCGCCTGACCCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTAcc                                                                                               >  6:295235/1‑149 (MQ=255)
                                                    gcCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGc                                                                                             >  8:55987/1‑148 (MQ=255)
                                                      cTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCaa                                                                                          <  4:11734/149‑1 (MQ=255)
                                                         aaCCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAaga                                                                                        <  2:235869/148‑1 (MQ=255)
                                                          aCCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGc                                                                                      >  8:65372/1‑149 (MQ=255)
                                                          aCCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGc                                                                                      <  2:361640/149‑1 (MQ=255)
                                                              gggCATCACAGGTCGCCAGCGCTTCAAATTTGTACGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTAGGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCaa                                                                                          >  8:255747/1‑141 (MQ=255)
                                                              gggCATCACAGGTCGCCAGCGCTTCAAATTTGTACGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTAGGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCaa                                                                                          >  8:250018/1‑141 (MQ=255)
                                                              gggCATCACAGGTCGCCAGCGCTTCAAATTTGTACGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTAGGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCaa                                                                                          <  7:255747/141‑1 (MQ=255)
                                                              gggCATCACAGGTCGCCAGCGCTTCAAATTTGTACGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTAGGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCaa                                                                                          <  7:250018/141‑1 (MQ=255)
                                                                  aTCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTAtt                                                                                                                       <  1:230257/108‑1 (MQ=255)
                                                                  aTCACAGGTCGCCAGCGCTTCAAATTTGTTCGGAGCGGTAACAAAAGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTAtt                                                                                                                       >  2:230257/1‑108 (MQ=255)
                                                                   tCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGc                                                                             >  5:234535/1‑149 (MQ=255)
                                                                    cacaGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGc                                                                             >  2:154502/1‑148 (MQ=255)
                                                                         gTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGAATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGt                                                                        >  8:286016/1‑148 (MQ=255)
                                                                                gcgcTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGaggcagg                                                                <  2:190841/149‑1 (MQ=255)
                                                                                  gcTTCAAATTTGTTCGGCGCGGTAACAAACGGTTCACAGGTAATGACTGGCAGTTCATCTGCTACGCGACGCGCATACTCCTGATGGGTATTTTGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCt                                                              <  6:206181/149‑1 (MQ=255)
                                                                                    ttCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCTACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGGGAGGCAGGCTGt                                                            >  2:99359/1‑149 (MQ=255)
                                                                                     tCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTa                                                           >  3:50240/1‑149 (MQ=255)
                                                                                         aTTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTc                                                        >  2:364247/1‑148 (MQ=255)
                                                                                                   cgcgGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTga                                             <  6:234535/149‑1 (MQ=255)
                                                                                                   cgcgGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTga                                             <  3:8348/149‑1 (MQ=255)
                                                                                                          aCAAACGGTGCACAGGTAATGACTGCCAGTGCATCTGCTACGCGACGCGCATACTCCGGAGGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATg                                      <  5:353297/149‑1 (MQ=255)
                                                                                                             aaCGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTc                                   <  7:55987/149‑1 (MQ=255)
                                                                                                              aCGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGTGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTGTGAGATTATGCTc                                   >  6:81603/1‑148 (MQ=255)
                                                                                                                   gCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCAt                             >  7:180443/1‑149 (MQ=255)
                                                                                                                       aGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCt                         <  4:302062/149‑1 (MQ=255)
                                                                                                                         gTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTAc                       >  2:210373/1‑149 (MQ=255)
                                                                                                                                   ccAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGa              >  4:111222/1‑148 (MQ=255)
                                                                                                                                   ccAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGa              >  3:101056/1‑148 (MQ=255)
                                                                                                                                   ccAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCACTACCCAGCCGGaa             <  5:295235/149‑1 (MQ=255)
                                                                                                                                    cAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGaaa            <  7:65372/149‑1 (MQ=255)
                                                                                                                                           tCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTg     >  7:37624/1‑149 (MQ=255)
                                                                                                                                              gCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGccc  >  7:174055/1‑149 (MQ=255)
                                                                                                                                              gCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGccc  >  1:242072/1‑149 (MQ=255)
                                                                                                                                               cTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGccc  >  2:27301/1‑148 (MQ=255)
                                                                                                                                               cTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGccc  >  4:278702/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCATTGGCGATTTTCATCAGTGACGCAGCCAACCCTTTCAACGCGCCGTGCGCCTGAACCAGGGCATCACAGGTCGCCAGCGCTTCAAATTTGTTCGGCGCGGTAACAAACGGTGCACAGGTAATGACTGCCAGTTCATCTGCTACGCGACGCGCATACTCCGGATGGGTATTTAGTCCAGTACCCACCGCTGTACCGCCAAGAGCCAGTTCCGCTACGTGAGGCAGGCTGTATTCGATATGTTTGAGATTATGCTCGAGCATCGCTACCCAGCCGGAAATCTCCTGCCC  >  NC_000913/1685707‑1685997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑