Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,902 C→G 6.8% intergenic (‑133/+19) abrB ← / ← ybgO regulator of aidB expression; inner membrane protein/putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,9020CG 6.8% 465.8 / 19.2 144intergenic (‑133/+19)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base C (70/64);  new base G (6/4);  total (76/68)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

Only 100 of 148 total aligned reads displayed.

CTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGT  >  NC_000913/747755‑748051
                                                                                                                                                   |                                                                                                                                                     
  gCAAAACTGGCATCCTTCGCTCCCTCATAATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAAGAAAAGGGGAAAATCAGCAAATTTCTGAAAGAGATGCCCTTTACGGCGGCAAAg                                                                                                                                                   >  2:974933/1‑150 (MQ=255)
   cAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                  >  1:619446/1‑150 (MQ=255)
   cAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAgg                                                                                                                                                  >  1:1853869/1‑150 (MQ=255)
      aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGTGCCAGGAAAAAATAAGAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGct                                                                                                                                               <  1:1083107/150‑2 (MQ=255)
       aCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcat                                                                                                                                              <  2:718150/150‑1 (MQ=255)
        cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatc                                                                                                                                             <  1:674274/150‑1 (MQ=255)
        cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatc                                                                                                                                             >  1:50571/1‑150 (MQ=255)
        cTGGCATCCTTCGCTCCCTCATTATTTTCACTTTGGTGTATCTATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAATGAAAATGGTAAAATCAGCTAATTTCAGAAAGAGATGCCCTTACCGGCGGCAAAGCTcatc                                                                                                                                             >  2:387302/1‑150 (MQ=255)
        cTGGCATCCTACGCTCACTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCAGCAAAGGGcatc                                                                                                                                             <  2:285697/150‑1 (MQ=255)
         tGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATATAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAATAGTTGAAAATCAGCAATTTTCTTAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                            <  1:1362119/150‑1 (MQ=255)
         tGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAACGCTGGGTCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                            <  1:262210/150‑1 (MQ=255)
         tGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGATAGGGcatct                                                                                                                                            <  1:1624827/150‑2 (MQ=255)
         tGGCATCCTTCGCGCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGACAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatca                                                                                                                                            <  1:72259/150‑1 (MQ=255)
             aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGt                                                                                                                                        >  2:2064552/1‑150 (MQ=255)
                cTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTGTCCGGCGGCAAAGGGCATCATGGTaaa                                                                                                                                     >  2:585883/1‑150 (MQ=255)
                  tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                                   >  1:57086/1‑150 (MQ=255)
                  tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                                   <  2:238437/150‑1 (MQ=255)
                  tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGCCGGCAAAGGGCATCATGGTAAATc                                                                                                                                   >  2:478629/1‑150 (MQ=255)
                    gCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAATGAGATGCCCTATCCGGCCGCAAAGGGCATCATGGTAAATCAg                                                                                                                                 >  2:1124210/1‑150 (MQ=255)
                    gCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAg                                                                                                                                 <  1:2058892/150‑1 (MQ=255)
                    gCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAg                                                                                                                                 <  2:923793/150‑1 (MQ=255)
                       cccTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGCCGGCAAAGCGCATCATGGTAAATCAGTta                                                                                                                              >  2:1750548/1‑150 (MQ=255)
                         cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                            <  2:2184504/150‑1 (MQ=255)
                         cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                            <  2:1089539/150‑1 (MQ=255)
                               attTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAAt                                                                                                                      >  2:810017/1‑150 (MQ=255)
                                  ttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCtt                                                                                                                   >  1:2043273/1‑150 (MQ=255)
                                      ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                               >  2:788925/1‑150 (MQ=255)
                                         tGGTGTATCGATACATCATAAAGAAAGCTGGGCAATGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGAAATTTTCTGAAAGAGATGGCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                            <  2:590445/150‑1 (MQ=255)
                                          ggTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                           >  1:770692/1‑150 (MQ=255)
                                          ggTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                           <  2:458903/150‑1 (MQ=255)
                                                tCGATACAGTATAAAGATATCTTGGCTAGTAAAAATTATCAGCAATCATTTATATGGTGAAAATAAGAAATTTTCTGAAAGAGATGCCCTTTCCGGTGGCAAAGGGCATCATGTTAAATCAGTTATAGATAATCTTAAAAAGCACCACTg                                                                                                     <  1:397960/150‑1 (MQ=255)
                                                tCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTg                                                                                                     <  1:2083892/150‑1 (MQ=255)
                                                 cGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGTGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGt                                                                                                    <  2:1200369/150‑1 (MQ=255)
                                                    tACAGCATAAAGAAAGCTGGTCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATc                                                                                                 >  2:1316298/1‑150 (MQ=255)
                                                    tACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATc                                                                                                 >  2:1275998/1‑150 (MQ=255)
                                                        gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGACAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                             >  1:328084/1‑150 (MQ=255)
                                                        gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                             <  1:1613087/150‑1 (MQ=255)
                                                         cATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTAAATCTTAAAAAGCACCACTGAATCGAAgg                                                                                            >  1:58286/1‑150 (MQ=255)
                                                         cATAAAGAAAGCAGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATAAAGGTAAATCAGTTATAGATAATCTTAAAAAGCACCACTGTATCGAAgg                                                                                            >  1:246783/1‑150 (MQ=255)
                                                          aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAggg                                                                                           >  2:993031/1‑150 (MQ=255)
                                                           tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAgggg                                                                                          >  1:601584/1‑150 (MQ=255)
                                                              agaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGccc                                                                                       >  2:1629169/1‑149 (MQ=255)
                                                                       ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                              <  1:296489/150‑1 (MQ=255)
                                                                        gCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTATCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                             >  1:52154/1‑150 (MQ=255)
                                                                          cAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                           >  2:83740/1‑150 (MQ=255)
                                                                           aGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATggggg                                                                          <  2:1173913/150‑1 (MQ=255)
                                                                            ggAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGt                                                                         >  1:1496539/1‑150 (MQ=255)
                                                                             gAAAAAATACCAGCAATCTAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATTGTAAATCAGTTATAGTTTATCTTAAATAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTg                                                                        <  2:1952687/150‑1 (MQ=255)
                                                                               aaaaaaaaCAGCAATCAAGGAAAAGCGGAAAATCAGCAATTATCTGAAAGAGATGCCCTTTCCGCCGGCAAAGGGCATCATGGTAAATCAGTAATAGTTAATATTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                      >  1:2127617/1‑150 (MQ=255)
                                                                                  aaaaaCAGCAATCTATGAAAAGCGGAAAATCAGCAATTTTCAGAATTAGATGCCCTTTCCGTCGGGAAAGGGTATCATGGTTAATCAGTTATAGCTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTGGATGGGGGTGCCTgg                                                                   <  1:387302/147‑1 (MQ=255)
                                                                                     accAGCAATCACGGACACGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAAt                                                                <  1:771155/148‑1 (MQ=255)
                                                                                     aaCAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAAt                                                                <  2:523725/150‑1 (MQ=255)
                                                                                       cAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                              >  1:724063/1‑150 (MQ=255)
                                                                                       cAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                              >  2:70201/1‑150 (MQ=255)
                                                                                        aGCAATCAAGGAAAAGGGGACAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                             <  1:1106981/150‑1 (MQ=255)
                                                                                          cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACg                                                           <  2:16803/150‑1 (MQ=255)
                                                                                            aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCAGCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGACAGCCTAGATTGGGGGGCCTGGCATGTACGaa                                                         >  2:691878/1‑150 (MQ=255)
                                                                                             tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGCCGaaa                                                        >  1:1129747/1‑150 (MQ=255)
                                                                                               aaGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGt                                                      >  1:838987/1‑150 (MQ=255)
                                                                                                 ggAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGtta                                                    >  1:779260/1‑149 (MQ=255)
                                                                                                  gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCa                                                   >  2:1408502/1‑150 (MQ=255)
                                                                                                      aGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCg                                               <  1:1370662/150‑1 (MQ=255)
                                                                                                       ggggAAAATCAGCAATTTTATGAAAGAGATGCCCTTTCCTGCGGCAAAGGGCATCATGGTAAATCAGATATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCAGt                                              >  2:653193/1‑150 (MQ=255)
                                                                                                        gggAAAATCAGCGATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGta                                             >  1:681206/1‑150 (MQ=255)
                                                                                                        gggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGta                                             <  1:782935/150‑1 (MQ=255)
                                                                                                          gAACATCCGCAATATTCTGAAAGAGCCTCCCATTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTAAAGTTAATCTTAAAAAGCACCACTGTATCTAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                           <  2:2081208/150‑1 (MQ=255)
                                                                                                          gAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                           >  2:252046/1‑150 (MQ=255)
                                                                                                           aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAg                                          >  2:1646012/1‑150 (MQ=255)
                                                                                                                 aGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttttt                                    >  2:56921/1‑150 (MQ=255)
                                                                                                                  gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGAAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGtttttt                                   >  2:897614/1‑150 (MQ=255)
                                                                                                                  gCAATTTTCTGAAAGAGATGCCCTATCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGtttttt                                   >  2:1687571/1‑150 (MQ=255)
                                                                                                                  ccacatttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGCGCATCATGGTAAATCACTTATAGTTAATCTTAAAACGCCCCACTGTATCGACGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGcttttt                                   <  2:210722/145‑1 (MQ=255)
                                                                                                                    aaTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCg                                 <  2:1863650/150‑1 (MQ=255)
                                                                                                                    aaTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCg                                 <  1:1275998/150‑1 (MQ=255)
                                                                                                                      ttttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                               <  1:2064552/150‑1 (MQ=255)
                                                                                                                       tttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAg                              >  1:981004/1‑150 (MQ=255)
                                                                                                                            gAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCaa                         >  2:966111/1‑150 (MQ=255)
                                                                                                                             aaaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCCCCACTGTATCGAAGGGGCCAGCCTCGATGGGGGGGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAAc                        <  2:912390/150‑1 (MQ=255)
                                                                                                                             aaaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAAc                        >  2:1391305/1‑150 (MQ=255)
                                                                                                                             aaaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAAc                        >  1:1572270/1‑150 (MQ=255)
                                                                                                                              ccatgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGGTAATCTTAAAACGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGATAGTTCAGCCGTAGAGGTTTTTTCGTAGACCcaca                       >  1:1177269/5‑150 (MQ=255)
                                                                                                                              aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCCCTGAATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa                       <  2:1075326/150‑1 (MQ=255)
                                                                                                                              aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa                       <  2:77325/150‑1 (MQ=255)
                                                                                                                               catgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCCGCCTCGATGGGGGTGCGTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACAt                      >  1:823210/4‑150 (MQ=255)
                                                                                                                                gagaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg                     >  2:1043356/1‑150 (MQ=255)
                                                                                                                                atgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg                     >  1:1903620/3‑150 (MQ=255)
                                                                                                                                 tgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAATAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGt                    >  1:1427862/2‑150 (MQ=255)
                                                                                                                                  gaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTc                   >  2:658695/1‑150 (MQ=255)
                                                                                                                                   aTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCa                  >  2:1759759/1‑150 (MQ=255)
                                                                                                                                    tGCCCTTTCCGGCGGCAAAGGGCATCATGGTTAATCAGTTTTAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATCTACGAAAGTTCAGCCCTATAGGTTTTTTCGTAGACCAACATGTCAc                 <  1:1264508/150‑1 (MQ=255)
                                                                                                                                     gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCCGCCTCGATGGGGGTGCCTGGAATCGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGtacct                >  1:1400805/1‑146 (MQ=255)
                                                                                                                                     gCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTAATAGTTAATCTTAAAAAGCACCACTGTATCGAAGCGGCCAGCCTCGAGGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACt                <  2:145161/150‑1 (MQ=255)
                                                                                                                                      cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTg               <  2:1599292/150‑1 (MQ=255)
                                                                                                                                      cccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTg               >  1:1167944/1‑150 (MQ=255)
                                                                                                                                         tttCCGGCGGCTAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTg            <  2:1444646/150‑1 (MQ=255)
                                                                                                                                          ttCCGGCGGCAAAGGTCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGa           <  2:601584/150‑1 (MQ=255)
                                                                                                                                             cggcggCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAAt        >  2:1531117/1‑150 (MQ=255)
                                                                                                                                                cggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCTTATACCAACATGTCACTGGTGaattagg     >  1:2138238/1‑145 (MQ=255)
                                                                                                                                                  gcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCg   <  1:2182634/150‑1 (MQ=255)
                                                                                                                                                   cAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGt  >  2:2091518/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAATCGGCGT  >  NC_000913/747755‑748051

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