Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,622,2500AG38.5% 56.4 / 10.6 39E161G (GAG→GGG) yjjJtoxin YjjJ
Reads supporting (aligned to +/- strand):  ref base A (4/20);  new base G (15/0);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.54e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATTG  >  NC_000913/4622016‑4622498
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
ggCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGc                                                                                                                                                                                                                                              <  2:135261/249‑1 (MQ=255)
         gAAGTTGTCTGGTAAAAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGCCCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGCTGGTTGGTTGGGGGGGGGaagtttaagcgctgggtttct                                                                                                                                                                                                                                   >  2:234857/1‑232 (MQ=255)
                                 aTGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCct                                                                                                                                                                                                           <  2:68406/251‑1 (MQ=255)
                                 aTGGCGACGAACAGTGGTTGGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATGCTGGCGGTTGGTTGGTCGGGGAGGGGAATGATCAGCGATGGATTACTGCACAACCCCCTTCTGTAATTCct                                                                                                                                                                                                           <  2:267449/251‑1 (MQ=255)
                                    gCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTTTTAGCGATGGGTTACTGCACAACACCCCTCGGAAATTCctctt                                                                                                                                                                                                        >  1:352893/1‑250 (MQ=255)
                                            cAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTTGGCGGGGGGGGGAAATTTCCGCGGTGGATTACTGCCCAACACCCCGCGGGAATTCCTCTGGATCaaa                                                                                                                                                                                                 >  1:412/1‑250 (MQ=255)
                                                    tGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc                                                                                                                                                                                         <  1:442745/250‑1 (MQ=255)
                                                      aTGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGAAATTCCTCTTGATCAAAAACTCccccct                                                                                                                                                                                      >  1:139174/1‑249 (MQ=255)
                                                       tGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGGCGGGGGGGGGGATTATCAGCGATGGATTACTTCACAACACCCTGCGGGAATTTCTCTTGATCAAAAACACACCCAtt                                                                                                                                                                                     >  1:117637/1‑251 (MQ=255)
                                                                                  ccGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGc                                                                                                                                                              <  2:46346/246‑1 (MQ=255)
                                                                                  ccGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTg                                                                                                                                                          <  1:316382/250‑1 (MQ=255)
                                                                                          aGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATTAGCAGTGGCTTACTCTACAAAACCCTCCGGAAAATCCTTCGGGTCAAAAAATCCCCCATTACCAACAAA‑TGGCAAGTgcagccccggcaggggg                                                                                                                                                  >  2:364112/1‑234 (MQ=255)
                                                                                            ggCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTTTCAGCGATGGAATAATGCCCAACACCCCTCGGGAATTCCTCTGGATCAAAAACCAACCCCTCACGAACAGC‑TGGCAAATGATGAACTGGCCGGCGaaa                                                                                                                                                >  1:191511/1‑250 (MQ=255)
                                                                                                ttttAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATTCGCGGTGGGTTACCGCACAAAACCCCGCGGGAAATCCTCTTGGGCAAAAAATCACCAATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGGAAAtgtg                                                                                                                                           >  1:14717/1‑251 (MQ=255)
                                                                                                   tAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCAATGGCAGGAGAAATTTTgggg                                                                                                                                        >  1:360783/1‑250 (MQ=255)
                                                                                                         gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGAAATATCAGCGGTGGATTACTGCACAACACCCTGCGGGAATTTCTCTGGATCAAAAAATCACCCATTACGAACAGCATGGCA‑GTGGTGCACTGGCAGGGGAAATTGTGGGTtcttct                                                                                                                                  >  2:9526/1‑251 (MQ=255)
                                                                                                               ggggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgggggc                                                                                                                            <  2:353326/251‑1 (MQ=255)
                                                                                                               ggggCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTC‑ACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGGTGGATTTATGCACAACACCCCGCGGGAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGGGATGCACTGGCCGGAGGAATTGTGGGTTCTTCTgcgggcg                                                                                                                           >  1:373341/1‑251 (MQ=255)
                                                                                                                            aGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACCCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑a                                                                                                               <  2:236826/250‑1 (MQ=255)
                                                                                                                                 cgcAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGGATTATCAGCGATGGGTTACTGCACAAAACCCCGCGGAAATTCCTCTGGGTCAAAAAATCACCCATTACGAACAGC‑TGGCAAGTGGTGCACTGGCAGGAGAAATTGTGGGGTCTTCTGCGGGCGGCGGGCAGCC‑CAAAttt                                                                                                         >  1:156711/1‑251 (MQ=255)
                                                                                                                                                  gatgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACaaa                                                                                         <  1:293836/250‑1 (MQ=255)
                                                                                                                                                     gatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTTCCTACTATGCACAAACGc                                                                                      <  2:307899/250‑1 (MQ=255)
                                                                                                                                                      atatACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCg                                                                                    <  2:56257/251‑1 (MQ=255)
                                                                                                                                                            gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTg                                                                                                                                                                  <  2:151737/168‑1 (MQ=255)
                                                                                                                                                            gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTg                                                                                                                                                                  >  1:151737/1‑168 (MQ=255)
                                                                                                                                                             tctcTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTGTTGGTCCGGGGGGGGAATTCTCACCAAACGATTACTGCACAACACCCTGCGGGAAATCCTCTGGGTCAAAAACTCACCCCATACGACCAGC‑CGGCAAATGTTGCACTGGCAGGCGAAAATGTTGGTTCTTTTGCGGGCGGGGAGGAGCC‑CAAAACTAACTAATATGCACAAACCACGTCAgga                                                                              >  2:88363/1‑249 (MQ=255)
                                                                                                                                                                 tGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACCCCCCGCGGAAATTCCTCTGGATCAAAAACTCACCCCTTACGAAAAGC‑TGGCAAGTGATGCCCCGGCAGGAGAAAATGTGGGGTCTTTTGCGGGGGGCGGGCCGCCCAAAATTAACC‑ACTATGTACAAACCCCGCCAgggcaaa                                                                          >  2:173372/1‑245 (MQ=255)
                                                                                                                                                                  ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATa                                                                          <  1:181959/249‑1 (MQ=255)
                                                                                                                                                                    cAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGa                                                                                                                                                                                                                                            <  1:287064/87‑1 (MQ=255)
                                                                                                                                                                             gaTGTGCTCTACGCCCTGACTGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACGCACCCATTACGAACAGC‑TGGCAAGTGAGGCACTGGCAGGAGAAATTGTGGGTTCGTCTGCGGGCGGCGAGCAGCC‑AAAATTTACGTACTATTCACAAACGCCGGCAGGCAATAAGCATGTGtttgt                                                             <  2:310824/251‑4 (MQ=255)
                                                                                                                                                                                gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGGTGGTCGGGGGGGGGAATTATCAGCGATGGGTTACTGCACAACACCCCGCGGGAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAAATGTGGGGTCTTCTGCGGGCGGCGGGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTgggga                                                           >  1:201024/1‑247 (MQ=255)
                                                                                                                                                                                           ccTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTAc                                               <  1:294491/251‑1 (MQ=255)
                                                                                                                                                                                              gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTAcc                                              <  1:301769/249‑1 (MQ=255)
                                                                                                                                                                                              gACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACca                                             <  1:336553/250‑1 (MQ=255)
                                                                                                                                                                                               aCCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACcaca                                           <  1:242223/251‑1 (MQ=255)
                                                                                                                                                                                                       tAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGGGGGGGATTATCAGCGATGGATTACTGCACAACACCCCGCGGAAATTTCTCTGGGTCAAAAACTCACCCCTTACGAACAGC‑TGGCAAGTGATGCACCGGCAGGGGAAATTTTGGGTTCTTTTGCGGGCGGGGAGGAGCC‑AAAATTTTCCTACTATGCACAAAAGCCGTCAGGCAATAAAAATGTGTTGGTGGAATTTACCGTACCACAGCAAAccc                                   >  1:171927/1‑250 (MQ=255)
                                                                                                                                                                                                        aaCGGTGAATATACTGGCGGTGGGTTGGTCGTGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAAGTCCTCTGGATCAAAATCGCACCCATTACTAAGATC‑TGGCAAGTGATGCGTTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGTGGCGATCAGCT‑AAAATTTACCTACTATGCACAAACTCCGGCATGCAATAGTCATGTGTTGGTGAATTTCACCGTACCACAGCAAACcg                                   <  1:443120/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggttggttGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTgggg                >  1:247633/1‑251 (MQ=255)
                                                                                                                                                                                                                            ttggttggTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGCTCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTg              >  1:447200/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGCCCAAAAACTCACCCATTACGAACAGC‑TGGCAACTAATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGAGGGCAGCGAGCAGCC‑AACATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGACAACGTTGGGGTGACCTGCTAAt    >  2:431891/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATTg  >  1:410595/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggaggggaATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATTg  >  1:348546/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
GGCCGCAGGGAAGTTGTCTGGTAACAGGCGCAGATGGCGACGAACAGTGGTTTGATGGTTTGCCCTGGTATTTGACCGATCTCCGACCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTGGTCGGGGAGGGGAATTATCAGCGATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGC‑TGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCC‑AAAATTTACCTACTATGCACAAACGCCGTCAGGCAATAAACATGTGTTGGTGAAATTCACCGTACCACAGCAAACCGCGGTCAGCCAACGTTGGGGTGACCTGCTAATTG  >  NC_000913/4622016‑4622498

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: ‑