Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,571,280 C→T E99K (GAG→AAG)  yjiR ← fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,571,2800CT97.6% 126.2 / NA 42E99K (GAG→AAG) yjiRfused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/31);  minor base G (0/1);  total (10/32)
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.99e-01

ATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTG  >  NC_000913/4571033‑4571526
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
aTCGTGTTAGGGCACCCGCTGGTAATGATGATGTCATCGGCGGTGCCCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAg                                                                                                                                                                                                                                                       <  2:528489/251‑1 (MQ=255)
aTCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAg                                                                                                                                                                                                                                                       <  1:253079/251‑1 (MQ=255)
      ttATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGGAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCCAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCCGGGCGCTTCGACATCCCGGGGGCTTTTTCTTAACGGAAAAATGGAAATGTCGCTTTGCGCCTTCCGCATATc                                                                                                                                                                                                                                                 >  1:88224/1‑251 (MQ=255)
        agggcggCCGCTGGTGATGATGATGTCATCGGGGGTGACCACCGAGCCGCTGTCGAGCACCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTGGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATcc                                                                                                                                                                                                                                                <  1:156321/244‑1 (MQ=255)
        agggCAGCCGCTGGAGATTATTATTTCATCGCGGGGGTCCCCCGAGCCGCTGCCGAGCATAAGGCGGTCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCAAAACCGAGACCGGTTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCa                                                                                                                                                                                                                                               <  1:320389/249‑1 (MQ=255)
         gggCAGCCGGTGGTGATTATGTTGTCTCCGGCGGTGCCCACCGAGCCGCTGTCGAGCATCAGGCGGCCAATCTGCTCCCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGTTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCTGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCa                                                                                                                                                                                                                                               <  2:426244/249‑1 (MQ=255)
          ggcagacgggggggaggaggatgTTATCGGGGGTGACCTCCGAGCAGCTGTCGGGCATCAGGGGGGAAATCGGCTCTGGCAATAGGCGTTGACGGGCTAACAAGGCATAACCGTGAACGGTTTGCGGATTATGCTGTACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTGAACGGAACAATGGAAAGGTCGCTATGCGCCTTCAGCTTATCCAGc                                                                                                                                                                                                                                             <  2:282697/232‑1 (MQ=255)
                cGCTGTTGATGAGGGTGTCAGCGCCGGTACCCCCCGAGCCGCGGTCGAGCATCAGTCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATACCCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATAGCCAGCACCTGa                                                                                                                                                                                                                                       <  2:390197/251‑1 (MQ=255)
                                   aTCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCa                                                                                                                                                                                                                    <  1:224819/251‑1 (MQ=255)
                                       gctgacacccccGAGCCGCTGTCGAGCATCAGGCGGCCAATCTGCGCACGCTATACGCGCTACCGGGCGAACAAGTCATACCCGAGACCGTTTTGCAGATTAGGCTGCACCCCCCGGCTTAGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCgg                                                                                                                                                                                                                <  2:482235/241‑1 (MQ=255)
                                            gACCCCCGAGCCGCTGTCGTGCATCCGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGAGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAGCAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCGGATCCCACTGGGTCATTTCCACCGGGCGCt                                                                                                                                                                                                           >  2:520480/1‑251 (MQ=255)
                                               caccGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGc                                                                                                                                                                                                         >  1:573912/1‑250 (MQ=255)
                                                             cGAGCACCAGGGGGGCAATCTGCTCCCGCAATACTCGCTGCCCGGCTAACAAGTCATACCCGAGAGGGGTTGGCAGATTATGCGGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCAt                                                                                                                                                                                           <  2:99336/250‑1 (MQ=255)
                                                                  aTCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACAGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGCCGCGTCATcggcgg                                                                                                                                                                                     >  2:1140/1‑251 (MQ=255)
                                                                       gCGGGCACTCTGCTCTCGCCATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGCACTg                                                                                                                                                                                >  2:241496/1‑251 (MQ=255)
                                                                                      ccgcaacccccGCTGACCGGCCAACAAGTCATAACCGAGAACGGTTTGCAGATTAGGCGGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTAc                                                                                                                                                               <  2:66165/241‑1 (MQ=255)
                                                                                             ccTCGCTGACCGGCTAACAAGCCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGc                                                                                                                                                        <  1:459246/250‑1 (MQ=255)
                                                                                                  cTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGGCTggcggacgggggtttcgggtgtcggcg                                                                                                                                                      >  2:131470/1‑228 (MQ=255)
                                                                                                             aaCAAGTCATAACCGAGAACGGTTTGCAGATTTGGCGGCCCCGCCCGGCTTGGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACAGCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAAccc                                                                                                                                           <  1:171799/249‑1 (MQ=255)
                                                                                                              aCAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCCGCCGGGCTTTTCGTTGTGTGGCAAAAAAAccccgc                                                                                                                                        >  1:86297/1‑248 (MQ=255)
                                                                                                               caggcaaTAACGGAGAACGGTTTGAAAATGATGGTGCACCCCCCGGCTACGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGGCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCAACGGCGGTACTGGCGGCTGGGCTTGACGTTGTGCGACAAAATAACCCGAAc                                                                                                                                       <  2:73348/245‑1 (MQ=255)
                                                                                                               cAAGTCATAACCGATAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGGGTTTCGTTGTGGGGCaaaaataccccac                                                                                                                                        >  1:553276/1‑241 (MQ=255)
                                                                                                                    aatacccGAGAACGGGTTGAAGATTATGCTTCACCCCCCGGCTTCGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCgg                                                                                                                                   <  2:195188/245‑1 (MQ=255)
                                                                                                                         ccGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCGTGGGGTTTCCGTTTGCGCCaaaaaaaacccaaaccggccgcgg                                                                                                                             >  1:296902/1‑231 (MQ=255)
                                                                                                                                 ggTTTGCAGATTATGCTGCCCCTCCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtgatga                                                                                                                     <  1:299304/251‑1 (MQ=255)
                                                                                                                                  ggttgcagagtaagctcgcccgcccGGCTTGGCCCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCGGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtgatg                                                                                                                      <  2:76424/227‑1 (MQ=255)
                                                                                                                                     tGCAGATTATGCGGCACCCCCGGGTTCAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCt                                                                                                                  <  2:109706/250‑1 (MQ=255)
                                                                                                                                               gctgcACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAATTGCGGGTGGGGCTTTCGGTTTGCGGGAAAATAAACCGAACGCGGCGGCGGGGTGAGGGGCTTCCTTGTCTcc                                                                                                       >  1:502903/1‑251 (MQ=255)
                                                                                                                                                                  gCTCCCGCCACAGCGGTGTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCTTCTGATACGCctgc                                                                                     <  2:236698/251‑1 (MQ=255)
                                                                                                                                                                                    tCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACt                                                                    <  1:254085/250‑1 (MQ=255)
                                                                                                                                                                                      aGGCTGGGCGCTTCTGCATCCGGTGTTCTTTTGCGTAACGGAACAATGGAACTGTCGCTATGCTCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGGTCCGGACGCGTAATCGGCGGTACTGGCTGCTGGTCTTTACTTTGTGCGACAAAATAACCCGAACGCGGCTTCGGAGTGATGAGCCTCAGCGTCTCCAGCGTCTGATTCGCCTCCTGCACGCTGCTGATACTGAc                                                                 <  2:190111/251‑1 (MQ=255)
                                                                                                                                                                                                         ccGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTaa                                               <  1:131470/250‑1 (MQ=255)
                                                                                                                                                                                                          cGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGTGCGCTGCGCCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGc                                             <  1:87428/251‑1 (MQ=255)
                                                                                                                                                                                                                    ttGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCg                                    <  1:77339/250‑1 (MQ=255)
                                                                                                                                                                                                                          aaCGGAACAATGGAACTGACGCTATGCGCCTGCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCAGCGGCTGTACTGGCGGCGGGGCTTTACGTTGTGCGACAAAATAACCCGAACTCGGCTGCGGAGTGATGAGCGTCATCGTCTCCAGCGTTTGATACGCCTGCTGCACGTTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCaa                             <  2:574376/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAAcgcg                                                                                                                                    <  1:310824/141‑1 (MQ=255)
                                                                                                                                                                                                                                 cAATGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGAc                                                   <  1:521995/222‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTGGAACTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTccc                     <  2:214487/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aCTGTCGCTATGCGCCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGAc               <  1:520616/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggcgcCTTCAGCATTTCCAGCACCTGATCCCACTGGGTAATTTCCGCCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTTATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGCCGGCAATTTCTCCCCTTGACGATACAGccc     <  2:376825/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    cgcCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTTACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCtt   <  1:241496/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     gcCTTCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTg  <  1:337453/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
ATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGACCGGCTAACAAGTCATAACCGAGAACGGTTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTCAGGCTGGGCGCTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTGGCGGCTGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCTTCATCGTCTCCAGCGTCTGATACGCCTGCTGCACGGTGCTGATACTGACGCCGTGCTCCTGACTTAAGCTGCGCACCGACGGCAATTTCTCCCCGTGACGATACAGCCCTTG  >  NC_000913/4571033‑4571526

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