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

CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCG  >  NC_000913/1376361‑1376654
                                                                                                                                                           |                                                                                                                                                     
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTc                                                                                                                                                  <  1:1027982/150‑1 (MQ=255)
ccgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCg                                                                                                                                                 <  2:1369189/151‑1 (MQ=255)
 cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGc                                                                                                                                                <  1:535295/151‑1 (MQ=255)
     gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                            >  1:1119449/1‑151 (MQ=255)
             cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgccg                                                                                                                                    >  2:65526/1‑151 (MQ=255)
               ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                                  >  1:131879/1‑151 (MQ=255)
                       ccTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                          >  1:374747/1‑151 (MQ=255)
                         tGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtcgtcg                                                                                                                         >  2:680953/1‑150 (MQ=255)
                                cGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                                 >  1:1220077/1‑151 (MQ=255)
                                   ccacaaaggggTGAGTCATTCCCGCCGCGTGCGGCAACGGGCCGCCGCGCGCCCCGAACCGGGGTTTTATCTCGAACCATAAACCCGTTATCAGCATTATATAACAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                                  <  1:40656/138‑1 (MQ=255)
                                       aaaaTGGGGAGTGCTCCCCTGCGTTTTCGGAAACGGGCCCCCGGGGGCCCCGAACCCTGGTTTTATCTCGAACGGTTAACCCGCTACCGGCATCATATGAAAAACCCCCTCGCCTCTGCCCGCCGTTACATCGTCGAAAACGATCTTaaa                                                                                                           <  1:1358247/150‑1 (MQ=255)
                                        aaaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAAc                                                                                                          >  1:5715/1‑150 (MQ=255)
                                           tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                      >  1:531589/1‑151 (MQ=255)
                                                agtgtTTCCGTGCGTTTTCGGAAACGGGCCGCCGCGGTTCCCGGACCCGTGGTGTATCAGGAACCGTTAACCCGCTAACAGGATCATATGACCACCCCCCCCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                                  <  1:1272019/150‑1 (MQ=255)
                                                agtgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                                    >  2:424804/1‑148 (MQ=255)
                                                 gttatccccTGCGCTTCCGGAAACGGGCCGCCGCGCGTCCCGGACCCGTGGTGTATCCCAACCCTTTAACCCGCTATAGGTATCATATGACAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                                  <  2:878039/143‑1 (MQ=255)
                                                 gtgttccccTGCGCTTCTGGAAACGGGCCCCCGCGCGTCCCGAACCGGGGGTTTATCTCAAACCTTAAACCCGCTATAGGGATAATATGACCACCCCCCCCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                        <  2:1341872/137‑1 (MQ=255)
                                                       cccggggcgTACGGAACCGGTCCGCCGCGCGCCCGGGACCGGTGTTGTACCCCAACCCGTAAACCCGCTATCAGAATAATATGACAACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                          <  2:1344882/145‑1 (MQ=255)
                                                       cccTGCGCGTCCTGAAAAGGGTCCCCGGGCGCACGGAACCGGGGGTGTACCCCGACCCTTAAACCCGCTATAAGGATAAAATAACAACCCCCCCCGCGGAGCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCgg                                                                                          <  2:371114/151‑1 (MQ=255)
                                                       cccTGCGCGTACTGAAACGGTTCGCCGCGCGTCCGGGAACGGTGGTTTATCTCGACCCGTTAACCCGCTATGAGGATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                             <  1:1244616/148‑1 (MQ=255)
                                                       cccGGCGCGTACGGAAACAGGCCGCCGCGCGTCCCGGACCGGTGTTGTCTCCCCACCCGTAAAACCGCTATCAGAATTATATGATCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                             <  2:1261498/148‑1 (MQ=255)
                                                         cTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                        >  1:1216828/1‑151 (MQ=255)
                                                          tGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCg                                                                                        >  2:362748/1‑150 (MQ=255)
                                                           gcgcGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                      >  2:729844/1‑151 (MQ=255)
                                                           gcgcGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                      >  1:574484/1‑151 (MQ=255)
                                                               gTACTGAACCGGGCCCCCCCGCGTCCGGAACCGGTTGTTTATCCCAACCCTTAAACCCCTAATAAGTATAATATAAACACCCCCCCCCCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                   <  2:1269527/150‑1 (MQ=255)
                                                                 aCTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                                >  2:46712/1‑151 (MQ=255)
                                                                     gaaaCCGGTCGCCGCGTGTGCGGGACGGGGGTTGTTTCCCGACCGTTAAACCCGTTATAGGATTCATATGACAAACCCCCGCCCGGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTAt                                                                                <  2:1309513/146‑1 (MQ=255)
                                                                        acggttCGCCGCGGGTCCGGGAACCGTGCTGTATCTCAACCCTTAAACCCGCTATCAGTATAATATGACCACCCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                            <  1:755374/144‑1 (MQ=255)
                                                                        aacagtcccccgcGTGTCCCGGACCCGGGGGGTCTCTCGACCCTTTAAACCGTTACAAGATTAATATGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTc                                                                                   <  1:1122635/134‑1 (MQ=255)
                                                                         acggtcccccccgcgTCCGGGACCTGTGGTGTTCCCCGAACCTTAAAACCCCTATAGGTTTCAAATGACAACCCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                           <  1:1105179/139‑1 (MQ=255)
                                                                          cgggcccccgcGGGTACCGGACCCGGGGTGTATCCCGCCCCTTTAACCAGCTATAAGATAAAAGTAACCACCACCCGGGCGAAACCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                         <  1:49178/144‑1 (MQ=255)
                                                                            ggcccccgccgcccccGGACCCGTGGTTTCTCTCGACCCTTAAACCCGCTATCAGAATAATATAAACACCCCCCTCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATg                                                                          <  2:245918/134‑1 (MQ=255)
                                                                               cgcccggcgccccGGACCGGTGGTGTATCTCGACCTTTAAACCGTCTATAAGAATCATATAACAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                       <  2:264370/137‑1 (MQ=255)
                                                                                 ccgcGCGTCCCGGACCCGTGGTGTTTCCCGACCCTTAAACCCGCTATCAGTATCATATACACAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                      <  1:927306/147‑1 (MQ=255)
                                                                                  cgcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                  >  2:360374/1‑150 (MQ=255)
                                                                                   gcgcgTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                                  >  2:1004444/1‑149 (MQ=255)
                                                                                    cgcgTCCGGGACCCGGGGTGTACCTCAACCCGTAAACCCCCTAAAAGGATCATATAACCAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                                 <  2:278397/149‑1 (MQ=255)
                                                                                      ggtcccGGGACCGGTGGTGTACCCCAAACCGTAAACCCTCTATCAGGATAATATGACCACCCCCCCCGCCGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                                <  2:566977/145‑1 (MQ=255)
                                                                                      cgTCCCGGAACCGTGGTGTCTCCCAACCCGTAAACCCGCTACAAGAATAATATAACCACCCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                                 <  2:1340057/147‑1 (MQ=255)
                                                                                           ccgAACCGGTGGTTTATCCCGAACTTTAAACCCGTATTAAGAATAATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                         <  2:47087/149‑1 (MQ=255)
                                                                                           cGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                        >  1:474435/1‑151 (MQ=255)
                                                                                           cGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                        >  1:436545/1‑151 (MQ=255)
                                                                                                 ctgtgttgTATCTCAACCTTTAAACCCGTTATCAGAATAATATAACCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTgg                                                    <  2:1243189/148‑1 (MQ=255)
                                                                                                 cGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcg                                                  >  1:1278325/1‑151 (MQ=255)
                                                                                                    tggtgTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                                 >  2:1063546/1‑150 (MQ=255)
                                                                                                    tggtgTACCTCAACCCTTAAACCCGTTATAAGGATAATATAACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                      <  2:1022003/145‑1 (MQ=255)
                                                                                                            tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                         >  1:759637/1‑151 (MQ=255)
                                                                                                            tctcGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTg                                         >  1:1274907/1‑151 (MQ=255)
                                                                                                               ccaaCCCTTAAACCCGTTATCAGAATCATAAACCAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                             <  2:787845/143‑1 (MQ=255)
                                                                                                                 acccGTTAACCCGCTATCGGGATCATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                     <  1:435022/148‑1 (MQ=255)
                                                                                                                  aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                      >  1:567815/1‑148 (MQ=255)
                                                                                                                  aCCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                   >  2:1241854/1‑151 (MQ=255)
                                                                                                                   cctttAACCCTCTATCAGAATAATATAATAACCCCCCGCCCGGGTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                     <  2:951239/145‑1 (MQ=255)
                                                                                                                    cGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                                  >  2:1241197/1‑150 (MQ=255)
                                                                                                                         aaCCGCCTATAAGGATCATATAATAACCACCCGCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                             <  1:382750/135‑1 (MQ=255)
                                                                                                                           acccccTATCAGGATCATGTGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                               <  1:1186858/143‑1 (MQ=255)
                                                                                                                                   aTCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGTGCAGGCGCTGCATGACAACCGCGACCTGc                      >  1:760316/1‑150 (MQ=255)
                                                                                                                                      cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                   >  2:513261/1‑151 (MQ=255)
                                                                                                                                        ggATCATATGAACAACCCCCCCGCCGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                        <  1:856671/144‑1 (MQ=255)
                                                                                                                                          aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                >  1:230109/1‑150 (MQ=255)
                                                                                                                                          aTCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTc                >  1:167037/1‑150 (MQ=255)
                                                                                                                                            tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt              >  2:382501/1‑151 (MQ=255)
                                                                                                                                            tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt              >  1:1241098/1‑151 (MQ=255)
                                                                                                                                             cATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg             >  1:1340862/1‑151 (MQ=255)
                                                                                                                                               tatGACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                   <  2:543839/143‑1 (MQ=255)
                                                                                                                                                  tGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGc          >  2:1119291/1‑150 (MQ=255)
                                                                                                                                                    aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat        >  2:107953/1‑150 (MQ=255)
                                                                                                                                                      cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt     >  1:1206697/1‑151 (MQ=255)
                                                                                                                                                        acacCCTCGCCGATGCCCGCCGTTACCTCCTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTg    >  2:519903/1‑150 (MQ=255)
                                                                                                                                                          cacCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCg  >  1:1046696/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                                     
CCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACCGTGGTGTATCTCGAACCGTTAAACCGCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATTGCG  >  NC_000913/1376361‑1376654

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: