Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,092,703 A→G F41L (TTC→CTC)  frvA ← putative PTS enzyme IIA component FrvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7030AG97.7% 148.9 / NA 43F41L (TTC→CTC) frvAputative PTS enzyme IIA component FrvA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/19);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGC  >  NC_000913/4092457‑4092940
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
gAATAATTTTGCGACACAGTGTGCCGATGATTTTGGCCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTCTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTgg                                                                                                                                                                                                                                           <  1:21532/251‑1 (MQ=255)
          gCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCCTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCgg                                                                                                                                                                                                                                 <  1:42097/251‑1 (MQ=255)
             acacAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTg                                                                                                                                                                                                                               >  1:105779/1‑250 (MQ=255)
               acaGTGGGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATa                                                                                                                                                                                                                            <  1:537916/251‑1 (MQ=255)
                                        aTCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCg                                                                                                                                                                                                   >  1:193806/1‑251 (MQ=255)
                                          cctcctCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACAGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGcc                                                                                                                                                                                                 <  2:528176/251‑1 (MQ=255)
                                               tCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACt                                                                                                                                                                                             <  2:156051/250‑1 (MQ=255)
                                                 gccgcTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGTCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTAATGTCGCCAACTGCt                                                                                                                                                                                          <  1:361670/251‑1 (MQ=255)
                                                                   gAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAgcgc                                                                                                                                                                        >  1:521463/1‑251 (MQ=255)
                                                                    aGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGcc                                                                                                                                                                       <  1:272990/251‑1 (MQ=255)
                                                                             cAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                              >  2:318663/1‑251 (MQ=255)
                                                                                  aTTGACGTCTTCGCCATCGCTGGCTTTCACGTCACTAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGCGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGCACGGAATAAGCGCCATTGCCCTGAatat                                                                                                                                                         >  2:89305/1‑251 (MQ=255)
                                                                                      aCGTCTTCGCCATCGCTGACTTTCCAGTCAATAGCCTGCGCTTTGAGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGAGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGCGTCGGTGATAAAACAGTTTTGTAACGCTATTGTCGCCAACTGCTTCAAAACGGAATAAGCGCCCTTGCCCTGAATATTc                                                                                                                                                       >  2:504163/1‑249 (MQ=255)
                                                                                              gCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATAc                                                                                                                                             <  2:106291/251‑1 (MQ=255)
                                                                                                cATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATAc                                                                                                                                             >  1:4084/1‑249 (MQ=255)
                                                                                                               agtcacgcgccagcgatttGCGGGCGTATAATACGAACGGTTGTTTAACGCCGGCGCTTTTACCGTGCGGAACGGCGACACCCGA‑CAAAACACCAGTAGAGCCCATTTTTTCGAGCAGCAGCCGGGTTTGCCGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAAAATTACGGTCAATACAACTTGCAGTacaggct                                                                                                                            >  2:67579/17‑245 (MQ=255)
                                                                                                                  cAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCTCGAAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTATCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCAtt                                                                                                                                                                    <  1:69854/208‑1 (MQ=255)
                                                                                                                  cAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAAAGCAGGAGCTTTTACCGTGCGGCACGGCGAAACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGAAGCAGCAGGGGTTGCAGGAGCTGGTGTGAGTCGGAGATAAAACCGTGTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCAta                                                                                                                                                                    >  2:69854/1‑207 (MQ=255)
                                                                                                                        ccTGCGCTTAGCGGGCGAATAATACGCACGGTTGTTTAACGCCGGCGCTTTTACCGTGAGGAACGGCGACACCGGCACAAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAACACCGTTTTGTAAAGCTATTGTCGCCAACTGCTTCAGAACGGAATAAACGCCATTGCCCTAAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTact                                                                                                                   >  2:231437/1‑249 (MQ=255)
                                                                                                                            cgcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCgg                                                                                                                >  1:91071/1‑250 (MQ=255)
                                                                                                                             gcTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGAc                                                                                                              >  1:300196/1‑251 (MQ=255)
                                                                                                                                  gCGGGCGAATAATACGAACTGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCAt                                                                                                         <  2:193806/251‑1 (MQ=255)
                                                                                                                                         aataatACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACG‑TATTGTCGCCAACCGCTTCCGGAAGGAATAAACCCCATTGGCCTGAATATTCAGGTCAATAAAACTTGCAGTAAGAGCTGCCATTAATTGGaattcctaaaatta                                                                                                  >  2:144417/1‑237 (MQ=255)
                                                                                                                                                  aaCGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGCGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATc                                                                                         >  1:180070/1‑251 (MQ=255)
                                                                                                                                                     ggTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAt                                                                                      >  1:106023/1‑251 (MQ=255)
                                                                                                                                                            tAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAt                                                                                <  1:115739/250‑1 (MQ=255)
                                                                                                                                                             aaCGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtg                                                                               <  2:300196/250‑1 (MQ=255)
                                                                                                                                                              aCGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggt                                                                             >  2:379567/1‑251 (MQ=255)
                                                                                                                                                                 cAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCACA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggtgg                                                                           >  1:197829/1‑250 (MQ=255)
                                                                                                                                                                 cAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCATTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGa                                                                          <  2:197829/251‑1 (MQ=255)
                                                                                                                                                                 cAGGCACTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggtgg                                                                           >  2:440458/1‑250 (MQ=255)
                                                                                                                                                                  aGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGaa                                                                                                                                                                                    <  2:151818/144‑1 (MQ=255)
                                                                                                                                                                  aGGCGCCTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGaa                                                                                                                                                                                    >  1:151818/1‑144 (MQ=255)
                                                                                                                                                                   ggCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGc                                                                        <  1:322655/251‑1 (MQ=255)
                                                                                                                                                                    gcgcTTTTACCGTGCGGCACGGCGACACCGGAACCAAC‑TCCAGTAGAGTGCATTTTTTCGCGCAGCCGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCCGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCg                                                                       >  1:348820/1‑251 (MQ=255)
                                                                                                                                                                                 gCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtgatga                                                           <  1:379567/250‑1 (MQ=255)
                                                                                                                                                                                                  gAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTa                                          >  2:517260/1‑250 (MQ=255)
                                                                                                                                                                                                       aaa‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGTCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                                    <  1:134286/251‑1 (MQ=255)
                                                                                                                                                                                                               gTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAtt                             <  2:105779/251‑1 (MQ=255)
                                                                                                                                                                                                                   agTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtga                                                              >  1:410885/1‑214 (MQ=255)
                                                                                                                                                                                                                   agTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtga                                                              <  2:410885/214‑1 (MQ=255)
                                                                                                                                                                                                                                        gcagggggtgaaGGAACTGGCGTGAGTCGGTGATACAACAGTTTTGGAACGGAATTGGAGCCGACTGATTCAGAACGGACTATGCGCCATTGCCCTGAATATTCAGGTCAATACAGATTGCAAAAAGAGCTACCTTTAATCGCACGCCCTAAGATATAATTTAATCAATTGCAATGGTGGAACGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATAAACATTTGAGACAGAGTGAAAGGTCAGAttt    >  2:468529/12‑251 (MQ=255)
                                                                                                                                                                                                                                          agGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACAGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGc  >  1:499196/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             
GAATAATTTTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTTTGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAA‑TCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGC  >  NC_000913/4092457‑4092940

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