Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,376,5050AC39.7% 123.0 / 19.8 75T158P (ACC→CCC) ycjRputative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (37/7);  new base C (0/29);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 4.39e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376355‑1376652
                                                                                                                                                        |                                                                                                                                                   
cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGACGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAaca                                                                                                                                                     <  2:523202/151‑1 (MQ=255)
  cTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACAcc                                                                                                                                                   >  2:1153295/1‑151 (MQ=255)
   ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAacac                                                                                                                                                    <  2:461229/149‑1 (MQ=255)
     aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACTCCCTc                                                                                                                                                <  2:209301/151‑1 (MQ=255)
        gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                             <  1:722207/151‑1 (MQ=255)
            aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                         >  1:1305887/1‑151 (MQ=255)
            aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                         <  2:530264/151‑1 (MQ=255)
                ccTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAGCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGccc                                                                                                                                      <  1:11010/150‑1 (MQ=255)
                     ccgcGTAGCCTGGACGGCGACCGCAAAATGGGGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCCCCGATGCCCGCCGtt                                                                                                                                >  1:1386100/1‑151 (MQ=255)
                       gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                              >  2:959154/1‑151 (MQ=255)
                                    ggCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                                 >  1:108364/1‑151 (MQ=255)
                                        aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                             >  2:1389325/1‑151 (MQ=255)
                                        aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                             >  1:1079784/1‑151 (MQ=255)
                                        aCCGAAAAAGGGTGAGTGTTTCCCTGCGCTTATTGAAACAGGTCGCCGCGTGTCGCGGAACGGTGGTGTATCTCGAACGGTTAACCCGCTATCAGGTTCATTTGATCAACACCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               <  2:43664/149‑1 (MQ=255)
                                         ccGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGAACCTTTAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCt                                                                                                             <  2:238607/150‑1 (MQ=255)
                                         ccGAAAAGTGGTGAGTGATCCCCTGCGCGTCCTGGAACGGGCCGCCGCGCGTACCGGAACGGGGGTGTATCTCGAACCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCtt                                                                                                            <  2:1244359/151‑1 (MQ=255)
                                             aaatgggTGAGTGATTCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTGTATCCCGACCCGTTAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                        <  2:36123/146‑1 (MQ=255)
                                              aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                        >  1:979741/1‑150 (MQ=255)
                                              aaaGGGTGAGTTGTCCCCGGCGCTTACTGAAACAGGGCGCCGCGCGTCAGGGACCCGTGGTTCATCCCGAACCGTCAACCCGTTATGAGGTAAATAACATCACCACCGTCGCGGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                               <  1:151918/143‑1 (MQ=255)
                                               aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                      >  2:11235/1‑151 (MQ=255)
                                                     gagtgaTTCCCTGCGCGTACTGGAAGAGGTCGCCGCGCGTCCGGGAACCGTGGTGTATCTCAAACCGTTAACCCGCTATCAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  <  2:87170/149‑1 (MQ=255)
                                                     gagtgaTCCCCTGCGCTTACTGAAACGGGTCGCCGCGCGTACGGAAACCGTGGTGTTTCTCAACCCTTAAAACCGCTATCAGAAAAATATGACCAACCCCCTCGCGGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                  <  2:987353/149‑1 (MQ=255)
                                                         gtttCCCTGGGCGTAGTGGAACAGGTCGCCGCGCGTACCGGACCCGTGGTTTATCACGACCCGTTAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                             <  1:783959/148‑1 (MQ=255)
                                                           tccccTGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGGACCCGTGGTTTATCTCGAACCTTTAACCCGCTATAAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                           <  2:585051/148‑1 (MQ=255)
                                                             cccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                        >  1:405462/1‑151 (MQ=255)
                                                             cccGGCGCTTACTGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTTTATCTCGAACCGTTAACCCGCTATCAGGATCATATGATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCg                                                                                         <  2:90393/150‑1 (MQ=255)
                                                                 gcgcGTACTGGAACAGGTCGCCGCGCGTCCCGGAACCGTGGTGTATCTCGACCCGTTAACCCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                    <  2:1307909/151‑1 (MQ=255)
                                                                 gcgcGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                    >  2:1075914/1‑151 (MQ=255)
                                                                   gcgTACTGGAACAGGCCGCCGCGCGTACCGGACCCGGGGTGTATCCCGACCCGTTAACCCGCTATCAGGATCATATGATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAttt                                                                                  <  2:324287/151‑1 (MQ=255)
                                                                        cTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                             >  2:748609/1‑151 (MQ=255)
                                                                          ggAACAGGTCGCCGCGCGTACCGGACCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                            <  2:298612/150‑1 (MQ=255)
                                                                              cAGGCCGCCGCGGACCCCGGACCCGGGGTGTACCCCGACCGTTAAACCCTCTAAAAGGATAATAGAACCACCCCCCTGGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATca                                                                            <  1:196715/146‑1 (MQ=255)
                                                                               acgTCGCCCCGCGTCCGGGAGCTGTGTTTTACCTCAAACTTTTAACCCGCTATAGGGATAATAGAAACACCCCCCTCCAAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                        <  1:222082/147‑1 (MQ=255)
                                                                                ggTCGCCGCGCGTCCGGGACCCGTGGTGTATCTCGACCTGTTAAACCGCTATCAGGATATTATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                     <  2:1107060/151‑1 (MQ=255)
                                                                                 gTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                     >  1:935910/1‑150 (MQ=255)
                                                                                     ccgcGCGTACCGGAACGGTGGTGTATCTCGACCCTTTAACCCGCTATCAGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                    <  2:526528/147‑1 (MQ=255)
                                                                                      cgcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                >  1:686828/1‑150 (MQ=255)
                                                                                       gcgctTCCCGGACCCGTGGTTTATCTAAACCCTTAAACCCCCTATCAGGATAATTTAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                 <  2:950340/148‑1 (MQ=255)
                                                                                           gTACCGGACCCGTGGTGTATCTCGACCCTTAAACCCGCTATCAGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGaagaag                                                           <  2:128182/150‑1 (MQ=255)
                                                                                           gTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  2:441101/1‑148 (MQ=255)
                                                                                            tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  1:1104099/1‑147 (MQ=255)
                                                                                            tACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                             >  2:220413/1‑147 (MQ=255)
                                                                                              ccGGACCCGTGGTGTATCGCGACCGGTTAACCCTCTAGCAGAATCATATGATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaagaa                                                            <  1:1049256/146‑1 (MQ=255)
                                                                                              ccGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                        >  2:13631/1‑150 (MQ=255)
                                                                                               cGGAACCCTTGTTGTATCTCGAACGTTAAACCCGTAACCAGGATAAGATGAAAAACCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                         <  1:765945/149‑1 (MQ=255)
                                                                                                ggaccggtggtgTATCTCCACCCGTTAACCCCCTACCAGATTCATATAATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGAt                                                         <  2:1223651/141‑1 (MQ=255)
                                                                                                ggAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                     >  1:914233/1‑151 (MQ=255)
                                                                                                ggAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                     >  1:1388347/1‑151 (MQ=255)
                                                                                                gaacccGTGGTTTATCTCGAACCTTTAAACCTCTATCAGGATAAAATGATCACCACCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                              <  1:953768/138‑1 (MQ=255)
                                                                                                  acccGTGTTTTATCTCAACCCTTTAACCCGCTATGAGGATCATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                      <  2:912279/146‑1 (MQ=255)
                                                                                                   aCCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                  >  1:2250/1‑151 (MQ=255)
                                                                                                      gtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               >  2:1016675/1‑151 (MQ=255)
                                                                                                      gtggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGCCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               >  1:919988/1‑151 (MQ=255)
                                                                                                         tggtttacctcaaCCCTTAAACCCTCTATCAGGATCATATGACCACCACCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                               <  1:45651/139‑1 (MQ=255)
                                                                                                            ttATCTCGACCCGTCAACCCGCTTACATGATCATATAACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                           <  1:1085049/149‑1 (MQ=255)
                                                                                                                cgcGAACCGTTAAACCGTTACCAGCATCAAATGATCACCCCCCTCGACGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgcg                                          <  1:1044594/145‑1 (MQ=255)
                                                                                                                       cGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                >  2:596655/1‑150 (MQ=255)
                                                                                                                            aacccGCTATCAGGATCATATGATAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                             <  2:1322966/146‑1 (MQ=255)
                                                                                                                            aaCCGCTATGAGCATCTTATAAACAACCCCCGCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           <  2:19535/150‑1 (MQ=255)
                                                                                                                            aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgcg                          >  1:1386084/1‑151 (MQ=255)
                                                                                                                            aaCCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcgc                           >  1:199861/1‑150 (MQ=255)
                                                                                                                             aCCGCTAGCAGGACCATTGGATAACCTCCCGCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAAGATCGAAGAAGATAGCCTGGCGCAGGCGCTGCATGACAAca                             <  1:616484/147‑2 (MQ=255)
                                                                                                                                cGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  1:81725/1‑151 (MQ=255)
                                                                                                                                 gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGAcc                       >  1:368701/1‑150 (MQ=255)
                                                                                                                                  cTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                     >  2:147176/1‑151 (MQ=255)
                                                                                                                                    aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGc                    >  1:406314/1‑150 (MQ=255)
                                                                                                                                       aGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGAGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                <  2:88185/151‑1 (MQ=255)
                                                                                                                                       aGGATAATATAATCAACCCCCCCGCCGAGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATAACGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg                <  1:1152488/151‑1 (MQ=255)
                                                                                                                                        ggATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc                  <  2:1306749/148‑1 (MQ=255)
                                                                                                                                         gATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc              >  2:396813/1‑151 (MQ=255)
                                                                                                                                           tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt            >  1:988781/1‑151 (MQ=255)
                                                                                                                                            cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg           >  1:643685/1‑151 (MQ=255)
                                                                                                                                             atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgt          >  1:923300/1‑151 (MQ=255)
                                                                                                                                                tGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc        <  2:225377/150‑1 (MQ=255)
                                                                                                                                                tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc        >  1:38916/1‑150 (MQ=255)
                                                                                                                                                tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc        >  1:569142/1‑150 (MQ=255)
                                                                                                                                                 gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCa       >  1:1253994/1‑150 (MQ=255)
                                                                                                                                                      acacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg  >  2:439316/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTG  >  NC_000913/1376355‑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: