Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC45.7% 84.8 / 36.8 70T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (34/4);  new base C (0/32);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 6.83e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376359‑1376652
                                                                                                                                                       |                                                                                                                                                   
tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                  <  2:1319083/150‑1 (MQ=255)
    gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             >  2:205919/1‑151 (MQ=255)
    gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             <  1:658430/151‑1 (MQ=255)
      cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                            >  2:1115824/1‑150 (MQ=255)
         tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                         <  2:54698/150‑1 (MQ=255)
          gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                       >  2:1262373/1‑151 (MQ=255)
             cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                    >  2:917898/1‑151 (MQ=255)
             cTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgc                                                                                                                                    >  1:1289286/1‑151 (MQ=255)
               cgccggGTTGCCTGTACGGCGACCCCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTCCCGGAC‑C‑T‑GTTGTGTACCGCGACCCGTTAAA‑CCGCTATCAGGATCATATGATCAACCCCCTCGCCTATGCccgcc                                                                                                                                   <  2:862636/150‑1 (MQ=255)
                gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCCCC‑GCGCCTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                   <  2:185403/149‑1 (MQ=255)
                  cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                >  2:1324785/1‑150 (MQ=255)
                   gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                              >  2:461932/1‑151 (MQ=255)
                      tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTAGCCGATGCCCGCCGTTACAtc                                                                                                                           >  2:951640/1‑151 (MQ=255)
                       aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                          >  1:207611/1‑151 (MQ=255)
                            ggACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                      >  2:632743/1‑150 (MQ=255)
                                ggCGCCCGAAAAATGGGGGGTGATCCCCTCGGCGTAGTGGAACGGGCCCCC‑GCGCTTCCGGGAC‑C‑C‑GGGGTATATCTCGAACCTTTAAC‑CCGCTAACAGAATCATAAGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                  <  1:554772/150‑1 (MQ=255)
                                 gCGACCGAAAAGTGGTGAGTGATCCCCTGCGCTTACTGGACCGGGTCGCC‑GCGCGCCCCGGAC‑C‑C‑GGGGTGTATCTCGAACCGTTAAC‑CCGCTATCAGGATCATAGGATCACCCCCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACg                                                                                                                 <  1:844279/150‑1 (MQ=255)
                                     ccGAAAATTGGTGATTGATTCCCGGCGCTTACGGAAACGGGTCGCG‑GCGCGTCCCGACA‑C‑G‑GTGGTGTACCGCGACCCGTAAAC‑CCGCTACCAGGATAATATGACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              <  2:655893/149‑1 (MQ=255)
                                     ccGAAAAAGGGTGAGTGATCCCCTGCGCGTACGGGAACGGGTCGCC‑GCGCGTACCGGAC‑C‑C‑GTGGTGTACCTCGACCCGTTAAC‑CCGCTATCAGGATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               <  1:1177619/148‑1 (MQ=255)
                                         aaaaGGGTGGGTTATCCCCTGGGCGTACGGAAACAGGTCGCC‑GCGCGTCCCGGAC‑C‑G‑GTGGTGTTTCTCGAACCGTAAAC‑CCCTTATCAGAATAATAGAACCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTa                                                                                                           <  1:848253/148‑1 (MQ=255)
                                         aaaaGGGTGAGTGATCCCCTGCGCGT‑CCGGGACAGGGCCCCGGCGCGT‑CCGGAA‑C‑CGGTGTTGTTTCTCAACCTTTTAAC‑CCGTTATCAGTATAATATGATAACCCCCCTCGCGGAGGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                              <  1:361300/145‑1 (MQ=255)
                                          aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                        >  2:529914/1‑150 (MQ=255)
                                              ggggggtttttCCCGGCGTGTTGTGGAGCGGGCCGCC‑GCGCCAACGGGAG‑C‑C‑GTGGTTTATCTGGATCCATAAAC‑CCGTTATAGGGATTATATAACAACCCCCCTCCCCGTCGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                   <  2:801600/142‑1 (MQ=255)
                                                 gagtgaTTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  >  1:436769/1‑149 (MQ=255)
                                                   gtgttccccGGCGCGTACGGGAACAGGTCGCC‑GCGCGTACCGGAC‑C‑G‑GGGGTGTATCTCGAACCGTAAAC‑CCGTTATCAGGATCATAGAACCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGa                                                                                               <  1:662569/144‑1 (MQ=255)
                                                      aTTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                           >  2:869087/1‑151 (MQ=255)
                                                       tccccTGCGCGTACTGAAACGGGTCCCC‑GCGCGTACCGGAA‑C‑G‑GGGGTTTTTCTCAACCCGTTAAC‑CCGCTATCAGGATAATATGACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             <  1:412576/146‑1 (MQ=255)
                                                         cccTGCGCGTCCGGAAACGGTGCCCC‑GCGCGT‑CCGGAACC‑C‑GGGGTTTTTCTCGACCCTTTAAC‑CCGCTATAAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                 <  1:730344/142‑1 (MQ=255)
                                                          ccTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                         >  1:237225/1‑149 (MQ=255)
                                                           cTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                      >  2:517663/1‑151 (MQ=255)
                                                                cgTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                 >  2:803634/1‑151 (MQ=255)
                                                                   aCTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              >  1:654530/1‑151 (MQ=255)
                                                                   aCTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                              >  2:809146/1‑151 (MQ=255)
                                                                    cgggAACGGGCCGCC‑G‑GGGTACCCGAA‑C‑CGGTGGTGTTCCCCGACCCGTAAAC‑CCCCTACAAGGATAATATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                    <  2:559111/142‑1 (MQ=255)
                                                                    cTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  2:404540/1‑151 (MQ=255)
                                                                       gAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                           >  2:95979/1‑150 (MQ=255)
                                                                                 cc‑gcGCGTACCGAAA‑C‑C‑GTGGTGTTTCTCGACCCGTTAAC‑CCGCTATCAGAATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  1:544938/147‑1 (MQ=255)
                                                                                 cc‑gcGCGT‑CCGGAA‑CGC‑GTGTTTTTCCCCGACCCTTAAAC‑CCGCTATAAGAATAAAATAATCACCCCCCCCGCGGTGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                                   <  2:830211/148‑1 (MQ=255)
                                                                                  c‑gcgcgT‑CCGGAA‑C‑CTGTGGTGTATCTCAACCCTTTAAC‑CCTCTATCAGGATCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                               <  1:594228/151‑1 (MQ=255)
                                                                                      gcgTCCCGGAA‑C‑G‑GTGTTGTTTCCCAACCCTTTAAC‑CCGCTATAAGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                            <  1:454609/151‑1 (MQ=255)
                                                                                      gcgTCCCGGAA‑C‑C‑GTGGTTTATCTCGAACCGTAAAC‑CCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                               <  1:390389/148‑1 (MQ=255)
                                                                                      gcgTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:265979/1‑150 (MQ=255)
                                                                                        gtcccGGAC‑C‑G‑GGGGTGTACCTCG‑ACCTTTAAACCCGCTATCAGAAAAAAAAGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  2:174093/144‑1 (MQ=255)
                                                                                        gccccGGAC‑C‑G‑GTGTTGTACCCCGAACCGTAAAC‑CCGCAATAAGAATAAAATAAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  2:409340/147‑1 (MQ=255)
                                                                                        gaaCCGAAC‑C‑C‑GTGGTTTATCTCAACCCGTAAAC‑CCGCTACAAGGATAATATAAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                               <  1:1298208/144‑1 (MQ=255)
                                                                                        gTACCGGAC‑C‑G‑GGGTTGTATC‑CCAACCCTTAAACCCGCTATCAGAACAATATGCCCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             <  2:1185154/148‑1 (MQ=255)
                                                                                         caCCGAAA‑C‑C‑GTGGTGTTTCCCG‑ACCGTTAAACCCGCTTAGAGAACCATATGATCACCCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  2:465404/139‑1 (MQ=255)
                                                                                           ccGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                        >  1:1053325/1‑150 (MQ=255)
                                                                                            cGGAA‑C‑CGGTGGTTTATCCCGACCGTTTAAC‑CCGTTATCAGGATCATATAATAAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                         <  1:1063700/149‑1 (MQ=255)
                                                                                                  c‑c‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 >  2:313912/1‑151 (MQ=255)
                                                                                                    cgttgttttATCTCAACCCTTTAAA‑CCGTTCACCGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  2:597202/143‑1 (MQ=255)
                                                                                                       tggttTATCTGAACCCTTTAAC‑CCGCTATCAGAATCTTATGATCAACCCCCCCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                 <  1:881941/148‑1 (MQ=255)
                                                                                                          tttaccccaacccGTTAAC‑CCCCTATCAGGATAATAAGACCCCCCCCCTCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCa                                              <  2:203050/137‑1 (MQ=255)
                                                                                                           gtATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTACACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                            >  2:956156/1‑149 (MQ=255)
                                                                                                                tcGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGc                                      >  1:1060573/1‑150 (MQ=255)
                                                                                                                      ctttaaagccGCTATAGGGCTCATATAAGCACCCCCCTCGTCGATTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 <  1:140191/142‑1 (MQ=255)
                                                                                                                      cGTTAAA‑CCGCTATCAGGATCATCTGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGTGCTGCATGACa                                >  1:278179/1‑150 (MQ=255)
                                                                                                                         taaacccGTTATAAGGATAAAATAAACAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                 <  1:1116762/142‑1 (MQ=255)
                                                                                                                          aaa‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGCTTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             >  2:464809/1‑149 (MQ=255)
                                                                                                                            a‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  1:748340/1‑150 (MQ=255)
                                                                                                                              ccGATATCAGAATCATAAGACACCCCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             <  2:52667/146‑1 (MQ=255)
                                                                                                                                   accAGGATAATATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       <  2:207605/145‑1 (MQ=255)
                                                                                                                                      agaaTCATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  <  1:18338/146‑1 (MQ=255)
                                                                                                                                          tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt            >  2:584889/1‑151 (MQ=255)
                                                                                                                                           cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg           >  1:1033755/1‑151 (MQ=255)
                                                                                                                                              ataaccAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCa             <  2:1063121/141‑1 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat    >  1:153977/1‑151 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat    >  2:586467/1‑151 (MQ=255)
                                                                                                                                                  tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCCCGGTCATGTGCatat    >  1:692364/1‑151 (MQ=255)
                                                                                                                                                    aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  1:432741/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCC‑GCGCGTACCGGAA‑C‑C‑GTGGTGTATCTCGAACCGTTAAA‑CCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376359‑1376652

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