Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,099 C→T intergenic (‑130/‑223) yifB ← / → ilvL putative magnesium chelatase YifB/ilvXGMEDA operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,0990CT97.8% 127.2 / NA 45intergenic (‑130/‑223)yifB/ilvLputative magnesium chelatase YifB/ilvXGMEDA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/27);  minor base G (0/1);  total (17/28)
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.96e-01

GTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTG  >  NC_000913/3949858‑3950342
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
gTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCaaaa                                                                                                                                                                                                                                               >  1:550921/1‑250 (MQ=255)
       caccaTCGTTAAGCCTGGTAGCCCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACACCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCAGTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCg                                                                                                                                                                                                                                       <  2:119694/251‑1 (MQ=255)
       caccaTCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCg                                                                                                                                                                                                                                       <  1:285104/251‑1 (MQ=255)
        accaTCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCgg                                                                                                                                                                                                                                      <  1:377520/251‑1 (MQ=255)
              gTTAAGCCGGGTAGACCTTTACTGATATGTTCCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAg                                                                                                                                                                                                                                <  2:213739/251‑1 (MQ=255)
                                  acagATATGTACCTCTACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCTCGGTTATGAACAATTGACATTGACATAAGCCCTTCTTGAGTCACCATTATGTGCAAAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTTGCTGATTCCTGTTTATTTTTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTACGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGa                                                                                                                                                                                                           <  2:329047/248‑1 (MQ=255)
                                   cTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGAATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGAAGCGGACAAAGGATGAGCTACgcgagaggga                                                                                                                                                                                                            >  1:559821/1‑241 (MQ=255)
                                     gATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGaaca                                                                                                                                                                                                        <  1:332372/251‑1 (MQ=255)
                                      atatGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGaacaa                                                                                                                                                                                                       <  1:94569/251‑1 (MQ=255)
                                            aCCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTc                                                                                                                                                                                                  >  1:307624/1‑250 (MQ=255)
                                            aCCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCa                                                                                                                                                                                                 <  1:276152/251‑1 (MQ=255)
                                                cAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATAc                                                                                                                                                                                              <  1:397708/250‑1 (MQ=255)
                                                     gggATCGGGGGCGCATTGACGCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTTGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAAt                                                                                                                                                                                        <  2:280229/249‑1 (MQ=255)
                                                           gggggCGCAGTTACTCCCAGTGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGATTCACCATTATGTGCATAAGATATCTCTGCTGTAGCCCGCTAATTCGTTAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAATTGAAGTTGAGTTGTTCTGGCGGTTGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACCACGAGGAAGGGAACAACATTCATACtgaatttgaatt                                                                                                                                                                                  <  2:335525/251‑1 (MQ=255)
                                                               gcgcATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                                                                                                                              <  1:398786/251‑1 (MQ=255)
                                                                  cATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTtca                                                                                                                                                                           <  1:473711/251‑1 (MQ=255)
                                                                   aTTTACTCCCAGGGCTGCGCGGGTATGGACAATTGACAGTGACATAAGCCCTCCTTGAGTAACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTACTGTTTATTTGTGCACGTGAAGTTGAGTTGTTCTGGCGGTGGAATGTTGTTCGCAAAAATGCAGCGGACAAAGGATGAAATACGAAGGAGGGGACCAAATTCATACTGAACATGAATTTTTTtcac                                                                                                                                                                          >  2:285611/1‑251 (MQ=255)
                                                                           ccAGGGCTGAGCGGGTATGAACTATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATAGCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTAATTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTGTCACTCACtata                                                                                                                                                                  <  2:6361/251‑2 (MQ=255)
                                                                                           aTGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACGACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAaaca                                                                                                                                                  <  1:515250/251‑1 (MQ=255)
                                                                                                     gACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTTGAATGATGGTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACTACATTCATACTGAAAATGAATTTTGTTCACTCACTATTTTATTTTTAAAAAACAACAATTtaa                                                                                                                                        <  2:239472/250‑2 (MQ=255)
                                                                                                           gACATTAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTTGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                                                                                                                                  <  2:199445/250‑1 (MQ=255)
                                                                                                               tAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAgg                                                                                                                                                                                                             <  2:263347/172‑1 (MQ=255)
                                                                                                               tAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAgg                                                                                                                                                                                                             >  1:263347/1‑172 (MQ=255)
                                                                                                               tAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGATGATTCCTGTTTATTTGTGCAAGTGAAATTGAGTTGTTCTGGCGGTGGAATGATGTTCGAAAAAATGCAGCGGACAAAGGATGAACTCCGGGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTTAAAAACAACCATTTATATTGAAattat                                                                                                                             >  1:246794/1‑251 (MQ=255)
                                                                                                                   cccgcctTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAAGGATGTTCGCAAAAATGCAGCGGACAAAGGATTAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTGTTTTCACTCACTATTTTATTTGTAAAAAACAACAATGTATATTGAAATGATTaaa                                                                                                                         <  2:524491/247‑1 (MQ=255)
                                                                                                                   cccTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGCATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGGTGAACTACGGGGAAGGGAACAAAATTCATACTGAAATTGAATTTTTTTAACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaagtttcga                                                                                                                          >  2:86862/1‑243 (MQ=255)
                                                                                                                   cccTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                                                                                                                                  <  1:86862/242‑1 (MQ=255)
                                                                                                                    catcctTGAGTCACCATTATGTGCATAAGATATCGCTGATGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGAGAAGTTGAGTTGTTCTAGCGGTGGAATGGTGTTCGCAAAAATGCAGCGTGCAAAGGCTGAAATACGAGGTAGGGAACAACATGCCTAATGAAAATGCAATTTTTTCACTCACAATTTTTTTTTTAAAAAAAACAAATTTATATTTAAATTATTaaa                                                                                                                         >  2:346588/3‑250 (MQ=255)
                                                                                                                           gagtcgtcagtatttgCATAAGATATCGTTGCTGTTGTCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATGTGTGCAAGTGAAGTTGAGTTGTTCTGTCTGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGATGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCCTTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                                                                                                                  <  2:293394/236‑1 (MQ=255)
                                                                                                                            aGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGAGGACAAAGGATGAACTACGAGGAAGGGAACAACATTTATATTGAAATTGAATTTTTTTCACTCAATATTTTATTTTTAAAAAACAAAAATTTATATTGAAATTTTTAAACGcataat                                                                                                                 >  2:374207/1‑247 (MQ=255)
                                                                                                                               cACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                                                                                                             <  1:556010/251‑1 (MQ=255)
                                                                                                                                  cATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                                                                                                                                     <  1:195176/224‑1 (MQ=255)
                                                                                                                                  cATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                                                                                                                                     >  2:195176/1‑224 (MQ=255)
                                                                                                                                                            ggaGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTa                                                                                <  2:550921/249‑1 (MQ=255)
                                                                                                                                                             tAGCCAGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGCATGATGTTCGCAAAAATGCAGCGGAACAATGATGAACTACGGGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACACACTATTTTATTTTTAAAAAACAACACTTTATATTGAAATTTTTAAACGCATCATAAAAATTGGCCAAAAAATATCTTGTAATATTTa                                                                                >  2:406613/1‑250 (MQ=255)
                                                                                                                                                              aGCCCGCTAATTCTTTAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTAc                                                                               <  2:106998/250‑1 (MQ=255)
                                                                                                                                                                            tGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAAc                                                                >  1:28702/1‑251 (MQ=255)
                                                                                                                                                                                       ggCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTAAGGGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGTAATTATTAAACGCATCATTAAAATCGGCCAAAAAATTTCTTGTACTTTTTACCAAACCTATGGTAACTCTTTTGGCAt                                                     >  2:276552/1‑251 (MQ=255)
                                                                                                                                                                                                          gtgCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCaa                                  <  1:409132/251‑1 (MQ=255)
                                                                                                                                                                                                              ccgtgcggttgcggggtgcgttctttGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGaa                               >  1:12191/26‑250 (MQ=255)
                                                                                                                                                                                                                       tgagttgTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACaaaa                      <  1:203503/250‑1 (MQ=255)
                                                                                                                                                                                                                       tgagttgTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAAt                     >  1:460351/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAAGAATGTAATAACATTCATCCTGAAATTGAAATTTTTTCTCTCACTATTTTATTTTTAAAAAACAACAATATATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAACATGACAGCCCTTCTAc       <  2:465556/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ggAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAACTGACAGCCCTTCTACg      >  2:269106/1‑249 (MQ=255)
                                                                                                                                                                                                                                           atgatgTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGtgtg  >  1:151333/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
GTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTG  >  NC_000913/3949858‑3950342

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