Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,321,672 C→T T425T (ACG→ACA trpE ← anthranilate synthase subunit TrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,321,6720CT98.2% 176.7 / ‑4.2 55T425T (ACG→ACAtrpEanthranilate synthase subunit TrpE
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/27);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGC  >  NC_000913/1321433‑1321910
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
gcgcGTGCTTTGTGACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATAtt                                                                                                                                                                                                                                     <  2:409699/251‑1 (MQ=255)
 cgcgGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTGCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAGCTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTc                                                                                                                                                                                                                                    <  1:358942/251‑1 (MQ=255)
   cgGGCTTTGTTACGGGTGTCGGCGGTTTCCGAGTGCGGAAAAGAATCAAGGACTACACCAGCAGCCGCTGGCACGGTGGCGATACCGTTTTCCACCAGCGGCGAGCGTATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGTCGCGGTAGGTGCCGCGGGGACTACCCTCCGCCTCGGCAATTAACTGCAGAGTGCGGATTCTCGGCTCGCCGCTTAATTTCCCCAAATTCAt                                                                                                                                                                                                                                  <  2:85782/251‑1 (MQ=255)
     ggCTTTGTTACGGTTTTCGTCGGCATCCGACTGCTGAACAGAATCAAGGACTACACCAGCCCCCTCTTGCACGGTGGCGATACCGTTTTCCACCAGCTCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCTGTGAAAGAACCGACCGCTCCGCCGTATCTGCCGCGGCGACGTGCTGCCGCCTCTGCAAGTAACTGCGCAGCGCGTACTGTCGGCGCCCCGCTTAATGTCCCCATAGTCAttc                                                                                                                                                                                                                                <  2:267593/251‑3 (MQ=255)
        tttGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGGCTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAgg                                                                                                                                                                                                                             >  1:240336/1‑251 (MQ=255)
                          ggCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAggg                                                                                                                                                                                                           >  1:182208/1‑251 (MQ=255)
                                             gAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGc                                                                                                                                                                                          >  1:163393/1‑249 (MQ=255)
                                              aaTCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGCCGCAGt                                                                                                                                                                                       >  2:564126/1‑251 (MQ=255)
                                               aTCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACAAGCGCCGAGCGGATCACAATGCCGGTGTCGAGATCGCCACGCGCGGTGAAATAACCTACCGCGCCGCCGTCGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTACTGTCCCCATATTCATACAGGAGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGtt                                                                                                                                                                                      >  1:373836/1‑251 (MQ=255)
                                                   aGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGcc                                                                                                                                                                                  >  1:474751/1‑251 (MQ=255)
                                                           cccAGCCCCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTAcgc                                                                                                                                                                          <  1:252764/250‑1 (MQ=255)
                                                                cACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCgagag                                                                                                                                                                     <  1:145216/251‑1 (MQ=255)
                                                                cACCCGCTTGCACGGGGGCGATACCTTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCAGAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCgagag                                                                                                                                                                     <  2:433203/251‑1 (MQ=255)
                                                                 aCCCGCTTGCACGGTGACGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCgag                                                                                                                                                                       >  2:67556/1‑248 (MQ=255)
                                                                          cACGGTGACGCTACCGTTTTCCACCAGCGCCGAGCGGATCACCATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCCCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCACACCGCTTCATGTCCCCATATTCCTCCAGGCGAGATCAGCGACCAGTGCGTCAAGATCGTGACGCCGTTCGCCGCCTACGCGAGAGACGAGGTGc                                                                                                                                                            >  2:256889/1‑250 (MQ=255)
                                                                             ggTGGCGATACCGTTTTCCACCAGCGCCGAGCGGAACACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTGCCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATAGTCATACAGGCGCGATAAGCGTGCAGGGCGGCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATca                                                                                                                                                        <  2:240336/251‑1 (MQ=255)
                                                                               tGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAATTAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGCCCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGCTCGTGACGCAGTTCGCCGACTACGCGAGAGACGCGGTGCATccc                                                                                                                                                       >  2:469669/1‑248 (MQ=255)
                                                                                  cGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGTCGCACCGCTTAATGTCCCCATATTCATACATGCGCGATAAGCGTTCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGAGCATCACAt                                                                                                                                                     >  1:117126/1‑249 (MQ=255)
                                                                                      aCCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAAt                                                                                                                                               <  1:8314/251‑1 (MQ=255)
                                                                                                  caGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACtt                                                                                                                                   <  1:209635/251‑1 (MQ=255)
                                                                                                   aGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGCAGCTGCCGCGGCGACGAACTTCAGCCCCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACttt                                                                                                                                  >  2:575511/1‑251 (MQ=255)
                                                                                                       ccGAGCGGATCACAATGCAGGGGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGGGCAGCACATAGGAATAACGGTCAACTTTGGTg                                                                                                                              <  2:108009/251‑1 (MQ=255)
                                                                                                         gAGCGGATCACAAGGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCTTGACGCAGTTCGCCGACTACGCGAGAGACGAGGGGCAGCACATAGGAATAACGGTCAACTTGGGTgag                                                                                                                            <  1:469669/251‑1 (MQ=255)
                                                                                                          aGCGGATCACAAGGCAGGTGTCGAGAGCGTCATGCGCGGAGAAATAACCTACCGAGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGTGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTCAGCGAGATAAGAGGTGCATCACATAGGAATAACGTGCAACTTTGGTgaga                                                                                                                           <  2:324165/251‑1 (MQ=255)
                                                                                                               aTCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGc                                                                                                                      >  2:444926/1‑251 (MQ=255)
                                                                                                                         aGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCgg                                                                                                            >  2:411718/1‑251 (MQ=255)
                                                                                                                               cGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCGCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTCCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCg                                                                                                      <  1:520086/251‑1 (MQ=255)
                                                                                                                               cGAGATCGCCATGCGCGGGGAAATAACCTACCGCGCCGCCGTAGCTGGCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCAGAGCGCGTACTTTCGGCGCACCGCTTAAGGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGTTGCATCACATAGGAAGAACGGTCAACTTTGGTGGGATCGGCGACGTAGCGGCTGcc                                                                                                       <  2:260341/250‑1 (MQ=255)
                                                                                                                                 agaTCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCggg                                                                                                    <  2:493290/251‑1 (MQ=255)
                                                                                                                                     cGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGt                                                                                                 <  1:96064/250‑1 (MQ=255)
                                                                                                                                      gCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGc                                                                                               >  1:328561/1‑251 (MQ=255)
                                                                                                                                         aTGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCaa                                                                                             >  2:455371/1‑250 (MQ=255)
                                                                                                                                         aTGCGCGGTGAAATAACCTACAGCGCCGCCGTAGCAGCCGCGGCGACGACATTTCGCCACGGCAATTAACTGCATAGCGCGTACTTTCGTCGCACCGCTTAATGTCCCCATATTCCTTCAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGGGAGACGAGGAGCATCACATAGGAATAACGGTCCACTTTGGTGAGATCGGCGCCGTAGCGGCTGCCggggggtcaac                                                                                            >  2:75353/1‑245 (MQ=255)
                                                                                                                                           gcgcGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATg                                                                                          <  1:564126/251‑1 (MQ=255)
                                                                                                                                                    actAACCTACCGCGCAGCCGTAGCTGCCGCGGTGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTACTGGCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTAAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGCTCGGCGACGTAGCGGCTGCCGGGGGTgaacaagggtgccagg                                                                                 >  2:501515/3‑236 (MQ=255)
                                                                                                                                                    aaTAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGcc                                                                                    >  1:180259/1‑248 (MQ=255)
                                                                                                                                                    aaTAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCACGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCACCATATTCATACAGGCGCGATAAGCGAGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATAGCCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGcc                                                                                    >  1:364871/1‑248 (MQ=255)
                                                                                                                                                              cgcgCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACggg                                                                       >  1:111361/1‑251 (MQ=255)
                                                                                                                                                               gcgcCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGAt                                                                                >  2:167850/1‑241 (MQ=255)
                                                                                                                                                               gcgcCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGGAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGAt                                                                                <  1:167850/241‑1 (MQ=255)
                                                                                                                                                                  ccgccgTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATAGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCgag                                                                   <  2:433622/251‑1 (MQ=255)
                                                                                                                                                                      cgTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCa                                                               <  1:272408/251‑1 (MQ=255)
                                                                                                                                                                       gTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGt                                                                                                                                                              <  1:581907/155‑1 (MQ=255)
                                                                                                                                                                       gTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCAACGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGt                                                                                                                                                              >  2:581907/1‑155 (MQ=255)
                                                                                                                                                                             gCCGCGGCGTCGTCCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATGCAGACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATTCGTGCCAGATCATTACGTGCGAGCCCAACCAg                                                          <  2:505450/249‑1 (MQ=255)
                                                                                                                                                                               cgcgGCGACGACCTTCCGCCTCTGCCATTACCTGCAGAGCGCGTACTTTCGGCGCTCCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCATGGCGTCAGGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGTTGCATCACATAGGAAGAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCACATGCGTGCCAGATCATTTCGGGCGATACCATCCAGCATc                                                      <  2:190183/251‑1 (MQ=255)
                                                                                                                                                                                 cggcgACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATc                                                      >  1:128543/1‑249 (MQ=255)
                                                                                                                                                                                 cggcgACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAg                                                    <  2:588197/251‑1 (MQ=255)
                                                                                                                                                                                                 ccTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGctct                                    <  1:423356/251‑1 (MQ=255)
                                                                                                                                                                                                                         gTACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGt            >  2:460561/1‑251 (MQ=255)
                                                                                                                                                                                                                          tACTTTCGGCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGtt           >  2:131230/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ggCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACg    <  1:411718/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ggCGCACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACg    <  1:482643/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgcACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCGTTATGATCGGTACGCATTTCCAGTTCAATACgg   <  1:455371/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cgcACCGCTTAATGTCCCCATATTCATACAGGCGCGATAAGCGCGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCCACTTTGGTGAGATCGGCGCCGTAGCGGCTGCCGGGGGTGCAACTGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGc  >  2:345021/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
GCGCGGGCTTTGTTACGGGTTTCGTCGGCTTCCGACTGCGGAACAGAATCAAGGACTACACCAGCACCCGCTTGCACGGTGGCGATACCGTTTTCCACCAGCGCCGAGCGGATCACAATGCAGGTGTCGAGATCGCCATGCGCGGTGAAATAACCTACCGCGCCGCCGTAGCTGCCGCGGCGACGACCTTCCGCCTCGGCAATTAACTGCATAGCGCGTACTTTCGGCGCACCGCTTAACGTCCCCATATTCATACAGGCGCGATAAGCGTGCAGGGCGTCAAGATCGTGACGCAGTTCGCCGACTACGCGAGAGACGAGGTGCATCACATAGGAATAACGGTCAACTTTGGTGAGATCGGCGACGTAGCGGCTGCCGGGGGTGCAAATGCGTGCCAGATCATTACGGGCGAGATCAACCAGCATCAGATGTTCAGACAGCTCTTTATGATCGGTACGCATTTCCAGTTCAATACGGC  >  NC_000913/1321433‑1321910

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