Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,390,012 Δ1 bp 16.1% intergenic (+48/+36) yfbP → / ← nuoN TPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,390,0120G.16.1% 473.0 / 92.4 161intergenic (+48/+36)yfbP/nuoNTPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base G (62/73);  new base . (13/13);  total (75/86)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-02

Only 100 of 166 total aligned reads displayed.

CAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACC  >  NC_000913/2389913‑2390111
                                                                                                   |                                                                                                      
cAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAg                                                                                                        <  1:767803/100‑1 (MQ=255)
cAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAg                                                                                                        >  1:986268/1‑100 (MQ=255)
cAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAg                                                                                                        <  2:2500905/100‑1 (MQ=255)
 aGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGc                                                                                                       >  2:1482255/1‑100 (MQ=255)
   gCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCAt                                                                                                     <  1:2272080/100‑1 (MQ=255)
    cAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATg                                                                                                    <  1:1359555/100‑1 (MQ=255)
    cAGTTATGTATAGCGTATTTTGGCGATCGTCGTAAACCGTTCACTTAACGCTCTCCGTATCTCTTCTCACCCTTCTTTTTCACAAAAGCCGCCGAGCATg                                                                                                    <  2:3039529/100‑1 (MQ=255)
      gTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTc                                                                                                 <  1:1483664/100‑1 (MQ=255)
          tGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcggcg                                                                                             <  1:338367/100‑1 (MQ=255)
           gTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTcggcgg                                                                                            >  1:1072952/1‑100 (MQ=255)
            tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAGAAACCGCCGACATGCTcggcggt                                                                                           <  1:2213057/100‑2 (MQ=255)
            tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCggcggc                                                                                          <  2:2824414/100‑1 (MQ=255)
             ttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTcggcggtt                                                                                          <  1:1534585/100‑3 (MQ=255)
               agagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTCGGCGGtttt                                                                                      <  1:2148358/100‑1 (MQ=255)
                gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTc                                                                                     <  1:3169085/100‑1 (MQ=255)
                gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTc                                                                                     <  1:1319955/100‑1 (MQ=255)
                 agTATTTTGGCGATAGTCTTACACATTTCACTTAACTCTCTGCGTCTATCTACTATACCTTCTTTTTCAGAAAAACCGCCGAGTATGTCGGCGGCTTTc                                                                                     <  2:2308547/99‑1 (MQ=255)
                 agTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCt                                                                                    <  2:401269/100‑1 (MQ=255)
                 agTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCt                                                                                    <  1:3259543/100‑1 (MQ=255)
                 agTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGACATGCTCGGCGGtttttc                                                                                    <  1:2727489/100‑3 (MQ=38)
                  gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTc                                                                                     <  2:307982/98‑1 (MQ=255)
                   tATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCt                                                                                    <  2:3555650/98‑1 (MQ=255)
                   tATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGa                                                                                  <  1:1535399/100‑1 (MQ=255)
                    aTTGTGGCGATAGTCGTACACAGTTCCCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGAc                                                                                 <  2:3230978/100‑1 (MQ=255)
                     ttttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACt                                                                                >  2:3578190/1‑100 (MQ=255)
                      tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACtt                                                                               <  1:1691541/100‑1 (MQ=255)
                      tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAGCCGCCGACATGCTCGGCGGtttttctgaaaa                                                                              <  1:2071804/100‑9 (MQ=18)
                      tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAGCCGCCGACATGCTCGGCGGtttttctgaaaa                                                                              <  1:1178273/100‑9 (MQ=21)
                       ttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTa                                                                              <  2:3393522/100‑1 (MQ=255)
                        tGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTAc                                                                             <  1:1735733/100‑1 (MQ=255)
                          gCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaa                                                                           <  2:2704286/100‑1 (MQ=255)
                           cGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaaa                                                                          >  1:2635758/1‑100 (MQ=255)
                           cGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACaaa                                                                          <  1:3125802/100‑1 (MQ=255)
                              tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                       <  1:2518403/100‑1 (MQ=255)
                              tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                       <  2:2669169/100‑1 (MQ=255)
                              tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                       <  2:1273005/100‑1 (MQ=255)
                              tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTa                                                                       <  2:1150408/100‑1 (MQ=255)
                              tAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTATTCAGCAAACCCGCCAAGCATGTCGGCGGCTTTCTAACTTACAAAGTa                                                                       >  2:700722/1‑100 (MQ=255)
                                gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTaga                                                                     >  2:1731081/1‑98 (MQ=255)
                                gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAAc                                                                     >  1:122922/1‑100 (MQ=255)
                                gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAAc                                                                     <  2:2180652/100‑1 (MQ=255)
                                gTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAAc                                                                     >  1:119219/1‑100 (MQ=255)
                                   gTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGa                                                                  >  2:686451/1‑100 (MQ=255)
                                    tACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAt                                                                 >  1:1101994/1‑100 (MQ=255)
                                     acacAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGAtt                                                                >  1:660280/1‑100 (MQ=255)
                                       acacTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTAc                                                              <  2:905452/100‑1 (MQ=255)
                                       acaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTAc                                                              <  1:1698318/100‑1 (MQ=255)
                                        caGTTCACTTAACGCTCTGCGTCTAGCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACa                                                             <  2:1623518/100‑1 (MQ=255)
                                              aCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCg                                                       <  1:802016/100‑1 (MQ=255)
                                                ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGc                                                     >  1:3606805/1‑100 (MQ=255)
                                                ttAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGc                                                     >  2:3276239/1‑100 (MQ=255)
                                                 tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                    <  1:1588515/100‑1 (MQ=255)
                                                 tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCACCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGc                                                     >  2:2985154/1‑99 (MQ=255)
                                                 tAACGCTCTGCGTCTATCTACTACCCCTTCTTTTCCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                    <  2:599730/100‑1 (MQ=255)
                                                     gCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGc                                                >  2:1450099/1‑100 (MQ=255)
                                                      ctctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGcc                                               >  1:959586/1‑100 (MQ=255)
                                                        ctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaa                                             >  2:2904520/1‑100 (MQ=255)
                                                          gCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAAc                                           >  2:2834748/1‑100 (MQ=255)
                                                            gTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGc                                         <  1:1931021/100‑1 (MQ=255)
                                                              ctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGcac                                       <  2:1798013/100‑1 (MQ=255)
                                                                atctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACaa                                     >  2:289067/1‑100 (MQ=255)
                                                                 tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAAt                                    <  2:2155875/100‑1 (MQ=255)
                                                                   tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGc                                  >  1:583112/1‑100 (MQ=255)
                                                                   tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGc                                  >  1:144282/1‑100 (MQ=255)
                                                                    actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                 >  2:352087/1‑100 (MQ=255)
                                                                    actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                 >  1:3223562/1‑100 (MQ=255)
                                                                       aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATACTAAt                              >  2:1328120/1‑100 (MQ=255)
                                                                       aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATACTAAt                              >  2:1326301/1‑100 (MQ=255)
                                                                          ccTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAg                           >  2:2783457/1‑100 (MQ=255)
                                                                           cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGTCATTGCTAAACGCACAATGCTAATCAg                           >  2:430868/1‑99 (MQ=255)
                                                                           cttcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGc                          >  1:3628456/1‑100 (MQ=255)
                                                                           cttcttTTTCAGAAAAACCGCCGAGCATGTAGGCGGCTTTCTAAATTACAAAGTAAAACATTACATCAGCGGCATTGACAAACAAAAAATGCTaaaaaac                          >  1:3039852/1‑95 (MQ=255)
                                                                             tcttTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCgg                        <  2:1349989/100‑1 (MQ=255)
                                                                               tttttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGACATTGCCAAACGCACAATGCTAATCAGCGGtt                      >  2:512842/1‑100 (MQ=255)
                                                                                ttttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtg                     >  1:2328708/1‑100 (MQ=255)
                                                                                  ttCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgtg                   <  2:1700498/100‑1 (MQ=255)
                                                                                  agtcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                    >  2:116329/11‑100 (MQ=21)
                                                                                  agtcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                    >  1:1968711/11‑100 (MQ=21)
                                                                                   tCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCgg                        <  2:2000393/94‑1 (MQ=255)
                                                                                   gtcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgtg                   >  1:889657/10‑100 (MQ=21)
                                                                                    tcagaaagccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                  >  1:1281711/9‑100 (MQ=21)
                                                                                    gaGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATAGCCAAACGCACAATGCTAATCAGCGGTCGTGGc                 >  2:648138/2‑100 (MQ=255)
                                                                                    cggAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                 >  1:2912965/3‑100 (MQ=255)
                                                                                    cAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                 >  1:2200413/1‑100 (MQ=255)
                                                                                    cAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                 >  1:3121419/1‑100 (MQ=255)
                                                                                     aGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGcc                <  2:2436750/100‑1 (MQ=255)
                                                                                     aGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                 <  2:1911299/99‑1 (MQ=255)
                                                                                      agaaagccgccgACATGCTCGGCGGTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCa               <  1:2074685/94‑1 (MQ=37)
                                                                                          agccgccgACATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATac            >  1:806171/3‑100 (MQ=37)
                                                                                          aaCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATaca           >  1:2329023/1‑100 (MQ=255)
                                                                                           aCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac          >  2:2326504/1‑100 (MQ=255)
                                                                                           aCCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATACTAATCAGCGGTTGTGGCCATacac          <  1:1328128/100‑1 (MQ=255)
                                                                                            ccgccgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAcc         <  1:304396/100‑1 (MQ=255)
                                                                                            ccgccgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAcc         >  1:1932121/1‑100 (MQ=255)
                                                                                             cgccgAGCATGTCGGCGGCTTTCTGACTTACCAAGTAACAGATTACGTCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAccc        <  2:2550528/100‑1 (MQ=255)
                                                                                             cgccgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACAccc        >  2:2282810/1‑100 (MQ=255)
                                                                                                cgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc     >  2:1821411/1‑100 (MQ=255)
                                                                                                cgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc     <  1:3566177/100‑1 (MQ=255)
                                                                                                cgAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCCGCATGGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc     <  2:619369/100‑1 (MQ=255)
                                                                                                   gCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAcc  <  1:3025179/100‑1 (MQ=255)
                                                                                                   |                                                                                                      
CAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATGTCGGCGGCTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACC  >  NC_000913/2389913‑2390111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: