Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC42.9% 92.9 / 18.9 63S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (33/3);  new base C (0/27);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 8.29e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTG  >  NC_000913/3556062‑3556347
                                                                                                                                                 |                                                                                                                                            
gCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCg                                                                                                                                         >  2:502484/1‑151 (MQ=255)
    ccaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAg                                                                                                                                     <  2:595708/151‑1 (MQ=255)
       ccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                  >  1:841254/1‑151 (MQ=255)
        cTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAgtgt                                                                                                                                  >  2:71470/1‑150 (MQ=255)
                tGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTc                                                                                                                         >  1:618146/1‑151 (MQ=255)
                 gACGCGTTTTTCCCCGACGAAAAAAAAGCGTTCGGTGATTTTTTCCCTTTCGCCTCAAGGCGAGGCGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTTTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCa                                                                                                                        <  1:924067/151‑1 (MQ=255)
                       cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                   >  2:772895/1‑150 (MQ=255)
                        ttttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTc                                                                                                                   >  2:453148/1‑149 (MQ=255)
                            cACCAACGAAAAAAAAGCGTTGGGTTATATTTTCACTTTCGACTGAAAGCGAGACGGTTTCCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCAGGCGGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAAc                                                                                                             <  2:686072/151‑1 (MQ=255)
                                 acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                        >  2:724540/1‑151 (MQ=255)
                                   gacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                      >  2:204220/1‑151 (MQ=255)
                                      aaaaaaGCGTTCGGTAATATTTCCCCTTCCGACTAAAAGGGAGACGTTTTTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                   <  2:835436/151‑1 (MQ=255)
                                        aaaaGCGTCGGGTGTTTTTTTCACTTTCTACTGCAAGCGGGCCGGTATTCCCCGCCCTGTGGCGGCGGGGCGGGTTAGGAATATTCTGTTCAGGGAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                 <  2:349006/151‑1 (MQ=255)
                                          aaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                               >  2:55831/1‑151 (MQ=255)
                                          aaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCacac                                                                                               >  1:332971/1‑151 (MQ=255)
                                            gCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                              >  2:199276/1‑150 (MQ=255)
                                               ttCGGTGATATTTTCACTTTCGCCTTCAAGCGAGAGGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGc                                                                                           <  1:795197/150‑1 (MQ=255)
                                                taggTGATATTTTCACTTTCGACTTCAAGCGAGCCGGTATTATCCGGCCACTTCTCCCGCGGCAGGTTATGTATATTCTCGCCGCCAAGCGAAAAACTCCCCGCCAGTGTAGCGATTTTACGCCCAGTAACATGACGCGGCAGAcc                                                                                              >  2:159977/3‑146 (MQ=255)
                                                tCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                          >  1:835169/1‑150 (MQ=255)
                                                    taattttttCCCTTCCGCCTTAAAGGGAGAGGGTATTCCCCGCCCCTTGCCCGCGGGGCGGTTTATAAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                        <  2:586804/143‑1 (MQ=255)
                                                       tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                  >  1:682614/1‑151 (MQ=255)
                                                         ttttCCCTTTCGCCTACAAGGGAGCCGTTATTCCCCGGCCCCGGGGCGCGGGGCGGTTTAGGAATTTTCGTTTCAGGCAGGGAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACt                                                                                 <  2:248434/150‑1 (MQ=255)
                                                         ttttCCCTTTCGACTTAAAGCGAGACGTTATCACCCGCCCCCTGGTCGCGCGGCAGTTTATGAATTTTCTTTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATa                                                                                      <  2:466824/145‑1 (MQ=255)
                                                            tCACTTTCGACTTCAAGCGAGACGGAATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTAATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtctc                                                                             >  2:641692/1‑151 (MQ=255)
                                                             cccTTTCGCCTTCAAGCGAGACGTTTTTCCCCGCCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtct                                                                              <  2:99107/147‑1 (MQ=255)
                                                              ccTTTCGCCTTAAAGCGAGCCGGTTTTCCCCGGCCCCGGGTCGCGGGGCAGTTTATGAATTTTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                           <  2:477983/150‑1 (MQ=255)
                                                               cTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                          >  1:449595/1‑151 (MQ=255)
                                                               cTTTCACCTTCAGGCGAGCCGTTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATAAATATTCTGTTCATGCAGGGAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtctc                                                                             <  2:877517/148‑1 (MQ=255)
                                                                   cGACTTAAAGCGAGACGTTATTCCCCGCCCCCGGGCCGCGCGGCAGTTTATGAATTTTCTGTTCATGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                      <  2:723993/151‑1 (MQ=255)
                                                                     aCTTAAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACTAGCTAATAGAACTTCTCCACGCATCt                                                                    <  1:391476/151‑1 (MQ=255)
                                                                       ttCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                    <  1:153952/149‑1 (MQ=255)
                                                                        tCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                  >  1:418567/1‑150 (MQ=255)
                                                                          aaGCGAGCCGGTATTCCCCGGCCCCTGGGCGCGCGGCAGGTAATAAATTTTCTTTTCTTGAGGGGAAAAGCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCa                                                                       <  1:737993/143‑1 (MQ=255)
                                                                            gCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                             >  1:704844/1‑151 (MQ=255)
                                                                             cGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATg                                                            >  2:79879/1‑151 (MQ=255)
                                                                              gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                           >  2:656615/1‑151 (MQ=255)
                                                                              gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                           >  1:12920/1‑151 (MQ=255)
                                                                                     tATTCCCCGGCCCCTGGCCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                     <  1:219474/150‑1 (MQ=255)
                                                                                          cccGGCCCCTGGTCGCGCGGCAGTTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                               <  1:555893/151‑1 (MQ=255)
                                                                                          cccGGCCCCTGGCCGCGCGGGAGTTTATGATTATTTTGTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                 <  2:654961/149‑1 (MQ=255)
                                                                                          cccGGCCCCGGGCCGCGCGGCAGTTTATAAATATTCTGTTAAGGAAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTc                                                 <  2:246665/149‑1 (MQ=255)
                                                                                           ccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                              >  2:897200/1‑151 (MQ=255)
                                                                                               cccgTGGCCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAAACCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCaa                                            <  1:833753/149‑1 (MQ=255)
                                                                                               ccccTGGTCGCGCGGCAGGTTATGAATATTCTTTTCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                          <  1:891278/151‑1 (MQ=255)
                                                                                                    ggccgcgcgGCAGTTTATAAATTTTCTGTTCATGAAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTg                                        <  2:564832/145‑1 (MQ=255)
                                                                                                    ggTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                      <  2:683616/150‑1 (MQ=255)
                                                                                                     gTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                      >  1:167025/1‑149 (MQ=255)
                                                                                                            ggCAGTTTATAAATATTCTGTTCATGCAGGGAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCCCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGtt                               <  2:427495/149‑1 (MQ=255)
                                                                                                            ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                              >  2:367346/1‑150 (MQ=255)
                                                                                                               aGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACg                          >  1:779857/1‑151 (MQ=255)
                                                                                                                   tATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                      >  2:761681/1‑151 (MQ=255)
                                                                                                                    aTGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                     >  2:577727/1‑151 (MQ=255)
                                                                                                                         ttttCTGTTCATGCAGGAAAAACCTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                 <  1:950259/148‑1 (MQ=255)
                                                                                                                           ttCTGTTCAGGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                 <  1:397657/148‑1 (MQ=255)
                                                                                                                             cTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg              >  2:476634/1‑149 (MQ=255)
                                                                                                                                ttCATGCAGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc          <  1:139238/150‑1 (MQ=255)
                                                                                                                                 tCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACtt        >  2:596482/1‑151 (MQ=255)
                                                                                                                                  cATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc       >  1:741006/1‑151 (MQ=255)
                                                                                                                                  cATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTCCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc       >  2:851625/1‑151 (MQ=255)
                                                                                                                                      cAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGtt   >  1:939113/1‑151 (MQ=255)
                                                                                                                                       aGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg  >  1:356206/1‑151 (MQ=255)
                                                                                                                                        ggaaaaacCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg  <  1:326789/147‑1 (MQ=255)
                                                                                                                                        ggaaaaacCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg  <  1:863736/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTG  >  NC_000913/3556062‑3556347

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