Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2000GA45.6% 95.1 / 52.4 93S217F (TCC→TTC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (40/10);  new base A (0/42);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 5.51e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556051‑3556343
                                                                                                                                                            |                                                                                                                                               
ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGa                                                                                                                                                <  1:521908/151‑1 (MQ=255)
     tgtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAAc                                                                                                                                           >  1:316055/1‑151 (MQ=255)
      gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACt                                                                                                                                          <  2:110324/151‑1 (MQ=255)
      gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACt                                                                                                                                          <  1:217201/151‑1 (MQ=255)
      gtgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACt                                                                                                                                          >  1:1312740/1‑151 (MQ=255)
        gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTcc                                                                                                                                        >  2:1262653/1‑151 (MQ=255)
        gCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTc                                                                                                                                         <  1:711049/150‑1 (MQ=255)
         ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTccc                                                                                                                                       <  1:220912/151‑1 (MQ=255)
            cAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGc                                                                                                                                    >  1:1082582/1‑151 (MQ=255)
              acccccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGcc                                                                                                                                   <  2:314239/150‑1 (MQ=255)
                 accTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGCGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAg                                                                                                                                 >  2:1226314/1‑149 (MQ=255)
                    tCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                              >  2:246688/1‑149 (MQ=255)
                                cttttttCCCCGCCAAAAAAAAAGCGTTGGGTGATATTTTCTCTTTCAACTCAAAGCGAGGCGGTATCCCCCGGCCCCTGGCGGCGCGGAGGGTTAAGAATTTTTTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACg                                                                                                                  <  1:1236951/146‑1 (MQ=255)
                                cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCATGc                                                                                                                 >  1:1300234/1‑150 (MQ=255)
                                 gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                               >  2:194475/1‑151 (MQ=255)
                                   tttcccccGACAAAAAAAAGGGTTTCGTTAATTTTTTCCTTTTCGACTTAAGGCGAGACGTTATTCCCCGGCCCCTGTTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCGGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                               <  1:736380/143‑1 (MQ=255)
                                    ttttccccaagaaaaaaaaGTGTTCGGGTGATATTTCCCCTTTCGCCTTCAACCGAGCCGTTATTCCCGGCCCCCTGGCCGCGCGGCAGGTTAAAAAAATTTTTTTCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                              <  1:146980/137‑1 (MQ=255)
                                    tttCCCCGACGAAAAAAAAGCGTTGGGTAATTTTTTCCTTTTCTACTCAAAGCGAGAGGGTTTTCCCCGGCCCCTGGTCGCGGGGCAGTTTATGAATATTTTTTTCATGAAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAg                                                                                                             <  1:59341/150‑1 (MQ=255)
                                    tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                            >  2:868161/1‑151 (MQ=255)
                                       cACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                         >  2:854480/1‑151 (MQ=255)
                                         ccGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACAt                                                                                                       >  1:191761/1‑151 (MQ=255)
                                         ccGACAAAAAAAAGGCGTGCGGTGATTTTTCCACTTCCGACTTCAAGCGAGCCGTTATTCCCCGCCCCTTGGTCGCGGGGCAGTTTATAAATATTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACa                                                                                                        <  1:869360/150‑1 (MQ=255)
                                          cgacaaAAAAAAAGGGTTGGTTGATTTTTTCGCTTCCGCCTAAAACAGAGAGGTTATTCCCCGCCCGCGGGGCGCGCGGAAGTTAATAAATATTTTGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                      <  2:434108/151‑1 (MQ=255)
                                              gacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                  >  1:648531/1‑151 (MQ=255)
                                               acAAAAAGGCGTTGGGTAATTTTTCCCCTTCCGCCTTCAAGCGAGCCGGTATCCCCCGCCCCCGGGTCGCGGGGCAGTTTATGAATTTTCTGTTCAGGAAGAGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                     <  2:151643/147‑1 (MQ=255)
                                                cAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                >  1:17437/1‑151 (MQ=255)
                                                 aaaaaaGTGTCGGGTAATATTTTCCCTTCCGACTCAAAGGGAGCGGGTATTCCCCGGCCCCGGGTCGCCCGGCAGGTTATGAATTTTCTGTTCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  1:255222/151‑1 (MQ=255)
                                                 aaaaaaGCTTTCGGTATTTTTTTCACTTGCGGCTGAAGAGGGGACGGTATTCCCCGCCCCCGGGCCGCGGGCGGGTTTTTGATTTTTCTGTTCACGGGGAAAAAAACACCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                 <  1:539026/149‑1 (MQ=255)
                                                 aaaaaaGCGTTCGTTGTTTTTTCCCCTTTCGACTAAAAGGGAGACGGTATTCCCCGCCCCTGGTTCCCGCGGCAGGTTAGAAATATTCGGTTCAGGCAGGAAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  2:1147321/151‑1 (MQ=255)
                                                 aaaaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               >  2:140756/1‑151 (MQ=255)
                                                 aaaaaaGCGTTCGGGAATATTTTCACTTTCGCCTCCAAGCGAGACGGTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  1:1064369/151‑1 (MQ=255)
                                                 aaaaaaGCGTCGGGTGATTTTTTCCCTTCCGTCTACAAGGGGGACGGTATTCCCCGCCCCCTGCCCGCGGGGCAGTTTATGAATTTTCTGTCCAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                               <  1:426650/151‑1 (MQ=255)
                                                  aaaaaGCTTTCGGTGATATTTTCTTTTCCGACTTCAAGCGAGCCGGTATCCCCCGCCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                              <  1:271383/151‑1 (MQ=255)
                                                   aaaaGCGTTGGGTGATTTTTCCCCTTTCGCCTTAAAGCGAGACGGTATTCCCCGCCCCGTGGTCGCGCGGAAGTTTATAAATTTTCTGTTAATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                             <  1:1494724/151‑1 (MQ=255)
                                                    aaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                             >  1:684368/1‑150 (MQ=255)
                                                     aaGCGTTCGGTGATTTTTTCCCTTTCGACTTAAAGCGAGACGGTTTTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTAATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                           <  1:480563/151‑1 (MQ=255)
                                                          ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCAAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAg                                                                                        >  2:84309/1‑149 (MQ=255)
                                                              ggTGATTTTTTCATTTTCGACTTCAAGCGAGACGTTATTCCCCGGCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                      <  1:1501601/148‑1 (MQ=255)
                                                               gTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                   >  1:1278163/1‑150 (MQ=255)
                                                                 gATTTTTTCCCTTTCGCCTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                  <  1:1159522/149‑1 (MQ=255)
                                                                  atatTTTCTCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCGGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                <  1:1183612/150‑1 (MQ=255)
                                                                   ttttttCTCTTCCGCCTACAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGa                                                                                <  1:70533/147‑1 (MQ=255)
                                                                   ttttttCCCTTTCGACTTAAAGGGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTTTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                 <  1:889041/146‑1 (MQ=255)
                                                                   ttttttCACTTCCGACTTCAAGCGAGACGGTTTTCCCCGGCCCCTGGCCGCGGGGCAGTTTAGGAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTa                                                                                     <  2:698479/142‑1 (MQ=255)
                                                                   tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                              >  2:1380111/1‑151 (MQ=255)
                                                                     ttttCCTTTTCGCTTTCAAGCGGGACGTTTTCCCCCGGCCCCTGGCCGCGCGGCGGTTTATAAATTTTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAg                                                                                 <  2:1005951/146‑1 (MQ=255)
                                                                     ttttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                            >  1:214451/1‑151 (MQ=255)
                                                                      tttCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                           >  2:964792/1‑151 (MQ=255)
                                                                         cACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                        >  1:1157525/1‑151 (MQ=255)
                                                                          aCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTcc                                                                        >  1:463437/1‑150 (MQ=255)
                                                                           cTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                      >  1:867618/1‑151 (MQ=255)
                                                                            tttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                     >  2:558435/1‑151 (MQ=255)
                                                                             ttCGCCTTCAAGCGAGGCGGTATTCCCCGCCCCCGGGTCGCGGGGCAGGTTATGAATTTTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                     <  1:457597/150‑1 (MQ=255)
                                                                             ttCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATTTTCTGTTAATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACg                                                                     <  1:1471712/150‑1 (MQ=255)
                                                                               cGACTCCAAGCGAGACGGTTTTCCCCGCCCCCTGGCCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                   <  1:295877/150‑1 (MQ=255)
                                                                                 aCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                >  2:844646/1‑151 (MQ=255)
                                                                                 aCTTAAAGCGAGGCGGTATTCCCCGGCCCCTGGCCGCGCGGCGGTTTATGAATTTTCTGTTCAGGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                      <  1:374326/145‑1 (MQ=255)
                                                                                  cTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                               >  1:598991/1‑151 (MQ=255)
                                                                                   ttaaaGCGAGCCGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                <  1:194664/146‑1 (MQ=255)
                                                                                   tcaaaGCGAGCGGTTTTTCCCGGGCCCCTGGTCGCGCGGCGGTTTATAAATTTTCTTTCCAGGCGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                   <  1:798265/143‑1 (MQ=255)
                                                                                   taaaggCGAGACGGTATTCCCGGGCCCCTGGCCGCGGGCAAGTTTATAAATTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                <  2:1356657/144‑1 (MQ=255)
                                                                                     cAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                            >  2:1424216/1‑151 (MQ=255)
                                                                                      aggCGAGACTGTATTCCCCGCCCCCTGGTCGCGGGGCAGTTTATAAATTTTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                            <  1:618276/148‑1 (MQ=255)
                                                                                      aaGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTTTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                           <  1:1000512/151‑1 (MQ=255)
                                                                                      aaGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                            >  2:1258614/1‑150 (MQ=255)
                                                                                          gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                       >  1:225779/1‑151 (MQ=255)
                                                                                           agaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                      >  1:943573/1‑151 (MQ=255)
                                                                                             gccgatttccccGGCCCTCTGCGCGCGCGGCGGTTGTTAAATTTTGTGTTCATGCGGCAAAACCCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                        <  2:531806/139‑1 (MQ=255)
                                                                                               ggTATCCCCCGGCCCCGGGTCGCGCGGCAGTTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  <  2:744205/151‑1 (MQ=255)
                                                                                                  tttgcccgggCCCTTGGCCGCGCGTCGTGTTAGAATTTTCCTGTTCATGCAGGCAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                   <  2:48184/140‑1 (MQ=255)
                                                                                                   tcccccGGCCCCGGGTCGCGCGGGAGTTTATAAATTTTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  <  2:1494489/144‑1 (MQ=255)
                                                                                                     ccccGCCCCCGGGCCGCGGGGAGGTTTATAAATTTTCTTTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                  <  2:1089043/144‑1 (MQ=255)
                                                                                                      cccGGCCCCTGGCCGCGCGGCAGGTTATAAATATTTGTTCCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  2:889506/149‑1 (MQ=255)
                                                                                                      cccGGCCCCGGGTCGCGCGGAAGTTTATAAATATTTTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  2:273097/149‑1 (MQ=255)
                                                                                                      cccGCCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                             <  1:1415862/149‑1 (MQ=255)
                                                                                                             cccTGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                     <  1:1262671/151‑1 (MQ=255)
                                                                                                                       gcgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                             >  2:1224937/1‑151 (MQ=255)
                                                                                                                            ggcggtttATAAATATTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                         <  2:1503162/145‑1 (MQ=255)
                                                                                                                            gCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                         >  1:311862/1‑150 (MQ=255)
                                                                                                                            gCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaa                         >  2:1482981/1‑150 (MQ=255)
                                                                                                                               ggTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                       <  1:868177/149‑1 (MQ=255)
                                                                                                                                  ttATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                   >  2:340545/1‑151 (MQ=255)
                                                                                                                                  ttATGAATATTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                     <  1:209819/149‑1 (MQ=255)
                                                                                                                                   tATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  2:574375/1‑151 (MQ=255)
                                                                                                                                   tATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  2:595155/1‑151 (MQ=255)
                                                                                                                                   tATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                  >  1:515550/1‑151 (MQ=255)
                                                                                                                                    aTGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                 >  1:609764/1‑151 (MQ=255)
                                                                                                                                       aaTTTTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                <  2:463431/149‑1 (MQ=255)
                                                                                                                                       aaTATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg              >  1:647125/1‑151 (MQ=255)
                                                                                                                                          atTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc           >  2:1061257/1‑151 (MQ=255)
                                                                                                                                            tCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg          <  1:1136124/150‑1 (MQ=255)
                                                                                                                                               tGTTCAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCgag        <  2:1275822/150‑1 (MQ=255)
                                                                                                                                                       gCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc  <  1:194084/149‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                               
CCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCC  >  NC_000913/3556051‑3556343

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: