Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,357,848 A→G pseudogene (33/186 nt) ypaA → protein YpaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,357,8480AG98.8% 262.6 / ‑4.6 83pseudogene (33/186 nt)ypaAprotein YpaA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (51/31);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTG  >  NC_000913/2357721‑2357983
                                                                                                                               |                                                                                                                                       
gcATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATtggc                                                                                                                             <  2:244709/140‑1 (MQ=255)
    gCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATtggctggc                                                                                                                         >  1:1054056/1‑140 (MQ=255)
        gAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAg                                                                                                                       <  1:110792/138‑1 (MQ=255)
        gAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAg                                                                                                                       >  2:110792/1‑138 (MQ=255)
         aGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGTTGCAGCAGGAAGTGCATCAAATTGGCTGGCAGGaa                                                                                                                    >  2:596194/1‑140 (MQ=255)
         aGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                                    >  2:1025087/1‑140 (MQ=255)
         aGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGaa                                                                                                                    >  2:1055257/1‑140 (MQ=255)
         aGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACAGCTGCGGGAGGAAGAGCATCAAATTGGCTGGCAGGaa                                                                                                                    >  1:413277/1‑140 (MQ=255)
          gCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGTAAGGACATCAAATTGGCTGGCAGGAAg                                                                                                                   <  1:1007902/140‑1 (MQ=255)
            aGGGCTTTGAACGTGAGATTGTTCTGGCGACTACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                 <  2:532263/140‑1 (MQ=255)
            aGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                 <  2:1255652/140‑1 (MQ=255)
            aGGGCTTTGAACGTGAGATTGTGCTAGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGACCGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                 <  2:1084456/140‑1 (MQ=255)
            aGGGCTTTGAACGTGAGATTGTCCTGGCGACTACCCAACCCACTGATGTTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGt                                                                                                                 <  2:413277/140‑1 (MQ=25)
              ggCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTaa                                                                                                               >  2:1074266/1‑140 (MQ=255)
              ggCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTaa                                                                                                               >  1:242787/1‑140 (MQ=255)
                  ttGAACGTGAGTTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTa                                                                                                           >  1:1328991/1‑140 (MQ=255)
                  ttGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCCGGAAGGGGATCAAAATGGCTGGCAGGAAGGTAAATTa                                                                                                           >  1:572918/1‑140 (MQ=255)
                   tGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAAGGGCTGCGGTAGGGAGGGCATCAAATGGGCTGGCAGGAAGGTAAATTAg                                                                                                          >  2:424674/1‑140 (MQ=37)
                   tGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATGGGCTGGCAGGAAGGTAAATTAg                                                                                                          >  1:632984/1‑140 (MQ=255)
                         tGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGtt                                                                                                    >  2:1487782/1‑140 (MQ=255)
                          gagaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGttt                                                                                                   >  1:228/1‑140 (MQ=255)
                          gagaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGttt                                                                                                   >  2:1241545/1‑140 (MQ=255)
                          gagaTTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGAATCAAATAGGCGGGTAGGAAGGTAAATTAGAAGGttt                                                                                                   >  1:433767/1‑140 (MQ=37)
                                 tgCTGGCGACAACCCTACTCACTGACGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGaa                                                                                            <  1:526898/140‑1 (MQ=255)
                                 tgCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGaa                                                                                            >  1:283308/1‑140 (MQ=255)
                                     ggCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGGTGGCAGGGAGGCAAGTTAGAACGTTTGCATGAACAAg                                                                                        >  1:713503/1‑140 (MQ=25)
                                      gCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCTTTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGc                                                                                       >  2:975752/1‑140 (MQ=255)
                                      gCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGc                                                                                       <  1:1074266/140‑1 (MQ=255)
                                      gCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCCGGCAGGAAGGTGAACTAGAAGGTTTGCATGAACAAGc                                                                                       >  1:588454/1‑140 (MQ=255)
                                         aCAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCAt                                                                                    <  2:1230460/140‑1 (MQ=255)
                                           aaCCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTa                                                                                  >  2:1072245/1‑140 (MQ=255)
                                              ccAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATAAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATAAAAAGGTTTGCATGAACAAGCCATTaaaa                                                                               <  2:328799/140‑1 (MQ=255)
                                               cAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAAt                                                                              <  1:850177/140‑1 (MQ=255)
                                                  cTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGc                                                                           <  1:1487782/140‑1 (MQ=255)
                                                     aCTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCttt                                                                        >  1:508406/1‑140 (MQ=255)
                                                       tGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGc                                                                                       <  1:1421695/123‑1 (MQ=255)
                                                       tGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGc                                                                                       >  2:1421695/1‑123 (MQ=255)
                                                        gATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTgcg                                                                     <  2:283308/140‑1 (MQ=255)
                                                          tGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCa                                                                   <  1:223701/140‑1 (MQ=255)
                                                             tGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGc                                                                <  2:508406/140‑1 (MQ=255)
                                                             tGATATTCCGAACTGTAGTTAATAGGAGAGTATAACGACTATTGCTGACCGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGc                                                                <  2:713503/140‑1 (MQ=25)
                                                               atatTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTg                                                              >  1:1158302/1‑140 (MQ=255)
                                                                  ttCGGACCTGTCATTAAAAGGAGAGATTAATAAGTATGGCTGAACGGCTGCGTCAGGAAGGGCATCAAATTGGCTGGTAGGAAGATAAATTAGAGGGTTTCCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                           <  2:1328991/140‑1 (MQ=25)
                                                                  ttCCGAACTGTCATTAAAAGGAGAGCTTATTGACTATGGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                           <  1:558920/140‑1 (MQ=255)
                                                                  ttCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                           <  2:228/140‑1 (MQ=255)
                                                                  ttCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGaa                                                           <  2:242787/140‑1 (MQ=255)
                                                                   tCCGAACTGTCATTACAAGGAGCGACTAATGACTATTCCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                          <  1:424674/140‑1 (MQ=37)
                                                                   tCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAAc                                                          <  1:838561/140‑1 (MQ=255)
                                                                   tCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTCGCAGGAAGGTAAATTATAAGGTTTGGATGGACAGGCCATTAACATTGCTTTGCGCATGCTGGAAc                                                          >  2:573420/1‑140 (MQ=25)
                                                                     cGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTCGCGCATGCTGGAACAg                                                        >  1:775817/1‑140 (MQ=255)
                                                                       aaCTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGTTTTGCATGAACAAGACCTTAAAATTACTTTACGCATGCTCGAACAggg                                                      >  2:754143/1‑140 (MQ=21)
                                                                       aaCTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg                                                      >  1:117175/1‑140 (MQ=255)
                                                                       aaCTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg                                                      >  1:960064/1‑140 (MQ=255)
                                                                       aaCTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAggg                                                      >  1:845068/1‑140 (MQ=255)
                                                                          ccacaacaaaaGGAGAGATTAATGACTATTGCTGAAGGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCtt                                                   <  2:53683/133‑1 (MQ=255)
                                                                           gTCATTAAAAGGAGAGAATAATCACTATTGCGGAAGGGCAGCGGCAGGGAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCttt                                                  <  2:632984/140‑1 (MQ=25)
                                                                            tCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTg                                                 >  2:971912/1‑140 (MQ=255)
                                                                               ttAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCACCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGAATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATc                                              <  1:1025087/140‑1 (MQ=255)
                                                                                 aaaaGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAAGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGGAGGTTTGCATGAACAAGCCATAAAAGTTGGCTTGCGCAGGCTGTAACAGGGCTTTGATCGt                                            >  2:448155/1‑140 (MQ=21)
                                                                                      gagagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATTAAATTGGCTGGCAGTAAGGTTAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGTGCATGCTGGAACAGGGCTTTGATCGTGACCa                                       >  1:377194/1‑140 (MQ=25)
                                                                                      gagagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTTCTTTGCGCATGATGAAACAGGACTTTGATCGTGACCa                                       >  1:648001/1‑140 (MQ=25)
                                                                                      gagagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                       >  1:1199042/1‑140 (MQ=255)
                                                                                      gagagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                       >  1:1331031/1‑140 (MQ=255)
                                                                                      gagagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCa                                       >  1:1133747/1‑140 (MQ=255)
                                                                                        gagaTTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAgg                                     >  1:1061245/1‑140 (MQ=255)
                                                                                            taaaTGACTATTGCTTTATGGCTGCGGAGGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGATGGTTGGCGTGACCAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCt                                 <  2:377194/138‑1 (MQ=25)
                                                                                                 gattATTGCTGACCGGCAGCGGAAGAAAGGCCATCAATTTGGCTGGCAGAAAGGTAATTAAGAGGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCgg                            <  2:1543469/137‑1 (MQ=25)
                                                                                                    tATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCa                         <  2:845068/140‑1 (MQ=255)
                                                                                                    tATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCa                         <  2:1061245/140‑1 (MQ=255)
                                                                                                    tATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCa                         <  1:975752/140‑1 (MQ=255)
                                                                                                    tATTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCa                         <  2:117175/140‑1 (MQ=255)
                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGTAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAc                        >  1:969653/1‑140 (MQ=255)
                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGGTTTCATCGTGACCTGGTGCTCGCGGCCAc                        >  1:710176/1‑140 (MQ=255)
                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAc                        >  1:377300/1‑140 (MQ=255)
                                                                                                     aTTGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAc                        <  2:1158302/140‑1 (MQ=255)
                                                                                                      ttGCTGAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCAcc                       <  2:960064/140‑1 (MQ=255)
                                                                                                          tGAACGGCTGCGGCATGAAGGGTATCAAAGTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGc                   <  1:754143/140‑1 (MQ=255)
                                                                                                           gAACGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCt                  >  2:958870/1‑140 (MQ=255)
                                                                                                            aaCGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTa                 >  1:1122593/1‑140 (MQ=255)
                                                                                                            aaCGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTa                 >  1:1174593/1‑140 (MQ=255)
                                                                                                            aaCGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATACTGGAACAGGGCCTTGATCGTGACCCGGGGATCGCGGCCACCCAGCTa                 >  2:616434/1‑140 (MQ=37)
                                                                                                              cGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCGTTGATCGTGACGAGGTGCTCGTGGCCACCCAGCTAAg               >  1:741715/1‑140 (MQ=255)
                                                                                                              cGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCGTTGATCGGGACAAGGTGCTCGGGGCCACCCAGCTAAg               >  1:416651/1‑140 (MQ=37)
                                                                                                              cGGCTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGCCTTCGATAGTGACCAGCTGCTGGCGGCCACCCAGCTAAg               >  2:473854/1‑140 (MQ=37)
                                                                                                              cGGCTGCGGCAGGAAGCCCATAAAATTGGCACGGACGAAGGTAAATTAAAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAg               <  2:741715/140‑1 (MQ=25)
                                                                                                                 cTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGa            >  2:300908/1‑140 (MQ=255)
                                                                                                                 cTGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGAACAGTTGTTCGCGGACACCCAGCTAAGCGa            >  1:535357/1‑140 (MQ=37)
                                                                                                                  tGCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGaa           <  2:1122593/140‑1 (MQ=255)
                                                                                                                   gCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg          >  2:830653/1‑140 (MQ=255)
                                                                                                                   gCGGCAGGAAGGGCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAg          >  1:165799/1‑140 (MQ=255)
                                                                                                                           aagggCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGGCACCCAGCTAAGCGAAGCCGATCTg  >  2:367432/1‑140 (MQ=255)
                                                                                                                           aagggCATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTGGATCGTGAGCAGGTGCTCGCGGCCAGCCAGGTAAGCGAAGCCGATCTg  >  2:437942/1‑140 (MQ=37)
                                                                                                                               |                                                                                                                                       
GCATGCTGGAGCAGGGCTTTGAACGTGAGATTGTGCTGGCGACAACCCAACTCACTGATGCTGATATTCCGAACTGTCATTAAAAGGAGAGATTAATGACTATTGCTGAACGGCTGCGGCAGGAAGGACATCAAATTGGCTGGCAGGAAGGTAAATTAGAAGGTTTGCATGAACAAGCCATTAAAATTGCTTTGCGCATGCTGGAACAGGGCTTTGATCGTGACCAGGTGCTCGCGGCCACCCAGCTAAGCGAAGCCGATCTG  >  NC_000913/2357721‑2357983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: