Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,468,768 C→T W362* (TGG→TGA chiA ← endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,468,7680CT96.9% 214.3 / NA 65W362* (TGG→TGAchiAendochitinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (38/24);  minor base G (0/2);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGT  >  NC_000913/3468524‑3469013
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
cccATGCGTCTTGATAGACAGCAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCg                                                                                                                                                                                                                                                     >  1:337564/1‑251 (MQ=255)
         cTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTa                                                                                                                                                                                                                                             >  1:310898/1‑250 (MQ=255)
            gATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGGACGCTCGAGACGCAATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGGTCTTGTAAAc                                                                                                                                                                                                                                          >  1:247194/1‑250 (MQ=255)
                 aCAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCGTGTAAACATTAcc                                                                                                                                                                                                                                    <  1:321868/251‑1 (MQ=255)
                  cAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTAAGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCACTTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGCGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTAcc                                                                                                                                                                                                                                    >  1:189144/1‑250 (MQ=255)
                   aGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCAt                                                                                                                                                                                                                                  <  1:315353/251‑1 (MQ=255)
                     aaGTTTTATAATGTGCAGTGTTATAATCTACGATGCGGTCCGCGCTGTACGTTTGGCGGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTg                                                                                                                                                                                                                                <  2:319004/251‑1 (MQ=255)
                      aGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCCCGCAGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACCTTACCATTg                                                                                                                                                                                                                                >  1:49046/1‑250 (MQ=255)
                             gagtttgactCGTTATATACTACGAAGCAGTCCGCGCTAGAGGCTGTGCTTGCCTGGAAATCGTCCGCAGCAAGGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGGTACCGTACTGACTCAAGTCCGCAGCGGTAGCGGTACGCTCGAGGCGCCAGGGGTTATTGGAGTTGTCAGCGTTTGTGCCCCGTGTGCGAGCTTTAGAGCCAACGCATCAGGCTTTCTGGGTAACAGTACCATCGAAAATAAc                                                                                                                                                                                                                        <  2:194682/242‑1 (MQ=255)
                             aaaggtgCACCGTTATAATCTACGATGCCGTACGCCTTTTACGCTTTGTTGTCCTGTAAACCATCCGCAACAAGGACGCCGCCGTTAGCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGCACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATGGAAAAGAAc                                                                                                                                                                                                                        <  2:322660/247‑1 (MQ=255)
                                 ttgCATCGTTATAATCTACGATGCGATCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATg                                                                                                                                                                                                                     <  2:104515/249‑1 (MQ=255)
                                     aTCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTaa                                                                                                                                                                                                                 >  1:341660/1‑250 (MQ=255)
                                       cgtaataATCTACGATGCTGTCCGGGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAtt                                                                                                                                                                                                             <  1:13281/247‑1 (MQ=255)
                                              aTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGCCGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGTCCCAACACATCCGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTAc                                                                                                                                                                                                      >  2:151504/1‑251 (MQ=255)
                                                  ccGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCt                                                                                                                                                                                                    <  2:341660/248‑1 (MQ=255)
                                                    gATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTCCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCt                                                                                                                                                                                               >  2:305777/1‑251 (MQ=255)
                                                     aTGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTg                                                                                                                                                                                              >  1:250240/1‑251 (MQ=255)
                                                           tCCGCGCTGTCCGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACTCCGCCGGTATCAATTTCACAGGTAGTCGGTTTACCGTGCTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGTCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGTCAATCTTTAGACCCAACCCAGCAGGCGTTCTTGTAAACATTCCCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGc                                                                                                                                                                                         <  1:47298/250‑1 (MQ=255)
                                                              gcgcTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTcc                                                                                                                                                                                      >  1:430362/1‑250 (MQ=255)
                                                               cgtggTACGCTTTGCTGGCCGGGAAACCATCCGCAACAAGGACGCCGCCGTTATCAATTTCACAGGGAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAATCATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCAc                                                                                                                                                                                    <  1:221569/247‑1 (MQ=255)
                                                                                    cgAAACCATCCGCCACAATGACGCCGCGGATATCACTTTCCCAGGAAGACGGGTTACCGAACTGACTCAATTCAGCCGCGATAGCTGTACGCCAGCGACGCCATGGGTTATTGGAGTTTTCTGCGTTGGTAACCCGTAGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTAAGACTAACATGGTAAATTTCGCTACCTTCCTGGCTGCTCCAGGCTTCCAATCAGTAATCGtt                                                                                                                                                                >  2:216682/2‑250 (MQ=255)
                                                                                     gAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGACAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTa                                                                                                                                                               >  2:135278/1‑250 (MQ=255)
                                                                                      aaaCCATCCGCAACAATGGCGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTTGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTTAAAATAACATGATAAATTTCGCTTCCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAAc                                                                                                                                                             <  1:33105/251‑1 (MQ=255)
                                                                                       aaCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAAGTACACGGGTTAACGTACTGACTCAAATCCGCAGCGGTAACTGTACGATCGAGACGCCATGGGATATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCCTGTAAACATTACCATTGAACATAACATGGTAAATTTCGCTACCTTCCTGGCAGCTCCAAGCTTGCAATGAGTAATCGCTAACg                                                                                                                                                            >  1:213788/1‑251 (MQ=255)
                                                                                        aCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGc                                                                                                                                                           >  2:81327/1‑251 (MQ=255)
                                                                                        aCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGc                                                                                                                                                           >  1:11298/1‑251 (MQ=255)
                                                                                               gCAACACTGACGACTCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCAATGCGTTATTGGAGTTTTAAGCGCTGGGACCCCGTGGGCAATCGTTAGACCCAACCCATAAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGGATATTTCGATCCCTTCCTGGACGCTCCACGCTTGCACTGCGTGATCGTTAAAGCTATCCGc                                                                                                                                                    >  2:184171/1‑251 (MQ=255)
                                                                                                caacaaTGCCGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTACTCGTTAACGCTATCCGCt                                                                                                                                                   >  2:365105/1‑251 (MQ=255)
                                                                                                    aaTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGgtg                                                                                                                                               >  1:15137/1‑251 (MQ=255)
                                                                                                     aTGACGCCGCCTTTATCAAGTTCACAGGTAGTCGGGTTACCGTACTGACTCAATGCCGCAGCGGTAGCTGGACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTGTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCGACCTTCCGGGCCGCTCCACGCTTGCAATGAGGAATCGGTAACGCTATCCGCTGgtgt                                                                                                                                              <  2:414291/251‑1 (MQ=255)
                                                                                                     aTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGgtgt                                                                                                                                              >  1:414291/1‑251 (MQ=255)
                                                                                                       gACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTGCCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTg                                                                                                                                            <  2:412354/251‑1 (MQ=255)
                                                                                                        ccgccgcCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGa                                                                                                                                           <  1:228170/250‑1 (MQ=255)
                                                                                                        aCGCCGCCGTTATCAATTTCACAGGTAGTCGGGATACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGAGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTgg                                                                                                                                           >  2:323073/1‑250 (MQ=255)
                                                                                                              ccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCCCTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTc                                                                                                                                     >  1:364224/1‑251 (MQ=255)
                                                                                                                    ccAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGtt                                                                                                                               <  2:337564/250‑1 (MQ=255)
                                                                                                                      aaTTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTAt                                                                                                                             >  2:379493/1‑251 (MQ=255)
                                                                                                                               ggTAGTAGGGCTACCCTCCTGCCTCCATTCCGCACCGGTAGCTGTACGCTCGAGCCGCCATGGCTTTTTGGTGTTTTCATCGCTGGTACCCCGTGGGCAATCTTTAGACCCCACCCATCTGACGTTCTTGTCAACATTACCATTGAAAATTACAt                                                                                                                                                                                                                      >  2:381851/1‑155 (MQ=255)
                                                                                                                               ggTAGGCGGGTTACCGTACTGACTCAATGCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTGATTGGAGTTTTTAGCACTGGTACCCCGGGGGCAATCTTAAGACCCAACCCAGCAGGCGTTATTGTAAACATTACCATTGAAAATAACAt                                                                                                                                                                                                                      <  1:381851/155‑1 (MQ=255)
                                                                                                                                       ggTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAcc                                                                                                           <  1:243101/251‑1 (MQ=255)
                                                                                                                                              gTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGttt                                                                                                    >  1:298862/1‑251 (MQ=255)
                                                                                                                                                   gACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCGATCCGCGGGTGTTGATGGCGCGCTGTTATCCGGAGTTGGGTTCACCGGTGTTTCCg                                                                                                 <  2:430362/249‑1 (MQ=255)
                                                                                                                                                       cAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGtt                                                                                            <  1:152924/250‑1 (MQ=255)
                                                                                                                                                        aaTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt                                                                                                                                                                                                              >  2:145786/1‑137 (MQ=255)
                                                                                                                                                        aaTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt                                                                                                                                                                                                              <  1:145786/137‑1 (MQ=255)
                                                                                                                                                             cGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCaa                                                                                     <  1:240798/251‑1 (MQ=255)
                                                                                                                                                                   ggTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCgg                                                                                 >  2:294611/1‑249 (MQ=255)
                                                                                                                                                                            aCGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGAACCCCGTGGGCAATCTTTAGACCCAACCCATCCGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATAAACCGGAGTCACAGCAc                                                                     >  1:408486/1‑251 (MQ=255)
                                                                                                                                                                                 cgagtggCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTCCCTTCCTGGCCGCTCCCCGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGAGGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCa                                                                <  2:364224/245‑1 (MQ=255)
                                                                                                                                                                                 cGAGACTCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCATCCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTGCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGcc                                                                 <  2:334257/250‑1 (MQ=255)
                                                                                                                                                                                    gaCGCCATGGGTTCTTGGAGTTTTCAGCGCTGGAACCCCGTGGGCAATCTTTAGACCCAAACCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGc                                                                                                                                                           >  2:413263/1‑159 (MQ=255)
                                                                                                                                                                                    gaCGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGc                                                                                                                                                           <  1:413263/159‑1 (MQ=255)
                                                                                                                                                                                               ttATTGGAGTATTAGGCGCTGGTAGCCCGTGGGTGATCTTTAGACCCAACCCCTCAGGCGCTCTTGTAACCATTACCCTTGACACTAACATGGTAACTTTCGCTACCTTCATGGCCGCTCCCCGCTCGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTGGATGACTGGCTGCTATCGGCTGTCGGGCTCACCGGTGTTTCCGGAGCTGGACCACCCGGGGTCAAAGCCCACCCATTATAATCCGGCtt                                                  >  2:410526/1‑251 (MQ=255)
                                                                                                                                                                                                       gTTTTCAGCGTGGGTACCCCGGGGGCATTTTTTAAACCCAACCCACCGGTTGTTCTTTAAAACATCACCATTGAAAAAAACATGATAATTTCCGCTACCTCCCGGGCCGCTCCCCGCTGGCAATGAGTAAGCTTTAACGCTACCCGCGGGTGTGTAGGGCGCGCTGTTATCCGGAGTGGGGGCCCCCGGTGTTTCCGGAGTTGAAGAAACGGTGCTCACAGCACCGCCATTATTATCCGGCTTAACGGaa                                           <  2:26105/250‑1 (MQ=255)
                                                                                                                                                                                                       gTTTTCAGCGCTGGTACCCCGAGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAACATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGTGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAAc                                          >  1:431381/1‑251 (MQ=255)
                                                                                                                                                                                                            cAGCGTTGGTCCCCCGTGGGCAATTTTTAGACCCACCCCATCAGGGTTTTTTGAAAACATTCCCATGAAAAAAAACATGATAATTTTCGCTCCCTCCCTGGCCGCCCCACGCTTGCAAGGAGTAACCGTTAACGTTACCCGCGGGTGTTGATGGCTCGCTGTTATCGGGAGTTGGGGTCACCGGTGTTCCGGGAGTTGGATCAACGGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGa                                      <  1:379493/250‑1 (MQ=255)
                                                                                                                                                                                                                   ggTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGt                              >  1:270509/1‑251 (MQ=255)
                                                                                                                                                                                                                      aCCCCGTGGGCAATCTTTAGACCCAACCCCTCAGGCGTCCTTGTACACATTACCATTGAAAATAACCTGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAAAGCTATCCGCTGGTGTTGATCGCTCGATGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCCCCGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAc                           >  1:403981/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccccgggtgcaatgttaacacccaacccaTCGGGCGTTTTTGAAAACACTCACATAAAAAATAACTGGAAAAATTTGGCTCCTTCCCTGGCCGCTCCCCGCTTGCAGAGAGTAATCGTAAACGCTCTCCGGGGGTGTGGAGGGCCCGCTGTTATCGGGGTTGGGGGCCCCCGGTGTTCCGGGGGTGGGACCACCCGGGGCCCCAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAcc                          <  2:89273/232‑1 (MQ=255)
                                                                                                                                                                                                                            tGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACt                     >  1:316019/1‑251 (MQ=255)
                                                                                                                                                                                                                                cAATCTTTAGACCCAAACCATCAGGAGTTCTTGTAAACATTACCATTGAAACTAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCt                 >  2:35940/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aTCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTgg               >  2:52065/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aTCTTTAGACCCAACCCATCAGGCGATCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTg                >  1:347956/1‑250 (MQ=255)
                                                                                                                                                                                                                                              aacccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccg   >  1:5485/1‑251 (MQ=255)
                                                                                                                                                                                                                                               acccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccgt  >  2:308697/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         
CCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGT  >  NC_000913/3468524‑3469013

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