Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA45.3% 76.1 / 49.5 75intergenic (‑330/+60)ybgD/gltAputative fimbrial‑like adhesin protein/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (30/11);  new base A (0/34);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 5.47e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCA  >  NC_000913/752977‑753273
                                                                                                                                                       |                                                                                                                                                    
aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCAtata                                                                                                                                                     >  2:213164/1‑150 (MQ=255)
aTTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATaa                                                                                                                                                    <  1:781794/151‑1 (MQ=255)
 ttttCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATaa                                                                                                                                                    >  2:758834/1‑150 (MQ=255)
     cATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAAt                                                                                                                                                <  1:669735/150‑1 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATc                                                                                                                                               >  2:313451/1‑150 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATc                                                                                                                                               >  1:456756/1‑150 (MQ=255)
      aTAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCa                                                                                                                                              <  2:267467/151‑1 (MQ=255)
       tAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCa                                                                                                                                              >  2:774988/1‑150 (MQ=255)
         gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAAcc                                                                                                                                           >  1:327806/1‑151 (MQ=255)
                     cGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAAc                                                                                                                               >  2:831332/1‑151 (MQ=255)
                         tACCATGAAGTTTTTAATTGAAA‑GTATTGGGATGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcg                                                                                                                           >  2:1035623/1‑151 (MQ=255)
                               gAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTa                                                                                                                      >  2:1025848/1‑150 (MQ=255)
                                aattttttAATTGAAA‑GTTTTGGGTTGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTA‑CAATTTTTTG‑TTTTAAGGAATATAAAAATCACCCCGCCATATGAACGGCGGGTTaaa                                                                                                                    <  2:771079/148‑1 (MQ=255)
                                aaGTTTTAAATTGAAA‑GAATGGGTTGGCTAATTGTTTGGGTTGTTCTTTTTTTTTAAAATATCTAAATAGGTCTTTAAATGGATTTTATTTTGA‑AAATTTTTGG‑TGTTTAGGCAAAAAAAAAACAACCCGCCATATGAACGGCGGGTTaag                                                                                                                    <  1:402045/151‑2 (MQ=255)
                                   tttttAATTAAAA‑GTATTGGGTTGCTGATATTTTGAGCTGTTCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTaaaa                                                                                                                   <  2:84850/149‑1 (MQ=255)
                                    ttttAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTaaaa                                                                                                                   >  1:1033388/1‑148 (MQ=255)
                                    tttaaatttaaatGTATTGGTTTAGTTATAGTTTGCGTTGTTTTATTCTTTTAAAAAATTTATATGGGTAT‑TTAAAGGTTTTTTTTTTTACCATTTTTTTT‑TTTTTGGGCTTATAAAAAACACCCCGCCATATGAACGGCGGGTTAAAat                                                                                                                  <  2:357760/150‑1 (MQ=255)
                                       tAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATAttt                                                                                                              >  1:131730/1‑150 (MQ=255)
                                                   aTTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTAAAATTAga                                                                                                    >  1:186001/1‑148 (MQ=255)
                                                     tggtttcttGATAATTTGGGTTTTTTTTTTTTTTTAAAATATTTTTATAGTTCTGTTAATGGTTTTTTTTTTTA‑CAATTTTTTG‑TTTTTGGGCATATAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                 <  2:427646/142‑1 (MQ=255)
                                                     tGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                 >  1:118117/1‑150 (MQ=255)
                                                      gggTTGTTGATATTTTGAGCTTTTCTTTTTTTTTAAAATATCTATATAGGTTTTTAAATGGTTTTTTTTTTTA‑CATTTTTTTG‑TTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAAt                                                                                                 <  2:971857/149‑1 (MQ=255)
                                                      gggTTGCTAATAATTGGGGTTTTTTTATTTTTTTTAAATATAAAAATAGGTCTTTTATTGATTTTTTTTTTTA‑CATTTTTTTG‑TTTTTAGACAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATc                                                                                                <  1:403730/150‑1 (MQ=255)
                                                         ttGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTTTTTTAATGGATTTTTTTTTTA‑CATTTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                                              <  2:739375/149‑1 (MQ=255)
                                                               aTAATTTGAGCTTTTTTATTTTTTTAAAATACTTATATGGTTTTTTAAATGATTTTTATTTTAA‑TATTTTTTTT‑TGTTAAGGCAAATAAAAATAACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                                        <  2:930654/149‑1 (MQ=255)
                                                                tAATTTGAGCTTTTCTTTTTTTTTAAAATATTTTTATGGGTCTGTTAAGGTATTTTATTTTTA‑AAATTTTTTG‑TTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                                       <  2:301790/149‑1 (MQ=255)
                                                                 attttGAGCTTTTTTATTTTTTTAAAATATCTAAATAGGTCTGTTAATGAATTTTTTTTTTA‑CATTTTTTTG‑TTTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTaa                                                                                       <  2:649721/146‑1 (MQ=255)
                                                                 aaTTTGAGTTGTTCTTTTTTTTTAAAATATCTATATAGTTCTGTTAATGTATTTTATTTTTA‑CAATTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                     <  2:198460/150‑1 (MQ=255)
                                                                   tttGTGCTTTTTTTTTTTTTTTAAATATTTAAATGGTTTTGTTAATGAATTTTATTTTTA‑AATTTTTTTG‑TTTTAAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                                          <  1:28513/143‑1 (MQ=255)
                                                                   tttGAGTTTTTTTTTTTTTTTAAAATTTTTATATAGTTTTTTTAATGGATTTTATTTTTA‑CATTTTTTTG‑TTTATAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAAc                                                                                      <  1:930371/147‑1 (MQ=255)
                                                                    ttggtttttttttttCTTTTAAAATTTTTATATATTTCTTTTAATGTTTTTTTTTTTTA‑CAATTTTTTG‑TTTGTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                                     <  1:680405/134‑1 (MQ=255)
                                                                      gagtttttCTATTTTTTTTAAATATTTATATAGTTTTTTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGAATATAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCtt                                                                                  <  2:906777/142‑1 (MQ=255)
                                                                          ttttttttttttttaAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTA‑AAATTTTTTG‑TGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGata                                                                              <  2:862263/138‑1 (MQ=255)
                                                                          tGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTTGCAAACAACAATTAACACTTGATCTCg                                                                           >  1:280681/1‑151 (MQ=255)
                                                                            tttttttttttttAAATATCTATATAGTTTTTTTAATGGATTTTATTTTTA‑AAATTTTTTG‑TGTTTAGGAAAAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                           <  1:911044/141‑1 (MQ=255)
                                                                            ttttttttttttaAAATATTTATATGGGTTTTTTAATGGATTTTTTTTTTA‑CAATTTTTTG‑TTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                                          <  2:665742/142‑1 (MQ=255)
                                                                            ttCTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCg                                                                           <  2:43424/149‑1 (MQ=255)
                                                                             tCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtt                                                                        >  2:1025076/1‑151 (MQ=255)
                                                                                 ttCTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAGTCAACCCGCCATATGAACGGCGGGGTTAAATTTTTTCAACTTAAAAATAAAACAATAACGCTTGAAATCTCttt                                                                       >  1:605357/1‑148 (MQ=255)
                                                                                  tttttttAAATATTTATATAGTTTTTTAAATGTTTTTTTTTTTTA‑CATTTTTTTT‑TTTTTAGAAAAAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                                    <  1:831288/148‑1 (MQ=255)
                                                                                  tCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTaa                                                                    >  2:1306/1‑150 (MQ=255)
                                                                                    tttttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCTCTTTTaaa                                                                   >  1:66696/1‑149 (MQ=255)
                                                                                       taaaaTATTTATATAGGTCTTTTAATGTATTTTTTTTTTA‑AATTTTTTTG‑TTTTTAGGCATAAAAAAATAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTc                                                                <  1:807044/147‑1 (MQ=255)
                                                                                          attatCTATATAGGTTTGTTAATGTATTTTTTTTTTA‑CAATTTTTTG‑TTTTTAGATATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTcg                                                               <  1:831546/145‑1 (MQ=255)
                                                                                            tatCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAAACATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                         >  2:83319/1‑151 (MQ=255)
                                                                                                 atataGGTTTTTTAATGTTTTTTTTTTTTA‑CATTTTTTTG‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                                         <  1:427330/146‑1 (MQ=255)
                                                                                                   ataGGTTTGTTAATGTTTTTTATTTTTA‑AATTTTTTTG‑TTTTTAGACTAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                                    <  1:112593/149‑1 (MQ=255)
                                                                                                    taGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTATAAATTAACAATCAACCATTAAAGCCTGATATCGCTTTTAAAGTCGCGtttttta                                                      >  1:118125/1‑144 (MQ=255)
                                                                                                     aGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTTCAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATcc                                                 >  1:925556/1‑150 (MQ=255)
                                                                                                      ggTTTGTTAATGGATTTTATTTTTA‑AATTTTTTTT‑TGTTTAGGCAAATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGttttt                                                        <  2:27957/142‑1 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCATTAACGCTTGGTATCGCTTTTTAAGTCGCGTTTTTCATATCCTGt                                              >  1:908381/1‑151 (MQ=255)
                                                                                                       gTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTCAAATTTTTTCAACTTAGCCATCAACCATTAACGCTTGGTATCGCTTTTAAAGTCGCGTTTTTCATATCCt                                                >  2:873310/1‑149 (MQ=255)
                                                                                                        tttGTTAATGGTTTTTTTTTTTA‑CAATTTTTTG‑TGTTTAGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGt                                              <  2:626371/148‑1 (MQ=255)
                                                                                                          ttttAATGGTTTTTATTTTTA‑CAATTTTTTT‑TTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtat                                            <  2:78944/148‑1 (MQ=255)
                                                                                                          tGTTAAGGGTTTTTTTTTTTA‑AATTTTTTTG‑TGTTTAGGAAAAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGta                                             <  1:1073724/149‑1 (MQ=255)
                                                                                                              aaTGGTTTTTATTTTTA‑AAATTTTTTG‑TGTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                        <  1:774945/150‑1 (MQ=255)
                                                                                                              aaTGGATTTTATTTTTA‑CAATTTTTTG‑TTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAg                                        <  2:999155/150‑1 (MQ=255)
                                                                                                              aaTGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAAATTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTTATATCCTGt                                              >  1:861523/1‑144 (MQ=255)
                                                                                                               aTGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAAATTAGCAATCAACCCTTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGc                                       >  1:404023/1‑150 (MQ=255)
                                                                                                                 ggATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCt                                      >  2:1010333/1‑149 (MQ=255)
                                                                                                                    ttttttttttA‑AATTTTTTTG‑TTTTTAGGTATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcg                                  <  2:728427/150‑1 (MQ=255)
                                                                                                                     tttttttttt‑TAATTTTTTGTTTTTTAGAAAAAAAAAAACAAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGa                                    <  1:742773/148‑1 (MQ=255)
                                                                                                                      ttattttTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCgg                               >  2:521783/1‑151 (MQ=255)
                                                                                                                          tttaA‑CATTTTTTTG‑TTTTTGGGCAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACg                            <  2:517147/150‑1 (MQ=255)
                                                                                                                               aaaaTTTTTTT‑TGTTTAGGAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCa                        <  1:518194/148‑1 (MQ=255)
                                                                                                                                 atttttttt‑tttttAGTAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCaa                       <  2:544624/147‑1 (MQ=255)
                                                                                                                                   ttttttt‑tGTTAAGGAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATc                     <  1:25563/149‑1 (MQ=255)
                                                                                                                                     ttttG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCtt                   >  1:189343/1‑149 (MQ=255)
                                                                                                                                      tttG‑TGTTTAGTAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCa                 <  1:520171/150‑1 (MQ=255)
                                                                                                                                       ttG‑TGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCAt                <  2:458414/150‑1 (MQ=255)
                                                                                                                                           tttttAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATa               <  1:368246/146‑1 (MQ=255)
                                                                                                                                           tgtTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAACAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAcc             >  1:47265/1‑150 (MQ=255)
                                                                                                                                             tttAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGt           <  1:758790/150‑1 (MQ=255)
                                                                                                                                             tttAGGAAAAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCg            <  1:313905/149‑1 (MQ=255)
                                                                                                                                                  gCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGGTTTTCATATCCTGTATACAGCTGa                                    >  2:344056/1‑120 (MQ=255)
                                                                                                                                                   cATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACtgt     >  2:520810/1‑150 (MQ=255)
                                                                                                                                                      ataAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgc                                 <  1:344036/119‑1 (MQ=255)
                                                                                                                                                      ataAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCa  <  1:724996/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
ATTTTCATAGGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAA‑GTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTA‑CAATTTTTTG‑TGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATACCGTCACTGTGCA  >  NC_000913/752977‑753273

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