Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC42.7% 188.8 / 32.9 104T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (48/11);  new base C (0/44);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 5.14e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376359‑1376651
                                                                                                                                                           |                                                                                                                                                  
tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAccc                                                                                                                                                 >  4:167065/1‑150 (MQ=255)
 aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                               >  2:1069782/1‑151 (MQ=255)
  ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                    <  2:91653/145‑1 (MQ=255)
  ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                               <  3:1373208/150‑1 (MQ=255)
  ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                               <  3:962647/150‑1 (MQ=255)
    gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                            <  2:457869/151‑1 (MQ=255)
     ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                           >  2:253808/1‑151 (MQ=255)
     ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                           >  2:1216754/1‑151 (MQ=255)
     ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                           <  2:37159/151‑1 (MQ=255)
      cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCg                                                                                                                                           >  1:274620/1‑150 (MQ=255)
           aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                               <  1:43899/141‑1 (MQ=255)
           aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                               >  4:43899/1‑141 (MQ=255)
           aCCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                     >  1:172018/1‑151 (MQ=255)
            ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                     <  4:172016/150‑1 (MQ=255)
               cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCGACACCCTCGCCGATGCccgc                                                                                                                                   <  4:327312/149‑1 (MQ=255)
               cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTAAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                  <  2:354457/150‑1 (MQ=255)
                gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                                >  2:255823/1‑151 (MQ=255)
                gccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGt                                                                                                                                >  3:179832/1‑151 (MQ=255)
                       aGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcg                                                                                                                         >  3:876055/1‑151 (MQ=255)
                           tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                      >  2:1320873/1‑150 (MQ=255)
                            gggcGGCGCCCGAAAAATGGGGAGTGTTCCCCTGCGCGTACGGGAACGGGTCGCCGCGCGTACCGGACCCGTGGTGTACCTCGACCCGTAAACCCGCTATCAGGATAATATGATCAACCCCCGCGCCGATGCCCGCCGTTACATCGTCGa                                                                                                                     <  2:1192119/147‑1 (MQ=255)
                            ggccGGCGCCCGAAAAAGGGTGAGTGTTCCCCTGCGCTTACTGGAACGGGTCGCCGCGCGTCCGGAACCCGGGGTTTATCTCGACCCGTTAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                      <  2:337717/146‑1 (MQ=255)
                            ggACGGCGACCCAAAAATGGTGGGTATTCCCCTGCTGTTGTGGAAACGGGTCGCCGCGCGTACGGGCACCGGGTTGTATCTGGAACCTTTAACCCGTTACCAGGATAATGTGACCAACCCCCTCCCCGATGCCCGCCGTTACAtcg                                                                                                                         <  4:406602/146‑1 (MQ=255)
                               cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                 >  3:207055/1‑151 (MQ=255)
                               cGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                 >  2:75238/1‑151 (MQ=255)
                                ggcgcccaaaaaaTGGTGTGTGATCCCGGGCGCGTGCTGGAACTGGTCGCGGCGCGTCCCGAAACTGTGGTGTCTCCCGACCCTTAAACCCGCTATAAGGATTTTATGACCACCCCCGCCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                 <  2:981707/141‑1 (MQ=255)
                                ggCGACCAAAAAATGGTGAGTGATTCCCGGCGCGTAGGGGACCGGGTCGCCGCGCGTCCGGGACCCGTGGTTTATCTCAACCCTTAAACCCCCTATCGGGATAATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                 <  2:912451/150‑1 (MQ=255)
                                  gcgcccaaaaaTGGGTGATTGATCCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCGGGACCCGTGGTGTATCTCGACCCATTAAACCGCTACCAGGATCATATGACCAACCCCCTCACCAATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                               <  3:202672/144‑1 (MQ=255)
                                   gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCGCGCCGTTACATCGTCGAAAACGATc                                                                                                             >  3:297602/1‑151 (MQ=255)
                                   gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                             >  2:145110/1‑151 (MQ=255)
                                    cccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                            >  2:1063236/2‑151 (MQ=255)
                                     cccaaaaaTGGTGAGGGATCCCCGGCGCGTACGGAAACGGGTCGCCGCGCGTACGGAACCCGGGGTGTATCTCGACCCTTTAACCCGCTATCAGGATCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                           <  4:84238/147‑1 (MQ=255)
                                     ccGAAAAAGGGTGAGTTTTCCCCTGCGCGTACGGGAACGGGTCGCCGCGCGTCCCGAACCCGTGGTGTATCTCGACCCTTAAACCCGCTATCAGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                            <  2:87743/150‑1 (MQ=255)
                                      cgaaaaaTGGTGAGTGTTCCCTTGCGCTTCCTGGAAGAGGTCCCGGCGCGTCCCGGAACGGGGGTATTTCTGGAACCTTAAACCCGTTATAGGGATCATATGAAAAACCCCCTCGCGAAGGCCCGCCGTTACATCGTCGAAAACGa                                                                                                               <  2:878828/143‑1 (MQ=255)
                                      cgaaaaaTGGTGAGTGATCCCCGCCGCGTACGGAAACGGGCCCCCCCGCGTCCCGGACCCGGGGTGTCTCTCGACCCTTAAACCCGCTCTCAGAATAATATGAAAACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                               <  3:627089/143‑1 (MQ=255)
                                        aaaagTGGGGAGTGTTCCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCCGTGTTGTATCCCGAACCGTAAACCCTCTATCAGGATAATATGAAAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                        <  3:1143132/150‑1 (MQ=255)
                                         aaagTGGTGAGTGATTCCCGGCGCGTACGGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTTTATCTCGACCCTTTAACCCGTTATCAGGATCATAGGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                        <  2:291465/150‑1 (MQ=255)
                                          aagagggtgggtgATTCCCGGCGGGTAGGGGAACGGGCCGCCGGGCGTCCCGGACCGGGGTTGTACCCCGAACCGTTAACCCGTTACAAGGATAAAATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                        <  4:91454/145‑1 (MQ=255)
                                           aaagggggagtacccccccGCGCGCACGGGAACCGGCCGCCGCGCGCCCCGGACCAGTGGTGTATCGCAACCCTTCAACCAGCGCGACGGATCATATGAACAACCCCCCCCCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACa                                                                                                      <  3:338936/136‑1 (MQ=255)
                                            aatgtggagtgtTCCCCTGCGGGTGCGGGAACAGGTCGCCGCGCGTCCGGAACCGGTGTTGAATCACGAACCGTAAACCCCTTTACAGGATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                     <  2:1135009/145‑1 (MQ=255)
                                               ggtgggtgATTCCCTGCGCGTACTGCACCGGGTCGCCGCGCTTCCCGGAACCGTGGTTTACCCCAACCCTTAAACCCGTTATCAGGACTGTATGTTAAACCCCCTCGCAGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTa                                                                                                  <  2:975396/151‑1 (MQ=255)
                                               ggTGAGTGTTTCCCTGCGCGTACGGGAACGGGCCGCCGCGGGTCCCGAAACGGTGGTGTATCCCGACCCTTAAACCCGTTATCAGGATAAAATAACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                     <  4:437669/148‑1 (MQ=255)
                                                 tgagtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                 >  2:166145/1‑150 (MQ=255)
                                                  gagtgaTCCCCTGCGCGTACGGGAACGGGTCGCCGCGCGTCCGGGACCCGGGGTTTATCTCAACCCGTTAACCCGTTATCAGGATCATATGACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                <  2:151801/150‑1 (MQ=255)
                                                     tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                            >  2:80240/1‑151 (MQ=255)
                                                     tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                            >  3:1002424/1‑151 (MQ=255)
                                                          cccTGCGCGTACGGGAACGGGTCGCCGCGCGTCCCGGACCCGTGTTGTATCTCGACCCGTTAACCCGTTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                       <  1:87026/151‑1 (MQ=255)
                                                           ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGTTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                      >  1:302321/1‑151 (MQ=255)
                                                               cgcgTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                  >  3:629537/1‑151 (MQ=255)
                                                                 gcgtagggaaaCAGGCCGCCGCGCGTACGGGACCCGTGTTTTATCTCAACCCTTAAACCCGTTAAAAGAATAATATAAACAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                 <  3:1273298/142‑1 (MQ=255)
                                                                  gTACTGGAATAGGCCGCCGCGCGTACCGTATCCGTGTTGTTTCCCGATCTTTTAACCCGCTATCAGGTTGATATAATCAACCCCCTCTCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                     <  2:120809/145‑1 (MQ=255)
                                                                     aatgaaaggggTCGCCGCGCGTCCCGAACCCGTGGTTTTCCTCAACCCGTAAACCCGTTATAAGGATCATAGAACAAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                            <  3:24976/127‑1 (MQ=255)
                                                                       ggAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                          >  1:292617/1‑151 (MQ=255)
                                                                        gAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcac                                                                          >  3:1143565/1‑150 (MQ=255)
                                                                         aaaCATGTCGCCGCGCGTACCGGACCCGTGGTGTACCTCGACCCTTTAACCCGCTATCAGAATCATATTATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                           <  1:380230/148‑1 (MQ=255)
                                                                          aaCAGGTCGCCGCGCGTACCGGACCGGTGTTGTATCTCAACCCGTTAACCCGCTATCAGGATCATATGATCAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                             <  1:197578/146‑1 (MQ=255)
                                                                              ggcccccccgccccccGGACCCGGGGTGTACCTCAATCCTTAAACCCCTTACGAGGTTCATATGATCACCCCCCTCGCCGACTCCCGCCGTTACCTCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGaa                                                                     <  4:106521/136‑1 (MQ=255)
                                                                              ggTGGCCGCGCGTCCGGGACCCGTGGTTTACCCCGACCCGTAAACCCGCTAGCAGGATAATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                      <  3:204055/149‑1 (MQ=255)
                                                                                      gcgcttccGGAACCGGTGGTGTATCTCGACCGTTAAACCCCTAATCAGGATCATATGACAAACCCCCCCGCCGTTGACCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                              <  1:90119/144‑1 (MQ=255)
                                                                                      cgcgTCCCGGACCCGGGGTTTATCCCGACCCGTTAACCCGCTATCAGGATCATATAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                             <  2:101480/150‑1 (MQ=255)
                                                                                      cgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                               >  2:9715/1‑148 (MQ=255)
                                                                                         gtccgggAACCGTGGTTTATCTCGACCCTTTAACCCGCTATCAGAATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                           <  1:478858/144‑1 (MQ=255)
                                                                                            cgcgggaccggttgtgtATCTCAACCCTTAAACCCTTTATCAGAATAATATCACCACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                           <  3:828585/135‑1 (MQ=255)
                                                                                            ccGGACCGGTGGTTTATCCGGACCTGTAAACCCGCTATAAGGTTAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                      <  4:454892/151‑1 (MQ=255)
                                                                                              ggAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                    >  3:261411/1‑151 (MQ=255)
                                                                                                aaCCGTGGTGTATCTCGACCCGTAAACCCTTTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                   <  2:637765/150‑1 (MQ=255)
                                                                                                aaCCGTGGTGTATCTCGAACCGTTAAACCGCTATCGGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                  >  4:95554/1‑151 (MQ=255)
                                                                                                aaCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTg                                                  >  2:704693/1‑151 (MQ=255)
                                                                                                  accggtgttttATCCCGACCCTTAAACCCGCTATAAGGATCATTTGATCAACCCCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                   <  3:54054/140‑1 (MQ=255)
                                                                                                   ccGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCACCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgc                                                >  1:319160/1‑151 (MQ=255)
                                                                                                          tttATCCCGACCCGTTAACCCGCTACCAGAATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                           <  3:1069828/148‑1 (MQ=255)
                                                                                                              tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                      >  2:1304471/1‑151 (MQ=255)
                                                                                                               ctgaaccctttAAACCGCTATAGGGATAAAATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                           <  4:345512/136‑1 (MQ=255)
                                                                                                                 caaaCCTTAAACCCGCTATCAGGATAATATGAACACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                   <  3:75248/149‑1 (MQ=255)
                                                                                                                 cGAACCTTTAACCCGCTATCAGAATCATATGACCAACACCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                         <  2:629511/145‑1 (MQ=255)
                                                                                                                  gAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                  >  2:919690/1‑151 (MQ=255)
                                                                                                                   aaCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                  >  1:370656/1‑150 (MQ=255)
                                                                                                                    aCCGTAAAACCGTTATCAGAATCATATGAAAAACCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                    <  3:1067301/147‑1 (MQ=255)
                                                                                                                    aCCCTTAAACCCGCTATCAGAATAATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCt                                       <  2:876020/145‑1 (MQ=255)
                                                                                                                     cctttAAACCCCTATAAGTATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                  <  3:50599/145‑1 (MQ=255)
                                                                                                                      cGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCACCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACaa                              <  4:319157/151‑1 (MQ=255)
                                                                                                                        ttAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                            >  3:281294/1‑151 (MQ=255)
                                                                                                                          aaaCCGTTATAAGAATAATAAGAACACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                          <  4:63238/151‑1 (MQ=255)
                                                                                                                              ccGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAc                       >  1:428259/1‑151 (MQ=255)
                                                                                                                                 gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                     >  2:907921/1‑151 (MQ=255)
                                                                                                                                  cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                    >  1:35287/1‑151 (MQ=255)
                                                                                                                                      aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                   >  3:351803/1‑150 (MQ=255)
                                                                                                                                       tCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                 >  1:336973/1‑151 (MQ=255)
                                                                                                                                           ggATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt              >  2:1291110/1‑151 (MQ=255)
                                                                                                                                              tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt           >  3:1208708/1‑151 (MQ=255)
                                                                                                                                                atatGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt         <  1:408859/151‑1 (MQ=255)
                                                                                                                                                atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt         >  4:156772/1‑151 (MQ=255)
                                                                                                                                                atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt         >  2:28453/1‑151 (MQ=255)
                                                                                                                                                  atGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                     <  3:80252/137‑1 (MQ=255)
                                                                                                                                                   taacaaaCCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                     <  3:151172/131‑1 (MQ=255)
                                                                                                                                                   tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc       >  3:617914/1‑150 (MQ=255)
                                                                                                                                                   tGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg        <  2:1002382/149‑1 (MQ=255)
                                                                                                                                                    gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat     >  3:1229223/1‑151 (MQ=255)
                                                                                                                                                      tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCatat   >  3:1352882/1‑151 (MQ=255)
                                                                                                                                                        aaCACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt  >  3:58060/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                                  
TACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376359‑1376651

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