Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,174,041 A→T E320D (GAA→GAT birA → DNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,174,0410AT98.3% 210.9 / NA 59E320D (GAA→GATbirADNA‑binding transcriptional repressor/biotin‑[acetyl‑CoA‑carboxylase] ligase BirA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (28/30);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCG  >  NC_000913/4173802‑4174283
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
aTGCTAATACGTGAATTACGGGCTGAGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGaa                                                                                                                                                                                                                                           >  2:136746/1‑249 (MQ=255)
 tGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAgg                                                                                                                                                                                                                                        >  1:419405/1‑251 (MQ=255)
        aCGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGTGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTAt                                                                                                                                                                                                                                 <  2:143356/251‑1 (MQ=255)
            gAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAGCCCTGGATGTGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGc                                                                                                                                                                                                                             <  1:557559/251‑1 (MQ=255)
            gAATGACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGc                                                                                                                                                                                                                             <  2:386935/251‑1 (MQ=255)
             aaTTACGTGCTGCGCTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGCTAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCATGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCt                                                                                                                                                                                                                            >  1:365337/1‑251 (MQ=255)
                tACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTcc                                                                                                                                                                                                                          <  2:496565/250‑1 (MQ=255)
                 aCGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTTCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTcccc                                                                                                                                                                                                                        <  1:273842/251‑1 (MQ=255)
                                   ctTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGCAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGt                                                                                                                                                                                                        >  1:31565/1‑249 (MQ=255)
                                             gaagGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACctct                                                                                                                                                                                            >  1:587832/1‑251 (MQ=255)
                                              aagGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTGTATTAATCGCTCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGTCGTTTCACGCGGAATAGACAAACAGTGGGCTTTATTACTTGAGCATGATGGAATAATAAAACTCTGGATGTGCTGTGAAATAGCCCTGCGTAGTGCAGATAAATAATAAAGGGGAGTTTTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCtctt                                                                                                                                                                                           <  1:236506/251‑1 (MQ=255)
                                                                 gTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCAc                                                                                                                                                                         <  2:16871/250‑1 (MQ=255)
                                                                     cgcTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTc                                                                                                                                                                    <  1:259418/251‑1 (MQ=255)
                                                                      gcTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTATCACTTTTCg                                                                                                                                                                   <  1:446145/251‑1 (MQ=255)
                                                                         gggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACGTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTc                                                                                                                                                                 <  2:192103/250‑1 (MQ=255)
                                                                          ggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTc                                                                                                                                                                 >  1:588625/1‑249 (MQ=255)
                                                                          ggAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCAGTGGTGATAAAGAAAGATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATGTGCGTAGTCTGACCTCTTCTACCGCATGATGAGCACTTGTCGTCAg                                                                                                                                                               <  2:223469/251‑1 (MQ=255)
                                                                                     aTAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTgg                                                                                                                                                    <  2:375571/251‑1 (MQ=255)
                                                                                       aaTTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGCAATATTTGGCATTTCACTCGGAATAGACATACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGCTTAGCACTTTTCGTCAGGATTAAACTGgcg                                                                                                                                                  >  1:104748/1‑251 (MQ=255)
                                                                                                aaTCGCCCAGTGAAACTTATCATTGGTGATAAGGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTTACCGCTGCTACCGCATGATTAGCACTTTTCGGCAGGATTAAACTGGCGCGCTCACg                                                                                                                                          <  2:585453/250‑1 (MQ=255)
                                                                                                    gCCCAGTGAACCTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgt                                                                                                                                       >  1:190590/1‑249 (MQ=255)
                                                                                                      ccAGTGAAACTTATCATTGGTGATAAAGCAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggt                                                                                                                                   >  2:357168/1‑251 (MQ=255)
                                                                                                        aGTGAAACTTATCATTGGTGCTAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggt                                                                                                                                   >  2:276864/1‑249 (MQ=255)
                                                                                                          tGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTAATTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGAACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGa                                                                                                                                >  2:590896/1‑250 (MQ=255)
                                                                                                           gAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGACTAGACAAACAGGGCGCTTTATTACTTGAGAAGGATGGAATAATAACACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAACTTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAat                                                                                                                              >  1:64932/1‑251 (MQ=255)
                                                                                                              aCTTATCATTGGTGATAAAGAAATATTTGGCATTTCGCGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGAGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACATCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTAACGAGTCGGTAGAATAtt                                                                                                                           >  2:321275/1‑251 (MQ=255)
                                                                                                                    cATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCt                                                                                                                      <  2:150121/250‑1 (MQ=255)
                                                                                                                        ggTGATATAGAAATATTTGGCATTTCACTCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGATTATTCGCTCCCCTGCAAATTATTTGCGTGGTCTGACCTCTTCTATCTCATGAGTAGCACTTTTCGTCAGGATGAAATTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTaa                                                                                                                  <  1:535897/250‑1 (MQ=255)
                                                                                                                             tAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGCGTAGGTAGAATATTTTGCTTTAAGTTc                                                                                                              >  2:498692/1‑249 (MQ=255)
                                                                                                                              aaaGAAATATTTGTCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAACAACACCCTGGATGGGCAGTGAAATATCCCTGCGTAGTGCAGATAAATAAGACCGGGGAGTATTCGCTCACCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCCTGATTAGCACTTTTCGTCAGGATTAATCTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGc                                                                                                           <  1:313930/251‑1 (MQ=255)
                                                                                                                                            cATTTCACGCGGAATAGACAACCAGCGGGGTTTATTTCTTGAGCAGGACAGAATAATAACCACATGGCTGGGCGGTGCAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTAGCTCCCCTGCAACTTATTTGAGTAGTCTGACCTCTTCTACCGCCTGATTAGCACTTTTCGTCAGGATTAAACTGGCACGATCCCGAGTAGGTAGAATATTTTGCTTTAAGTTCAGTCAGTTGATCTCttt                                                                                             >  2:359725/1‑251 (MQ=255)
                                                                                                                                                   cgcgGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATTGAATAACAAAACCCTGGATTGGCGGGGAAATATCACTGAGTAGCGCAGATAAATCAGAAAGGGGAGTCTTCGCTCCCCTGCAAAAGCTTTGCGTAGTCTGACCTCTTCTACCGGATGATTAGCACTTATCGTCATGCTTACACTTGCTCGCGCTCGAGTAGCTAGAATATTGTGCTTTAAGTCCACCCAGTAGATCTCATTCCACaa                                                                                       <  2:343517/250‑1 (MQ=255)
                                                                                                                                                             aCAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAg                                                                            <  1:412103/251‑1 (MQ=255)
                                                                                                                                                                 aCAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCtt                                                                        <  2:419405/251‑1 (MQ=255)
                                                                                                                                                                            tattaCTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGcttct                                                              <  1:389025/250‑1 (MQ=255)
                                                                                                                                                                                   tGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGGTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAg                                                       <  1:357168/250‑1 (MQ=255)
                                                                                                                                                                                       cAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAAt                                                  >  1:125104/1‑251 (MQ=255)
                                                                                                                                                                                           aTGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTcg                                              <  1:196025/251‑1 (MQ=255)
                                                                                                                                                                                           aTGGAATAATAAAACCCTGGATCGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTcg                                              >  1:222955/1‑251 (MQ=255)
                                                                                                                                                                                              gAATAATAAAGCCCTGGATGGGCGGTGAAATAACCCTGCGGAGTGCAGATAAATAAAAAAGAGTAGTAATCTCTCCCCTGCAAATTATTTGCTTAGTCAGACCTGTTCTACCGCATTATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACTAGTAGGTATAATATTTTGCTTGAAGGTCAGCCATTTGTTTCCTTGCCACAATGTCTTGGCAGTCTTAATCGTTTCGTCTTTAGTTAAGTTCGCGt                                           <  2:245585/251‑1 (MQ=255)
                                                                                                                                                                                                ataataCAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGt                                           >  1:421999/1‑249 (MQ=255)
                                                                                                                                                                                                ataataAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTTCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCGGTCTTAATCGCTGCTTCTTTAGTTAATTTCGCGTAg                                         <  1:225144/251‑1 (MQ=255)
                                                                                                                                                                                                ataataAAACCCGGGATAGGCCGTGAAATATCCCTGCGTAGTTCAGCTAAATAACAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTCCCGAATTATTAGCACTTTTCGTCAGGATTAATCCGGCTCGCTCACGAGTAGGTAGAATATTTTGCTGTAAGTTCAGCCAGTTGATCTCTTTCCACAATGACAGGGCAGTCTTAATCGCTTATTCTTTatcaaagttcgcgtag                                         <  2:53834/251‑16 (MQ=255)
                                                                                                                                                                                                 taataaAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGt                                        >  1:160083/1‑251 (MQ=255)
                                                                                                                                                                                                  aataaAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGtt                                       >  1:165050/1‑251 (MQ=255)
                                                                                                                                                                                                      aaaCCCTGGATGGGCGGTGACATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATg                                    >  2:376930/1‑250 (MQ=255)
                                                                                                                                                                                                      aaaCCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATg                                    >  1:411715/1‑250 (MQ=255)
                                                                                                                                                                                                          ccTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGTATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAAt                               <  2:412970/251‑1 (MQ=255)
                                                                                                                                                                                                               aTGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCCGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGa                           >  2:144776/1‑250 (MQ=255)
                                                                                                                                                                                                                tGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTc                                                                                                                                                                    >  2:165891/1‑112 (MQ=255)
                                                                                                                                                                                                                tGGGCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCAGGATTAGCACTTTTc                                                                                                                                                                    <  1:165891/112‑1 (MQ=255)
                                                                                                                                                                                                                 gggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAATAAATGGGAGTATTCGCTCCCCTGCAAATTATTTGCAGAGTCTGACCTCTTCTACCCCATGAATAACACGTTGCATCAGGAtt                                                                                                                                                           <  2:553561/120‑1 (MQ=255)
                                                                                                                                                                                                                 gggCGGTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGCAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGCTTAGCACTTTTCGTCAGGAtt                                                                                                                                                           >  1:553561/1‑120 (MQ=255)
                                                                                                                                                                                                                     ggTGAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCggg                    >  2:244111/1‑251 (MQ=255)
                                                                                                                                                                                                                        gAAATATCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCg                 <  2:411715/251‑1 (MQ=255)
                                                                                                                                                                                                                           atatCCCTGCGTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAGTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTGTAGTTAATTTCGCGTAGGTATGAAAATAGGAATCCGGGTCGGTa              <  2:285169/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gTAGTGCAGCTAAATAAGAACGGGGTGTATTAGCTCCCCCGCAACTTATTTGCGTACTCTGCCCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCCGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCtt    >  2:307753/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gTAGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCtt    >  2:155102/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aGTGCAGATAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTcg  >  1:493186/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  
ATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAATTTTATTAATCGCCCAGTGAAACTTATCATTGGTGATAAAGAAATATTTGGCATTTCACGCGGAATAGACAAACAGGGGGCTTTATTACTTGAGCAGGATGGAATAATAAAACCCTGGATGGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTTTCCACAATGTCATGGCAGTCTTAATCGCTTCTTCTTTAGTTAATTTCGCGTAGTTATGAAAATAGGAATCCGGGTCGGTAAAAGCCCCTTCG  >  NC_000913/4173802‑4174283

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