Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA24.6% 146.3 / 19.6 65intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (39/10);  new base A (0/16);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 8.20e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTC  >  NC_000913/752986‑753257
                                                                                                                                            |                                                                                                                                    
gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                                                      <  1:616645/140‑1 (MQ=255)
gggTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCatat                                                                                                                                      <  1:1137530/140‑1 (MQ=255)
   ttATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATaaa                                                                                                                                   >  2:1484944/1‑140 (MQ=255)
           gcgcttACAATGAATTTTTTAATTGAAATTTTTGGTTGGTTAATAATTTGAGCTTTTTTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTTTTTTTAAATTTTTTTGTGTTTAGGCATAAAAAAACCAccc                                                                                                                           <  2:698900/135‑1 (MQ=255)
           gCGACTACCATGAATTTTTTATTTAAAAGTATTGGGTTGCTGATAATTTGAGTTGTTCTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGTTTTTATTTTTAAATTTTTTTGTGTTTAGGAATAAAAAAATCAAcc                                                                                                                           <  2:824612/140‑1 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGTGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                           >  1:507759/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                           >  2:1460056/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                           >  1:698900/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTAAAAGTATTGGGTTGCTGATAATTTGAGCTTTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                           <  1:176379/140‑1 (MQ=255)
            cGACTACCATGAAGTTTTTAATAAAAAGTATGGGGTTGCTGATAATTTGAGTTTTTCTTTTTTTTTTAAATATATATATAGGTTTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAATCAAccc                                                                                                                          <  2:1225011/140‑1 (MQ=255)
             gACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCg                                                                                                                         >  2:604767/1‑140 (MQ=255)
                  ccATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTGTGTTTAGGCATATAAAAATCAACCCGCCatat                                                                                                                   >  1:493545/1‑140 (MQ=255)
                         ttttttAATTGAAATTATTGGGTTGTTAATAATTTGGGTTGTTTTTTTTTTTTTAAATATCTATATAGGTCTGTAAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAAcgg                                                                                                             <  2:1086017/139‑1 (MQ=255)
                             ttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCggg                                                                                                         >  1:569667/1‑140 (MQ=255)
                                  taaaaGTATTGTGTTGCTGATAATTTGAGCTTTTCTTTTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTaaa                                                                                                    <  2:507759/138‑1 (MQ=255)
                                    aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                                  >  1:1376351/1‑140 (MQ=255)
                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                       >  2:1613043/1‑140 (MQ=255)
                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                       >  2:1121148/1‑140 (MQ=255)
                                               ttGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACtt                                                                                       >  1:1580776/1‑140 (MQ=255)
                                                 gCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAg                                                                                     <  1:1121148/140‑1 (MQ=255)
                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                    >  1:1205549/1‑140 (MQ=255)
                                                  cTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGc                                                                                    >  1:201005/1‑140 (MQ=255)
                                                       attttGAGTTTTTTTTTTCTTTTTAAATATCTATATAGGTTTTTAAATGTTTTTTTTTTTTACAATTTTTGTTTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                              <  2:493545/138‑1 (MQ=255)
                                                        ttttGAGTTTTTCTATTCTTTTTAAATATCTATATAGGTCTTTTAATGGTTTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                              <  1:23212/139‑1 (MQ=255)
                                                        ttttGAGTTGTTCTATTTTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCaa                                                                              <  1:1010564/139‑1 (MQ=255)
                                                            gAGCTTTTTTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGCATATAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                          <  2:228711/140‑1 (MQ=255)
                                                            gAGCTTTTCTATTTTTTTTAAATATTTATATGGGTTTGTTAATGGATTTTATTTTTACAATTTTTTGTTTTTAGGAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                          <  2:1378352/140‑1 (MQ=255)
                                                            gAGCTGTTCTTTTCTTTTTAAATATTTATATAGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                          <  2:923939/140‑1 (MQ=255)
                                                              gCTGTTTTTTTTTTTTTAAATATCTATATAGGTTTTTAAATGGATTTTTTTTTTAAATTTTTTTTTTTTAAGTAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                        <  2:201005/140‑1 (MQ=255)
                                                              gCTGTTCTATTTTTTTTAAATCTCTATATAGGTTTGTTAATGGTTTTTATTTTTACAATTTTTTGTTTTTAGGCATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTa                                                                        <  1:1460056/140‑1 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                          >  1:1452376/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                          >  2:1193927/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                          >  2:552651/1‑140 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGGTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATCTAAAAATCAACCCGCCATATGAACGGCGGGTTTAAAAATTTACAACTTAGCAATCAACCATTAACGCTTGGTATAGc                                                          >  2:607853/1‑140 (MQ=255)
                                                                            tttAAATATATATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                          >  1:1451900/1‑140 (MQ=255)
                                                                              tAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCtt                                                        <  1:1095149/140‑1 (MQ=255)
                                                                                          tgggTCTTTTAATGGTTTTTATTTTTACTTTTTTTTGTTTTTAGAAATAAAAAAACCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                            <  2:1022837/138‑1 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                         >  2:214473/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGtttt                                         >  2:1373238/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                            <  1:1282162/137‑1 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                            >  2:1282162/1‑137 (MQ=255)
                                                                                               ttGTTTATGGATTTTATTTTTACAATTTTTTGTTTTTAGTTATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                       <  2:569667/139‑1 (MQ=255)
                                                                                               ttGTTAATGGATTTTTTTTTTACATTTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                       <  1:1194531/139‑1 (MQ=255)
                                                                                               ttGTTAATGGATTTTTTTTTTAAATTTTTTTGTTTTTAGGAAAAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                       <  2:714126/139‑1 (MQ=255)
                                                                                               ttGTTAATGGATTTTATTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                       <  1:1012830/139‑1 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATTTTTACAACTTAGCAATCAACCATTAAAGCTTGATATCGCCTTTTAAGTCGCGTTTTTCa                                      >  2:1061483/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTct                                      >  1:1022837/1‑139 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                      >  1:150686/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                      >  2:1039242/1‑140 (MQ=255)
                                                                                                tGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGGTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCa                                      >  2:777482/1‑140 (MQ=255)
                                                                                                 tttAATGGTTTTTTTTTTTAAAATTTTTTGTGTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                     <  1:604767/139‑1 (MQ=255)
                                                                                                 gtaaaTGTATTTTTTTTTTAAAATTTTTTGTTTTTAGGAAAAAAAAAAACACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                     <  2:1562984/137‑1 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTTAAGTCGCGTTTTTCata                                    >  1:710483/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                            >  2:309193/1‑132 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                            <  1:309193/132‑1 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                    >  1:1030826/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                    >  2:709270/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                    >  1:115649/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                    >  2:1068597/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGGCGCGTTTTTCata                                    >  1:487127/1‑140 (MQ=255)
                                                                                                    attagcttttattttTCAATTTTTTAAAGTTTTGATATACAAAAAAAATCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAGCCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                  <  2:864020/134‑1 (MQ=255)
                                                                                                    aatgttttttttttttAAAATTTTTTTTTTTTAGAAATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                  <  2:487127/134‑1 (MQ=255)
                                                                                                    aatgttttttttttttAAAATTTTTTGTTTTTAGGAATAAAAAAAACAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATc                                  <  2:1205549/134‑1 (MQ=255)
                                                                                                       ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAGCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                               >  2:385243/1‑140 (MQ=255)
                                                                                                       ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTTAACTATTTACAACTTAGCAATCAACCATTAACGCTTGATATTGCTTTTAAAGTCCCGTTTTTTATATCCTg                               >  2:39600/1‑140 (MQ=255)
                                                                                                       ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                               >  2:126583/1‑140 (MQ=255)
                                                                                                       ggATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTg                               >  1:714126/1‑140 (MQ=255)
                                                                                                           tttattttTTAAATTTTTTTTTTGAAGGAAAAAAAAATATAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGtata                           <  2:913255/140‑1 (MQ=255)
                                                                                                                      aaTTTTTTGTGTTTAGGCATATAAAAATCAGCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGAcgcg                <  1:385243/140‑1 (MQ=255)
                                                                                                                           tttGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACgg           <  1:709270/140‑1 (MQ=255)
                                                                                                                             tgttTTTAGGCATAAAAAAATCCACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc         <  2:710483/140‑1 (MQ=255)
                                                                                                                             tgttTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc         <  1:894309/140‑1 (MQ=255)
                                                                                                                             tgttTTTAGACATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc         <  2:1451900/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc         <  1:1193927/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc         <  1:1039242/140‑1 (MQ=255)
                                                                                                                             tgtgTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGc         >  2:1426173/1‑140 (MQ=255)
                                                                                                                                     aGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTc  >  1:394519/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                    
GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTC  >  NC_000913/752986‑753257

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