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

CCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATC  >  NC_000913/1376200‑1376487
                                                                                                                                                      |                                                                                                                                          
ccaccGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTAc                                                                                                                                           >  2:212081/1‑151 (MQ=255)
   ccGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCtt                                                                                                                                        >  1:787751/1‑151 (MQ=255)
     gCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTcc                                                                                                                                      <  1:156790/151‑1 (MQ=255)
      ccTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTGCCTTCCg                                                                                                                                     <  2:384954/151‑1 (MQ=255)
           ggTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAc                                                                                                                                >  1:7916/1‑151 (MQ=255)
            gTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGAAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAcc                                                                                                                               >  1:1163609/1‑151 (MQ=255)
                cTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAccgccg                                                                                                                           >  1:970356/1‑151 (MQ=255)
                cTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGCGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTAccgccg                                                                                                                           >  1:1203089/1‑151 (MQ=255)
                    gACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATTTTTACCTTCCGCTTACCGCCGATg                                                                                                                        <  1:676151/150‑1 (MQ=255)
                      cGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGAcc                                                                                                                     >  2:1123270/1‑151 (MQ=255)
                       ggatggatGGGCGGCTTTAGGGAGGAGCGTCGTCTTAATGGCTTAAAGCAGAGCGGGCGCATTCGCGAAGGGGTGGGAGAAGTGGGGGGTAAAGGCATCGTGTGCCCGGGTGGGTGGGGGATGTTTCCCTGCGGGTTGCCGCCGATGACCt                                                                                                                    <  2:204525/151‑1 (MQ=255)
                          tggatTGGCGTTTTTATCGAGGAGCGTCGTTTTAATGCCTTAAAGCAGATCGGCCGCATTCTCGAAGCGCTGGAGGAGGTGGGGGGTAAGGCCATCTTCGTTCCGGCTGCGGGGGGAATGTTTACCTTCCGCTTACCGCCGATGACCTcg                                                                                                                  <  2:805228/150‑1 (MQ=255)
                          tggatTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTcgc                                                                                                                 >  2:836475/1‑151 (MQ=255)
                           ggatTGGCGGTTTTATCGAAGGGCGTCGTTTTAATGGCTAAAAGCAGATCGAGCGCATTCACGAAGGGCGGACAGAGGGGGGCGGTAAGGGCATCGCCGTTCGGGCTGGGTGGGGCATGTTTACCTTCCCCTTCCCGCCGATGACCTcg                                                                                                                  <  2:903357/149‑1 (MQ=255)
                                ggCGACTTTATCGAAGAGCGTCGTCTTAATGGTTAAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTa                                                                                                           <  1:1151925/151‑1 (MQ=255)
                                   gACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGcc                                                                                                        >  1:380083/1‑151 (MQ=255)
                                   gACTTTATCGAAGAGCGTCGTCTAAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGc                                                                                                         <  1:387860/150‑1 (MQ=255)
                                      tttATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTgg                                                                                                     >  1:356368/1‑151 (MQ=255)
                                       ttATCGAAGGGCCTCGTCTTAATGGTTTAACGCAGTCCGCCCCCCTCCTCGAGGCCATGCAAAACGTGGCAGGTAAAGGCAGCCTCTTCCCGGTTGCGTGGGGATTTTTTACTTCCCGTTTACCGCCGATGACCTCGCCGCGTAGCCt                                                                                                       <  2:164009/148‑1 (MQ=255)
                                       ttATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTgg                                                                                                     >  2:645454/1‑150 (MQ=255)
                                         aTCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACg                                                                                                  >  1:806675/1‑151 (MQ=255)
                                         aTCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACg                                                                                                  >  1:1087477/1‑151 (MQ=255)
                                         aTCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCGGGCAGAGGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGAc                                                                                                   <  1:223626/150‑1 (MQ=255)
                                          ccGAAGGGCGCCTTCTTAAGGGCTAAAAGCAGATGGGGCGCATTCCCGAGGCGTTGGCAGAAGTGGGCGAAAAGGCCTTCGTCGTCCCGCCTGCGTGGGGCATTTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGAc                                                                                                   <  2:348353/148‑1 (MQ=255)
                                            gAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCg                                                                                               >  1:648967/1‑151 (MQ=255)
                                            gAAGAGCGTCGTCTTAATGGCTAAAAGCAGATCGAGCGCTTTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGa                                                                                                    <  1:96438/146‑1 (MQ=255)
                                             aaGGGCGGCGTCTTAATGGCTGAAAGAAGATCGCGCGCATTCTGGAAGCGCGGGCAGAAGTGGGCGGTAAGGGCACCGTCGTTCCGGCTGCGTGGGGATGGTTTCCCTCCCGCTTCCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGa                                                                                              <  2:1163511/151‑1 (MQ=255)
                                             aaGAGCGTCGTCTTAATGGCTTAAAGAAGATCGAGCGCTTTCTCGAAGGGCGGGCAGAAGTGGGCGGAAAGGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACg                                                                                                  <  2:490732/147‑1 (MQ=255)
                                                 gCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGc                                                                                          >  1:889687/1‑151 (MQ=255)
                                                 gCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGc                                                                                          >  1:1172606/1‑151 (MQ=255)
                                                  cgtcgtcgtAATGGCTTAAAGCAGATCGAGCGAATTCTCGAAGCGCTGGCAGAAGTGGGAGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGATTACCTTCCGCTTACCGCCGATCGCCTCGCCGCGGAGCCTGGACGGCGAcc                                                                                            >  2:437049/1‑148 (MQ=255)
                                                   gtcgtcTTAAGGGCTAAAGGAAGATCGGGCGCTTCCTCGAAGCGCTGGCAGAAGGGGGCGGAAAGGGCATCGTCGTCCGGGTTGGGTGGGGCATGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGAcc                                                                                            <  2:897146/147‑1 (MQ=255)
                                                      gtcTTAATGGCTAAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCaaaa                                                                                      <  1:886792/150‑1 (MQ=255)
                                                            aTGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGgtg                                                                                 >  2:1082435/1‑149 (MQ=255)
                                                             tGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGgtgagt                                                                              >  1:653137/1‑151 (MQ=255)
                                                                 ttAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGAtt                                                                          >  2:367235/1‑151 (MQ=255)
                                                                  tAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTc                                                                         >  1:733598/1‑151 (MQ=255)
                                                                   aaaGCAGGCCGCGCGTTTTCGCGAAGCGCGGGCAGAAGTGGGCGGTAAAGACCTCGTCGTTCCGGCTGCGGGGGGCATTTTTCCTTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGAt                                                                           <  2:437386/148‑1 (MQ=255)
                                                                   aaaGCAGCTCGAGCGCATTCTCGAAGCGCTGGAAGAAGTGGGCGGTAAGGGCATCGTCGTTCCGGCTGCGGGGGGTTTGTTTCCCTTCCCCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGgtgagtga                                                                            <  1:90054/147‑1 (MQ=255)
                                                                   aaaGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGGGGGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGAtt                                                                          <  2:877646/149‑1 (MQ=255)
                                                                    aaGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTccc                                                                       >  1:885073/1‑151 (MQ=255)
                                                                        aGATCGGGCGCATTCTCGGAGCGCTGGCAGAGGGGGGCGTTAAGGGCTTCGTCTTTCCGGCTGTGTGGGGATTGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCt                                                                      <  2:196690/148‑1 (MQ=255)
                                                                        aGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGAAAAGGCATCGTCGTCCCGGCTGCGTGGGGCATGTTTACCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTgcg                                                                   <  2:180351/151‑1 (MQ=255)
                                                                         gagcgggcgcATTCTCGAGGCGCGGGCAGAAGTGGGCGGTAAAGGCATCTTCGTTCGGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTgcg                                                                   <  2:1182075/144‑1 (MQ=255)
                                                                           acggccgtaTTCTCAAAGCGCTGGCTGAAGTGGGGGAAAATGGTGGTTCAGTTCGGCAGGCGTGGGAAATTTTTCCCTACGATTTGCCGCCGATGGCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGcgcg                                                                 <  1:1111595/142‑1 (MQ=255)
                                                                             gAGCGCATTCTCGAGGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCTTCCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGt                                                                <  2:69643/149‑1 (MQ=255)
                                                                             gAGCGCATTCCCGAAGCGCTGGAAGAAGTGGGCGGTAAAGGCACCGTCGTTCCGGCTGCGTGGGGCATTTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGt                                                                <  1:572949/149‑1 (MQ=255)
                                                                              aGCGCATTCTCGAGGGCTTGGCAGAGTTGGGCGATAAAGGCATCGTCGTTCGGGTTGCGGGGGAAGTGTTTCCCTCCCGTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTa                                                               <  2:476209/149‑1 (MQ=255)
                                                                              aGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACt                                                             >  1:570179/1‑151 (MQ=255)
                                                                              aGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTCCGTGGGGCATGTTTACCTTCCGATTACCGCCGACGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGCGAGTGATTCCCTGCGCGTAc                                                              >  1:357252/1‑150 (MQ=255)
                                                                                 acATTCTCAAAGCGTTGGAAGAAGTGGGCGGAAAAGGCATCGTCTTCCCGGCTGCGGGGGGCATTTTTTCCTTCCGTTTCCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGa                                                          <  1:400237/150‑1 (MQ=255)
                                                                                     tctcGAAGCGCTGGCAGAGGTGGGCGGTAAAGGCATCGGCGTTCCGGCGGGGTGGGGCATTTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACa                                                       <  2:325483/150‑1 (MQ=255)
                                                                                     tctcGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAAc                                                        >  1:1179100/1‑149 (MQ=255)
                                                                                       tcGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAgg                                                     >  1:800642/1‑150 (MQ=255)
                                                                                       tcGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAgg                                                     >  2:215094/1‑150 (MQ=255)
                                                                                         gaggcgcTGGCAGAAGTGGGGGGTAAAGCCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGt                                                    <  1:687988/146‑1 (MQ=255)
                                                                                         gAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTcg                                                  >  1:1063085/1‑151 (MQ=255)
                                                                                                 ggaagaggggggTGGTAAAGGCACGGTCTTTCCGGCTGCGTGAGCCATTTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTcgccg                                               <  1:567315/137‑1 (MQ=255)
                                                                                                 ggCAGAAGGGGGCGGAAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGt                                          <  1:907969/151‑1 (MQ=255)
                                                                                                     gAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCg                                      >  1:1138586/1‑151 (MQ=255)
                                                                                                     gAAGTGGGCGGTAAAGGCATCGTCGTCCCGGCGGCGGGGGGCATGTTTCCCTTCCGTTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTa                                         <  1:564808/148‑1 (MQ=255)
                                                                                                       aGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCgg                                     >  2:265137/1‑150 (MQ=255)
                                                                                                        gTGGGCGGTAAGGGCATCGTCGTTCCGGCGGCGTGGGGCATGTTTCCCTCCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGaa                                   <  2:866933/151‑1 (MQ=255)
                                                                                                           ggCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCg                                <  1:902274/151‑1 (MQ=255)
                                                                                                           ggCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCg                                <  1:811920/151‑1 (MQ=255)
                                                                                                            ggggTAAGGGCTTCGTCGTCCCGGCGGCGTGGGGCTTTTTTCCCTCCCGGTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCgt                               <  2:112992/149‑1 (MQ=255)
                                                                                                            gCGGTAAAGGCATCGTCTTTCCGGCTGGGTGGGGCATTTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCgt                               <  1:229116/151‑1 (MQ=255)
                                                                                                                 aaaGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCgtggtg                           <  2:1087424/150‑1 (MQ=255)
                                                                                                                 aaaGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGgtgt                          >  1:416552/1‑151 (MQ=255)
                                                                                                                  aaGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGCACCGTGgtgt                          >  2:946788/1‑150 (MQ=255)
                                                                                                                  aaGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGgtgt                          >  1:46694/1‑150 (MQ=255)
                                                                                                                     ggcttcgccttccgggCTGCGTGGGGCATGTTTCCCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTAt                        <  1:851021/136‑1 (MQ=255)
                                                                                                                      caccgtcgtCCCGGCTGCGTGGGGCATGTTTCCCTTCCCCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTAtct                      <  2:363179/147‑1 (MQ=255)
                                                                                                                        accgtggtccgggCTGGGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTAtct                      <  2:1061549/138‑1 (MQ=255)
                                                                                                                         cgtcgtTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGa                   >  2:997686/1‑150 (MQ=255)
                                                                                                                          gtcgtTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAAc                 >  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)
                                                                                                                                                      |                                                                                                                                          
CCACCGCCTGCGGTGGCTATGACGGATGGATTGGCGACTTTATCGAAGAGCGTCGTCTTAATGGCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATC  >  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: