Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,897 G→C 6.3% intergenic (‑128/+24) 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,8970GC 6.3% 424.4 / 19.1 143intergenic (‑128/+24)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (67/64);  new base C (6/3);  total (75/68)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-02

Only 100 of 151 total aligned reads displayed.

CCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAAT  >  NC_000913/747750‑748045
                                                                                                                                                   |                                                                                                                                                    
ccccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAACGAAAGCTGGGCCAGGAAAAAATAACAGCTATGAAGGAAAAGGGGAAAATCAGCAATTTTATGAAAGAGATGCCCTTTCcggc                                                                                                                                                    >  1:375785/1‑150 (MQ=255)
ccccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATTCCCTTTCcggc                                                                                                                                                    >  2:643065/1‑150 (MQ=255)
ccccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTGTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCcggc                                                                                                                                                    >  1:1379892/1‑150 (MQ=255)
 cccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCcggcg                                                                                                                                                   <  1:657062/150‑1 (MQ=255)
  ccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTGTCTCTAAGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCATTTCcggcgc                                                                                                                                                  >  1:1104118/1‑149 (MQ=255)
     cTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTATCCGGCGGCaa                                                                                                                                               >  2:1045646/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)
             cTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAATTGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGAAAGGGcatc                                                                                                                                       <  2:1266335/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)
              tGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTGCCGCCGGATAGGGcatct                                                                                                                                      <  1:1624827/150‑2 (MQ=255)
               ggCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGcatcat                                                                                                                                     <  2:451925/150‑1 (MQ=255)
                gCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATg                                                                                                                                    <  2:1764541/150‑1 (MQ=255)
                  aTCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGt                                                                                                                                  >  2:2064552/1‑150 (MQ=255)
                       tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                             <  2:238437/150‑1 (MQ=255)
                       tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATc                                                                                                                             >  1:57086/1‑150 (MQ=255)
                       tCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGCCGGCAAAGGGCATCATGGTAAATc                                                                                                                             >  2:478629/1‑150 (MQ=255)
                         gCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAATGAGATGCCCTATCCGGCCGCAAAGGGCATCATGGTAAATCAg                                                                                                                           >  2:1124210/1‑150 (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)
                                     ttttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                               <  2:1675044/150‑1 (MQ=255)
                                     ttttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                               <  2:2085100/150‑1 (MQ=255)
                                       ttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCtt                                                                                                             >  1:2043273/1‑150 (MQ=255)
                                       ttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCtt                                                                                                             >  1:256149/1‑150 (MQ=255)
                                           ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                         >  2:1370662/1‑150 (MQ=255)
                                           ctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaa                                                                                                         >  2:788925/1‑150 (MQ=255)
                                              tGGTGTATCGATACATCATAAAGAAAGCTGGGCAATGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGAAATTTTCTGAAAGAGATGGCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                      <  2:590445/150‑1 (MQ=255)
                                               ggTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                     >  1:979341/1‑150 (MQ=255)
                                                 tgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcacc                                                                                                   <  2:1602270/150‑1 (MQ=255)
                                                     tCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTg                                                                                               <  1:810017/150‑1 (MQ=255)
                                                      cGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGCAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCACTTATAGTTAATCTTAAAAAGCACCACtct                                                                                              >  2:1743873/1‑148 (MQ=255)
                                                       gATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTa                                                                                             >  2:59370/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)
                                                             gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAg                                                                                       >  1:453332/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)
                                                                  aagacagCTGGGCCAGCAAAAAATAACACCAATCAAGGAAAACGCGAAAACCAGCAATTTTCTGAAATAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGcc                                                                                  <  1:2089909/150‑1 (MQ=255)
                                                                   agaaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGccc                                                                                 >  2:1629169/1‑149 (MQ=255)
                                                                       agCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGCGGAAACTCAGCAATTTTCTGAAAGAGATGCCCTTTCCTGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCt                                                                             <  1:557003/150‑1 (MQ=255)
                                                                            ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                        <  1:296489/150‑1 (MQ=255)
                                                                            ggCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATg                                                                        >  1:1319757/1‑150 (MQ=255)
                                                                             gCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATTCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCTCCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                       >  2:14193/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)
                                                                                  taaaaaaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTg                                                                  <  2:1961228/149‑1 (MQ=255)
                                                                                    aaaaaaaaCAGCAATCAAGGAAAAGCGGAAAATCAGCAATTATCTGAAAGAGATGCCCTTTCCGCCGGCAAAGGGCATCATGGTAAATCAGTAATAGTTAATATTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGcc                                                                >  1:2127617/1‑150 (MQ=255)
                                                                                      aaaTAACAGCAATCAAGGAAAATGGGAAAATCAGCAATATTCTGAATGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTATAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                              <  1:701655/150‑1 (MQ=255)
                                                                                       aaaaaCAGCAATCTATGAAAAGCGGAAAATCAGCAATTTTCAGAATTAGATGCCCTTTCCGTCGGGAAAGGGTATCATGGTTAATCAGTTATAGCTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTGGATGGGGGTGCCTgg                                                             <  1:387302/147‑1 (MQ=255)
                                                                                         tAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCATGGaa                                                           >  2:1520006/1‑150 (MQ=255)
                                                                                          accAGCAATCACGGACACGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAAt                                                          <  1:771155/148‑1 (MQ=255)
                                                                                           aCAGCAATCAAGGAAAGGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATg                                                         <  1:2005099/150‑1 (MQ=255)
                                                                                            cAGCAATCAAGGAACAGCGTTAAATCAGCAATTTTCTGACAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGTGGTTCCTGGAATgg                                                        <  1:114300/150‑1 (MQ=255)
                                                                                            cAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                        >  1:724063/1‑150 (MQ=255)
                                                                                            cAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATgg                                                        >  2:70201/1‑150 (MQ=255)
                                                                                             aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                       <  1:59370/150‑1 (MQ=255)
                                                                                               cAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACg                                                     <  2:16803/150‑1 (MQ=255)
                                                                                                 aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCTGCGGCAAAGGGCATCATGGTAAATAAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGTGGTGCCTGGAATGGCCGaa                                                   >  1:122251/1‑150 (MQ=255)
                                                                                                 aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCAGCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGACAGCCTAGATTGGGGGGCCTGGCATGTACGaa                                                   >  2:691878/1‑150 (MQ=255)
                                                                                                  tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGCCGaaa                                                  >  1:1129747/1‑150 (MQ=255)
                                                                                                    aaGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGt                                                >  1:838987/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)
                                                                                                               gAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGtata                                     >  2:252046/1‑150 (MQ=255)
                                                                                                                      aGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGttttt                              >  2:56921/1‑150 (MQ=255)
                                                                                                                       gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGtttttt                             >  1:774179/1‑150 (MQ=255)
                                                                                                                       gCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGAAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGtttttt                             >  2:897614/1‑150 (MQ=255)
                                                                                                                       ccacatttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGCGCATCATGGTAAATCACTTATAGTTAATCTTAAAACGCCCCACTGTATCGACGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGcttttt                             <  2:210722/145‑1 (MQ=255)
                                                                                                                         aaTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCg                           <  1:1275998/150‑1 (MQ=255)
                                                                                                                         aaTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCg                           <  2:57086/150‑1 (MQ=255)
                                                                                                                         aaTTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCg                           <  2:1863650/150‑1 (MQ=255)
                                                                                                                          aTTTTCTGAAAGAGATGCCCTTTCAGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGt                          <  1:643065/150‑1 (MQ=255)
                                                                                                                           ttttCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                         <  1:2064552/150‑1 (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)
                                                                                                                                   aaGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa                 <  2:77325/150‑1 (MQ=255)
                                                                                                                                     gagaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg               >  2:1043356/1‑150 (MQ=255)
                                                                                                                                     atgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg               >  1:1903620/3‑150 (MQ=255)
                                                                                                                                     atgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATg               >  1:2061713/3‑150 (MQ=255)
                                                                                                                                      tgaTGCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAATAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGt              >  1:1427862/2‑150 (MQ=255)
                                                                                                                                       gaTGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTc             >  2:658695/1‑150 (MQ=255)
                                                                                                                                          gCCCTTTGCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCCGCCTCGATGGGGGTGCCTGGAATCGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGtacct          >  1:1400805/1‑146 (MQ=255)
                                                                                                                                            ccTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg        >  1:747380/1‑150 (MQ=255)
                                                                                                                                              tttCCGGCGGCTAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTg      <  2:1444646/150‑1 (MQ=255)
                                                                                                                                               ttCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGa     <  2:50571/150‑1 (MQ=255)
                                                                                                                                                  cggcggCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAAt  >  2:1531117/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTTCCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTGAAAT  >  NC_000913/747750‑748045

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: