Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,3490AC39.2% 130.2 / 16.1 81T106P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (36/12);  new base C (0/31);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 2.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATC  >  NC_000913/1376200‑1376487
                                                                                                                                                      |                                                                                                                                          
ccaccGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTAc                                                                                                                                           >  2:212081/1‑151 (MQ=255)
   ccGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCtt                                                                                                                                        >  1:787751/1‑151 (MQ=255)
     gCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTcc                                                                                                                                      <  1:156790/151‑1 (MQ=255)
      ccTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTGCCTTCCg                                                                                                                                     <  2:384954/151‑1 (MQ=255)
           ggTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAc                                                                                                                                >  1:7916/1‑151 (MQ=255)
            gTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGAAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAcc                                                                                                                               >  1:1163609/1‑151 (MQ=255)
                cTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAccgccg                                                                                                                           >  1:970356/1‑151 (MQ=255)
                cTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGCGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAccgccg                                                                                                                           >  1:1203089/1‑151 (MQ=255)
                    gACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATTTTTACCTTCCGCTTACCGCCGATg                                                                                                                        <  1:676151/150‑1 (MQ=255)
                      cGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGAcc                                                                                                                     >  2:1123270/1‑151 (MQ=255)
                       ggatggatGGGCGGCTTTAGGGAGGAGCGTCGTCTTAATGGCTTAAAGCAGAGCGGGCGCATTCGCGAAGGGGTGGGAGAAGTGGGGGGTAAAGGCATCGTG‑TGCCCGGGTGGGTGGGGGATGTTTCCCTGCGGGTTGCCGCCGATGACCt                                                                                                                    <  2:204525/151‑1 (MQ=255)
                          tggatTGGCGTTTTTATCGAGGAGCGTCGTTTTAATGCCTTAAAGCAGATCGGCCGCATTCTCGAAGCGCTGGAGGAGGTGGGGGGTAAGGCCATCTTC‑GTTCCGGCTGCGGGGGGAATGTTTACCTTCCGCTTACCGCCGATGACCTcg                                                                                                                  <  2:805228/150‑1 (MQ=255)
                          tggatTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTcgc                                                                                                                 >  2:836475/1‑151 (MQ=255)
                           ggatTGGCGGTTTTATCGAAGGGCGTCGTTTTAATGGCTAAAAGCAGATCGAGCGCATTCACGAAGGGCGGACAGAGGGGGGCGGTAAGGGCATCGCC‑GTTCGGGCTGGGTGGGGCATGTTTACCTTCCCCTTCCCGCCGATGACCTcg                                                                                                                  <  2:903357/149‑1 (MQ=255)
                                ggCGACTTTATCGAAGAGCGTCGTCTTAATGGTTAAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTa                                                                                                           <  1:1151925/151‑1 (MQ=255)
                                   gACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGcc                                                                                                        >  1:380083/1‑151 (MQ=255)
                                   gACTTTATCGAAGAGCGTCGTCTAAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGc                                                                                                         <  1:387860/150‑1 (MQ=255)
                                      tttATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTgg                                                                                                     >  1:356368/1‑151 (MQ=255)
                                       ttATCGAAGGGCCTCGTCTTAATGGTTTAACGCAGTCCGCCCCCCTCCTCGAGGCCATGCAAAACGTGGCAGGTAAAGGCAGCCTC‑TTCCCGGTTGCGTGGGGATTTTTTACTTCCCGTTTACCGCCGATGACCTCGCCGCGTAGCCt                                                                                                       <  2:164009/148‑1 (MQ=255)
                                       ttATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTgg                                                                                                     >  2:645454/1‑150 (MQ=255)
                                         aTCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACg                                                                                                  >  1:806675/1‑151 (MQ=255)
                                         aTCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACg                                                                                                  >  1:1087477/1‑151 (MQ=255)
                                         aTCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCGGGCAGAGGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGAc                                                                                                   <  1:223626/150‑1 (MQ=255)
                                          ccGAAGGGCGCCTTCTTAAGGGCTAAAAGCAGATGGGGCGCATTCCCGAGGCGTTGGCAGAAGTGGGCGAAAAGGCCTTCGTC‑GTCCCGCCTGCGTGGGGCATTTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGAc                                                                                                   <  2:348353/148‑1 (MQ=255)
                                            gAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCg                                                                                               >  1:648967/1‑151 (MQ=255)
                                            gAAGAGCGTCGTCTTAATGGCTAAAAGCAGATCGAGCGCTTTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGa                                                                                                    <  1:96438/146‑1 (MQ=255)
                                             aaGGGCGGCGTCTTAATGGCTGAAAGAAGATCGCGCGCATTCTGGAAGCGCGGGCAGAAGTGGGCGGTAAGGGCACCGTC‑GTTCCGGCTGCGTGGGGATGGTTTCCCTCCCGCTTCCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGa                                                                                              <  2:1163511/151‑1 (MQ=255)
                                             aaGAGCGTCGTCTTAATGGCTTAAAGAAGATCGAGCGCTTTCTCGAAGGGCGGGCAGAAGTGGGCGGAAAGGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACg                                                                                                  <  2:490732/147‑1 (MQ=255)
                                                 gCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGc                                                                                          >  1:889687/1‑151 (MQ=255)
                                                 gCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGc                                                                                          >  1:1172606/1‑151 (MQ=255)
                                                  cgtcgtcgtAATGGCTTAAAGCAGATCGAGCGAATTCTCGAAGCGCTGGCAGAAGTGGGAGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGATTACCTTCCGCTTACCGCCGATCGCCTCGCCGCGGAGCCTGGACGGCGAcc                                                                                            >  2:437049/1‑148 (MQ=255)
                                                   gtcgtcTTAAGGGCTAAAGGAAGATCGGGCGCTTCCTCGAAGCGCTGGCAGAAGGGGGCGGAAAGGGCATCGTC‑GTCCGGGTTGGGTGGGGCATGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGAcc                                                                                            <  2:897146/147‑1 (MQ=255)
                                                      gtcTTAATGGCTAAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCaaaa                                                                                      <  1:886792/150‑1 (MQ=255)
                                                            aTGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGgtg                                                                                 >  2:1082435/1‑149 (MQ=255)
                                                             tGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGgtgagt                                                                              >  1:653137/1‑151 (MQ=255)
                                                                 ttAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGAtt                                                                          >  2:367235/1‑151 (MQ=255)
                                                                  tAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTc                                                                         >  1:733598/1‑151 (MQ=255)
                                                                   aaaGCAGGCCGCGCGTTTTCGCGAAGCGCGGGCAGAAGTGGGCGGTAAAGACCTCGTC‑GTTCCGGCTGCGGGGGGCATTTTTCCTTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGAt                                                                           <  2:437386/148‑1 (MQ=255)
                                                                   aaaGCAGCTCGAGCGCATTCTCGAAGCGCTGGAAGAAGTGGGCGGTAAGGGCATCGTC‑GTTCCGGCTGCGGGGGGTTTGTTTCCCTTCCCCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGgtgagtga                                                                            <  1:90054/147‑1 (MQ=255)
                                                                   aaaGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGGGGGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGAtt                                                                          <  2:877646/149‑1 (MQ=255)
                                                                    aaGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTccc                                                                       >  1:885073/1‑151 (MQ=255)
                                                                        aGATCGGGCGCATTCTCGGAGCGCTGGCAGAGGGGGGCGTTAAGGGCTTCGTC‑TTTCCGGCTGTGTGGGGATTGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCt                                                                      <  2:196690/148‑1 (MQ=255)
                                                                        aGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGAAAAGGCATCGTC‑GTCCCGGCTGCGTGGGGCATGTTTACCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTgcg                                                                   <  2:180351/151‑1 (MQ=255)
                                                                         gagcgggcgcATTCTCGAGGCGCGGGCAGAAGTGGGCGGTAAAGGCATCTTC‑GTTCGGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTgcg                                                                   <  2:1182075/144‑1 (MQ=255)
                                                                           acggccgtaTTCTCAAAGCGCTGGCTGAAGTGGG‑GGAAAATGGTGGTTCAGTTCGGCAGGCGTGGGAAATTTTTCCCTACGATTTGCCGCCGATGGCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGcgcg                                                                 <  1:1111595/142‑1 (MQ=255)
                                                                             gAGCGCATTCTCGAGGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑TTCCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGt                                                                <  2:69643/149‑1 (MQ=255)
                                                                             gAGCGCATTCCCGAAGCGCTGGAAGAAGTGGGCGGTAAAGGCACCGTC‑GTTCCGGCTGCGTGGGGCATTTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGt                                                                <  1:572949/149‑1 (MQ=255)
                                                                              aGCGCATTCTCGAGGGCTTGGCAGAGTTGGGCGATAAAGGCATCGTC‑GTTCGGGTTGCGGGGGAAGTGTTTCCCTCCCGTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTa                                                               <  2:476209/149‑1 (MQ=255)
                                                                              aGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACt                                                             >  1:570179/1‑151 (MQ=255)
                                                                              aGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTCCGTGGGGCATGTTTACCTTCCGATTACCGCCGACGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGCGAGTGATTCCCTGCGCGTAc                                                              >  1:357252/1‑150 (MQ=255)
                                                                                 acATTCTCAAAGCGTTGGAAGAAGTGGGCGGAAAAGGCATCGTC‑TTCCCGGCTGCGGGGGGCATTTTTTCCTTCCGTTTCCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGa                                                          <  1:400237/150‑1 (MQ=255)
                                                                                     tctcGAAGCGCTGGCAGAGGTGGGCGGTAAAGGCATCGGC‑GTTCCGGCGGGGTGGGGCATTTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACa                                                       <  2:325483/150‑1 (MQ=255)
                                                                                     tctcGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAAc                                                        >  1:1179100/1‑149 (MQ=255)
                                                                                       tcGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAgg                                                     >  1:800642/1‑150 (MQ=255)
                                                                                       tcGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAgg                                                     >  2:215094/1‑150 (MQ=255)
                                                                                         gaggcgcTGGCAGAAGTGGGGGGTAAAGCCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGt                                                    <  1:687988/146‑1 (MQ=255)
                                                                                         gAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTcg                                                  >  1:1063085/1‑151 (MQ=255)
                                                                                                 ggaagaggggggTGGTAAAGGCACGGTC‑TTTCCGGCTGCGTGAGCCATTTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTcgccg                                               <  1:567315/137‑1 (MQ=255)
                                                                                                 ggCAGAAGGGGGCGGAAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGt                                          <  1:907969/151‑1 (MQ=255)
                                                                                                     gAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCg                                      >  1:1138586/1‑151 (MQ=255)
                                                                                                     gAAGTGGGCGGTAAAGGCATCGTC‑GTCCCGGCGGCGGGGGGCATGTTTCCCTTCCGTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTa                                         <  1:564808/148‑1 (MQ=255)
                                                                                                       aGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCgg                                     >  2:265137/1‑150 (MQ=255)
                                                                                                        gTGGGCGGTAAGGGCATCGTC‑GTTCCGGCGGCGTGGGGCATGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGaa                                   <  2:866933/151‑1 (MQ=255)
                                                                                                           ggCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCg                                <  1:902274/151‑1 (MQ=255)
                                                                                                           ggCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCg                                <  1:811920/151‑1 (MQ=255)
                                                                                                            ggggTAAGGGCTTCGTC‑GTCCCGGCGGCGTGGGGCTTTTTTCCCTCCCGGTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCgt                               <  2:112992/149‑1 (MQ=255)
                                                                                                            gCGGTAAAGGCATCGTC‑TTTCCGGCTGGGTGGGGCATTTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCgt                               <  1:229116/151‑1 (MQ=255)
                                                                                                                 aaaGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCgtggtg                           <  2:1087424/150‑1 (MQ=255)
                                                                                                                 aaaGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGgtgt                          >  1:416552/1‑151 (MQ=255)
                                                                                                                  aaGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGCACCGTGgtgt                          >  2:946788/1‑150 (MQ=255)
                                                                                                                  aaGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGgtgt                          >  1:46694/1‑150 (MQ=255)
                                                                                                                     ggcttcgccttccgggCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTAt                        <  1:851021/136‑1 (MQ=255)
                                                                                                                      caccgtc‑gtCCCGGCTGCGTGGGGCATGTTTCCCTTCCCCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTAtct                      <  2:363179/147‑1 (MQ=255)
                                                                                                                        accgtggtccgggCTGGGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTAtct                      <  2:1061549/138‑1 (MQ=255)
                                                                                                                         cgtc‑gtTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGa                   >  2:997686/1‑150 (MQ=255)
                                                                                                                          gtc‑gtTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAAc                 >  1:442123/1‑151 (MQ=255)
                                                                                                                                 ccggtggCGTGGGGTATGTTTCCCTCCCCTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGtt             <  2:1169598/143‑1 (MQ=255)
                                                                                                                                   ggCTGCGTGGGGCATGTTTCCTTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCg                    <  2:1087532/140‑1 (MQ=255)
                                                                                                                                    gCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAAcc        >  2:1253886/1‑151 (MQ=255)
                                                                                                                                           ggggCATTTTTCCCTCCCGCTTCCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTAt   <  2:187010/149‑1 (MQ=255)
                                                                                                                                           ggggAATTTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATc  <  1:1934/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
CCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTC‑GTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATC  >  NC_000913/1376200‑1376487

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