Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,401,893 G→A Q455* (CAA→TAA)  abgB ← p‑aminobenzoyl‑glutamate hydrolase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,401,8930GA97.8% 148.7 / NA 47Q455* (CAA→TAA) abgBp‑aminobenzoyl‑glutamate hydrolase subunit B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/19);  minor base C (1/0);  total (27/20)
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

GAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCAT  >  NC_000913/1401654‑1402129
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
gAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAGCCTTGTTACCAATTCTGTCGACCCAGCCATAGATCTTCCCGGATTGGGACGACGACGGTATGGATTACATACTCATTGGCATTCCACGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTTTCACGTTTTTCTGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGc                                                                                                                                                                                                                                     <  1:90162/250‑1 (MQ=255)
       aTAAGCAAAGAGATAAAGAGCAGAAAAGGATGGGGTACCTTGTTTCGAATTCGTTCGAACCAGCCATAGAGCTTCCCGGATTGGGAGTACGACGGTATGGGTGACATACTCATGGGCTTTTCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTGTAAAGGTGACGGTGTCACGTTTTTCTGGTTAGTGCAGTGATACGGTTGCGTGTCCGTTACGTGCTGATGCTCTTGTTAGCATTCTTGCAATAGccc                                                                                                                                                                                                                             <  2:223927/251‑1 (MQ=255)
          agcAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGcc                                                                                                                                                                                                                              <  2:183202/245‑1 (MQ=255)
          aTCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGcc                                                                                                                                                                                                                              >  1:183202/1‑247 (MQ=255)
               aTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCaa                                                                                                                                                                                                                     <  1:123306/251‑1 (MQ=255)
                   aTAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGa                                                                                                                                                                                                                  >  1:430415/1‑250 (MQ=255)
                    tAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCGGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGGATCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGACTCAAGGa                                                                                                                                                                                                                  >  2:433261/1‑249 (MQ=255)
                       agagCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAag                                                                                                                                                                                                                >  1:312984/1‑248 (MQ=255)
                         agCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGAGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGAtt                                                                                                                                                                                                          >  1:362914/1‑251 (MQ=255)
                               aaaGGATGGGGAACCTTTTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACtgt                                                                                                                                                                                                    <  2:312984/251‑1 (MQ=255)
                                  ggATTGGGAACCTTGTTTCCAATTCGTTCTACCCAGCCATAGAGCTTCCCGGGTTGGGAGGTCGACGTTATGGATGACATACTCATGGGCATTCCTCGGTTTTTGTGTTTGTGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGTGATAGGGCAGTGATACGGGTGCGTGTCCGTTACTTGCGGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGt                                                                                                                                                                                                 <  2:312434/251‑1 (MQ=255)
                                   gatggaaaacattGTTAACAATTCTTTCGACCCAGCCATAGAGCTTCCGGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGAAGGGGTCACGTTTTTCGGGATAGGGCTGTGATCCGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCACTAGCCCTGATTCAAGGACGAGATTGGCTGTGGtt                                                                                                                                                                                                >  2:529200/12‑251 (MQ=255)
                                           aCCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGCGATTGACTGTGGTTGCTGCCAt                                                                                                                                                                                        >  2:295061/1‑251 (MQ=255)
                                                 ttACCAATTCTTTCGACCAAGCCATAGAGCTTCCCGGATTGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGACTCAAGGAAGAGATTGACTGTGGTTGCTGCCCTAGTTTTc                                                                                                                                                                                 >  2:260370/1‑251 (MQ=255)
                                                   aCCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTc                                                                                                                                                                                 <  1:79721/250‑1 (MQ=255)
                                                     cAATTCTTTCGACCCAGCCATAGACCTTACCGGATTGGGCGGACGAAGGTCTGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGATTGCGTTGATGTTGTAATTTTAAAGGAGACGGTGTCACGATTTCCGGGATAGGACAGTGATAAGGTTGAGTGTCCGTTACTTGCTGCTGCTCTTGTTCGCATTCTTGCAATAGCCCTGACTCAAGGAAGGAGTTGTCTGTGTTTGCTGCCATAGTTTTcac                                                                                                                                                                               >  2:534714/1‑248 (MQ=255)
                                                      aaTTCTTGCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTcgccg                                                                                                                                                                             <  1:223260/251‑1 (MQ=255)
                                                                 cccAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCAtt                                                                                                                                                                  <  1:413352/251‑1 (MQ=255)
                                                                      ccATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCttt                                                                                                                                                             <  2:193049/251‑1 (MQ=255)
                                                                          agagCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATg                                                                                                                                                          <  1:367686/250‑1 (MQ=255)
                                                                                        tGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTc                                                                                                                                            <  1:295061/250‑1 (MQ=255)
                                                                                         gggAGGACGACGGTCTGGATGACATACTCATGGGGATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTTATGCTCTTGTTAGCATTCTTGCAATAGCCCTGCATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCg                                                                                                                                           <  2:383794/250‑1 (MQ=255)
                                                                                                               cATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACg                                                                                                                    >  1:145366/1‑251 (MQ=255)
                                                                                                                 tACTCATGGGCATTCCTCGGTTGTTGTGTTTGAGTTGTTGTTGTTATTTTCAAGGTGACGGTGTCACGTTTTTCGGGATAGAGCAGTGATACGGTTGATTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGTGAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCATCAGCATTCCTTTATGAGCAATATATGTTCGCCCCTGGCTAACCCGTTGCCACGtt                                                                                                                  >  2:135376/1‑250 (MQ=255)
                                                                                                                                 gcGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTGCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTAcc                                                                                                    <  2:430415/248‑1 (MQ=255)
                                                                                                                                                gCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAgggg                                                                                    <  1:401815/250‑1 (MQ=255)
                                                                                                                                                               ttttAAAGGCGACGGAGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGGATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCACCTGGCTAACCAGTTGCCACGTATGACGCGGTGTACCGACGGCAAAACAGGGGTTGAACCACTGGGCa                                                                     >  2:553013/1‑250 (MQ=255)
                                                                                                                                                                        tGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAgcagca                                                                                                                                                                    >  2:199377/1‑146 (MQ=255)
                                                                                                                                                                        tGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAgcagca                                                                                                                                                                    <  1:199377/146‑1 (MQ=255)
                                                                                                                                                                         gACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAACACAGGGGCTGAAACACTGGGCACAAGGCAGttt                                                          >  2:557379/1‑251 (MQ=255)
                                                                                                                                                                            ggTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTcc                                                        >  1:294271/1‑250 (MQ=255)
                                                                                                                                                                                   cGTTTTTCGGGATAGGGCATTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAAGCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACg                                                <  1:367223/251‑1 (MQ=255)
                                                                                                                                                                                      ttttCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCg                                             >  2:157477/1‑251 (MQ=255)
                                                                                                                                                                                         gcGGGATTGGGCAGTAATCGGTTTGCGGTTCCGTTACTTGCTGATGCTCATGTTTGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGGTGCTGCCATAGTTTGCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCTACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGtcgtcgtcg                                          <  2:440267/250‑4 (MQ=255)
                                                                                                                                                                                               tAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGATGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATc                                     >  1:420287/1‑250 (MQ=255)
                                                                                                                                                                                                  ggcagttagtcGGTGGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCGGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAGCAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTc                                 <  2:188278/241‑1 (MQ=255)
                                                                                                                                                                                                    cAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCg                                >  1:429320/1‑250 (MQ=255)
                                                                                                                                                                                                     aGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGAt                              >  1:19832/1‑251 (MQ=255)
                                                                                                                                                                                                      gTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTTTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACCTCAGTCGATg                             >  2:92558/1‑251 (MQ=255)
                                                                                                                                                                                                           aCGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGc                         >  1:108731/1‑250 (MQ=255)
                                                                                                                                                                                                                     tgtCCGTTACATGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGACTCAAGGAAGAGATTGACTGTGGTTGCTGCCCAAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCATTAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTCGCGGTGTACCGGCGGCAAAACCAGGGCTGAAACACTGGGCACCATGAGGTTTCCAACTGACGCCGCCAACATCAATCGATGCCGCAAGCACGta                >  2:188709/1‑249 (MQ=255)
                                                                                                                                                                                                                           ttACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTg         >  1:236762/1‑251 (MQ=255)
                                                                                                                                                                                                                            tACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTgg        >  1:457007/1‑251 (MQ=255)
                                                                                                                                                                                                                            tACTTGCTGATGCTCTTGTTAGCATTATTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCAGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTgg        >  1:444570/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCt                                                                                                                                                               >  1:376067/1‑95 (MQ=255)
                                                                                                                                                                                                                                 gCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCGGCCATAGTGTTCGCCGCCAGCAGCATTCCt                                                                                                                                                               <  2:376067/95‑1 (MQ=255)
                                                                                                                                                                                                                                   tGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCat  >  1:230520/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
GAGTACGATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCAT  >  NC_000913/1401654‑1402129

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