Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,978,166 T→A 21.0% Y245F (TAT→TTT)  lysA ← diaminopimelate decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,978,1660TA21.0% 88.6 / 4.9 46Y245F (TAT→TTT) lysAdiaminopimelate decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base T (19/16);  new base A (2/8);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 7.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACC  >  NC_000913/2977932‑2978400
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
tCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCataata                                                                                                                                                                                                                                        <  1:739386/239‑1 (MQ=255)
    cTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATAAGCCAGACTGCGCTACCCGGGAGCGACCCGGTTCAATTTCCCGTTTCACAGGGTGGCCCAAAAAGCGGGAGGATTGTTCACCCCCGGCGTTTCACAGACCCtaacaaatt                                                                                                                                                                                                                                    >  1:661931/1‑236 (MQ=255)
     tGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGtt                                                                                                                                                                                                                                   <  1:179834/239‑1 (MQ=255)
        cAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCgg                                                                                                                                                                                                                                <  1:1011610/239‑1 (MQ=255)
                ttatgtggtgtTACCTCCCGTTCTTTGCGGGGCGCATCAGATTTTGGAACCCGGCATAACCCAGAACAAAGGGGGGGCCCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGGGCTACCAGGAAGCGCCCCGGTTCAATTTCCAGTTTCACAGGGGGGCCCAAATGGCGGGCGTTTTGCTCCCGCGCGGCATCCCAAAGACCATAATAAGGTGCGGTATCAAcc                                                                                                                                                                                                                        <  2:376965/234‑1 (MQ=255)
                       gtggtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCCGGAAGCGACCCGGTTCAATTTCCAGTTTTACAAGGGGGCCCACATTGCGGGCGATTTGCTCACCCGCGGCCTTTCCCCAACCATAATAATGTTCGGGATCCACccccctt                                                                                                                                                                                                                  >  1:172804/1‑235 (MQ=255)
                            aaCTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATTGCGGGCGAATTTCTCACCCGCGGCATTCCCCCAACCAACATAATATTCGGGATCCAACCCCTCtttcac                                                                                                                                                                                                             >  1:315889/1‑234 (MQ=255)
                                  cGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGACGCGACCCCGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTTCTCCCCCGCGGCCTTTCACAGACCCTAATAATTTTCGGGATCACACGCCCCTTCAccccctttg                                                                                                                                                                                                      >  1:187825/1‑233 (MQ=255)
                                    tACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGGGGCCCCAATGGCGGGCGAGTTTCTCACGCGCGGCATTCCCCAGACCCTAATAAAGTTCGGTGTCAACCCCCCCTTTAcccctttgtat                                                                                                                                                                                                    >  1:678883/1‑231 (MQ=255)
                                         tGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGtt                                                                                                                                                                                                                                   <  1:72148/203‑1 (MQ=255)
                                                                 cccGGCACCAACCAGCACAAAGGGGCCGCTCCCAATTTGGTGGACCCCCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                                                       <  1:297543/239‑1 (MQ=255)
                                                                 caccgCATCACCCACAAAAAGGGGGGGCCTCCCCTTTTGTTTGAGGCTCGGACACTGATTACTTAACCCGCCAGCCGGCGCACCCGGAGAGCGACCGGGTTCATTTCCGGGTTTAAAGGGGTGGCCTAAATGGCGGGCGTTTTGCCCCCGCGCGGCTTCCCACAGACAATAATAATGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCCCCGCCCGCAGaaa                                                                                                                                                                       <  2:661931/235‑1 (MQ=255)
                                                                    ggcaccccccAGCCCAAAGTGGCGCCCCCCATTTTGTTGGCGGCCCGCCACATATGAAATTAAGCCGCCAGTCGGGCCACCAGGGAAGCGCCCCGTTTAAATTCCCAGTTCCCCGGGGGGCCCAAAATGGCGGCCGTTTGCCTCCCGCGCGGCTTCCCACAGACAATAATAATGTCCGTTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAg                                                                                                                                                                    <  2:644915/231‑1 (MQ=255)
                                                                    ggCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCCACGCACCTGAGTAATTAATAAGCCAGACTGCGCTACCAGGAAGCGACCCGGCTCAATTTCCAGTTTTACAGGGGGGCCCCAAATGCGGGCGATTTTGCCGCCCGCGGCATTCCCCACACCAAAAAAAAGTTCGGTCTCAACACCCCCTTCACCCCGTTGATAAGAAACAAAAAAACCACCCCCCCCACAAATa                                                                                                                                                                     >  2:118239/1‑238 (MQ=255)
                                                                                                     ttttttggccccTCCCCCCCGGAGTAATTAATCCGCCGGGCTGCGCTCCCAGGAAGCGACCCGGTTCAATTCCCAGTTTCACAGGGTGGCCAAAATGGGGGGCGATTTCCTCCCGCGGGCAATTCCAAGACCCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGa                                                                                                                                    <  2:1014075/228‑1 (MQ=255)
                                                                                                         tttgACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGTCGATTTGCTCACCCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGc                                                                                                                                 >  1:761320/1‑238 (MQ=255)
                                                                                                            gccGCTCCGCCCCTGAGTAATTAAGCCGCCAGACTGCCCTACCAGGAAGCGCCCCGGTTCAATTTCCGGTTTCACGGGGTGGCCCAAAGGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCa                                                                                                                             <  1:503129/237‑1 (MQ=255)
                                                                                                             aCGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCCCCAt                                                                                                                            >  1:1079525/1‑239 (MQ=255)
                                                                                                              cGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGaaa                                                                                                                                                                       >  1:93959/1‑195 (MQ=255)
                                                                                                                     caaacgcgtatatttatCCGCCAGCCGGCGCTAGCAGGAGGCGCCCGTGTCAAACTCCCCTGTCCCAAGGGGGCCCAAATGGGCGGGCAGTTGGCCCACGCCGGCAATTCCCAAGCCCAGTATAATGGCCGGGAGCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGcacca                                                                                                                     <  2:591351/228‑1 (MQ=255)
                                                                                                                            aGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACcacaca                                                                                                                 >  1:1093246/1‑236 (MQ=255)
                                                                                                                               aattaacCCCCCAGGCTGCGCTCCCAGGGAGCGCCCGTTTCTAATTCCCCGTTCCCAGGGTGGCCCAAAGGGGCGGGCGTTTTGCTCCCCGGCGCCTTCCCCAAGCCCATAAAATTGTCCGGTAACACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGtt                                                                                                           <  2:231284/239‑1 (MQ=255)
                                                                                                                                   aatACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCGTCAATAGCCCGTAAATCCTGACCGAATTCGCCAACCCGACGCACCATAACACAACACACCCCTTCCaa                                                                                                       >  2:34279/1‑238 (MQ=255)
                                                                                                                                   aatACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAAAAAGCCCACCGCCCGCCGAAACAGCCTGTAAATCCTGACCGAATTCGAGGACCTGAAGCACCAAAACACCACAAACCCGGTCCaa                                                                                                       >  2:365347/1‑238 (MQ=255)
                                                                                                                                        gCCAGACTGCGCTACCAGGAAGGGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGAATTCCACAGACCATAATAATGATCGGCATCAACCTCCTCTTCACCCTGTTGATAAGGAACAAAAAGACCACCCCCCGCAGAAATAGCCGGTAAATCCTGACCGAATTCCACGACCTTAAGCACCATAGCACCACACACCTCGTCCCGAttg                                                                                                   >  2:164202/1‑236 (MQ=255)
                                                                                                                                               gcgccccccAGGAGGCGCCCCGGTTAATGTACCAGTTCCCCAGTGTGCACAAAAGGGCGGGCTCTTTGCCCACGGCCGGCATTCCACAAACAAAAAAAATTTGCGGATCCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTTACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATc                                                                                           <  2:761320/232‑1 (MQ=255)
                                                                                                                                                 ggcttcccgggAAGCGGCCCGTTTCAATTTCCAGTTTCCAGGGGTTGCCAAAATGGCGGGCGATTGTCTCACGCGGGGCATTCCAAAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCa                                                                                          <  1:100210/230‑1 (MQ=255)
                                                                                                                                                               gCGACCCGGTTCAATTTCCAGTTTCAGAGGGTGCCCAAAATGGCGGGCGATTTGCTCACGCGGGGCATTCCACAGCCCAAAATAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCa                                                                              <  1:161912/238‑1 (MQ=255)
                                                                                                                                                                  gccgcgTTTAATTTCCCGGTTGCACAGGGTGCCCAAAATGCGGGGCGTTTTGCCCAGGCGGGGCATCCACCAGACCAAAAAAAGGTTGGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAt                                                                            <  2:172804/233‑1 (MQ=255)
                                                                                                                                                                  aCCGGGTTCAATTCCCATTTTCACAGGGTGGCCCAAAGGGCGGCCGTTTGGCTCCCGCGCGGCATCCCACAGACCAAAAAAATTTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAtgt                                                                          <  1:1086993/239‑1 (MQ=255)
                                                                                                                                                                  aCCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCATACACCTGTTCCAGATGGGCATAAt                                                                                            >  1:660954/1‑221 (MQ=255)
                                                                                                                                                                       gTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACCACCCATAGCACCACATACCTGTTCCAGATGGGCATAATCAAACCCAAAACCacagcgcagg                                                                     >  1:1074225/1‑230 (MQ=255)
                                                                                                                                                                                              ggggCCAAAATGGCGGGCGTTTTGCTCCCCGCCGGCATCCCCAAAACCAAAAAAATGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCa                                                  <  1:402422/234‑1 (MQ=255)
                                                                                                                                                                                                ggCCCAAAGGGGGGGCGATTGCCACACGGCTGGCATCCCCCAGACCATAATAATGTTCGGTATCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatga                                             <  1:494481/239‑1 (MQ=255)
                                                                                                                                                                                                 gCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatgat                                            >  1:220610/1‑239 (MQ=255)
                                                                                                                                                                                                  cccAAATGGCGGGCGATTTGCTCACGCGCGGCCTTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTTTAAATCCTGACCGAATTCGCTGACCTGACGCACCATAGCACCCCACCCCTGGTCCAGATGGGCATAATCAACCCCAGAACCAATGAGCATGTGAATGACGACAAGATTCAgatgaaa                                           >  2:256723/1‑237 (MQ=255)
                                                                                                                                                                                                   caaaaGGGCGGGTGATTTGCTCCCGCGCGGCATCCCACTGACCATAATAAGGTGCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGatgatg                                           <  1:255281/236‑1 (MQ=255)
                                                                                                                                                                                                     aaaGGGCGGGCGATTTCCCCCCGCGGGGCATTCCACAGACCAAAATAAGGTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACg                                        <  1:293236/239‑1 (MQ=255)
                                                                                                                                                                                                       aTGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGt                                                                                                                                                               >  2:358408/1‑118 (MQ=255)
                                                                                                                                                                                                       aTGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGt                                                                                                                                                               <  1:358408/118‑1 (MQ=255)
                                                                                                                                                                                                       aTGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGCCGtt                                      >  1:434318/1‑239 (MQ=255)
                                                                                                                                                                                                                            cgcgcgGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                      <  1:99215/156‑1 (MQ=255)
                                                                                                                                                                                                                            cgcgcgGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGa                                                                                                      >  2:99215/1‑156 (MQ=255)
                                                                                                                                                                                                                                   cttatcccaaGACAATAATAATTTTGGGTATCAACCGCCTCTTCCCCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCa              <  1:520897/229‑1 (MQ=255)
                                                                                                                                                                                                                                         cacaGACCATAATAATGTTCGGTATCAACCGCCTATTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt       <  1:365347/239‑1 (MQ=255)
                                                                                                                                                                                                                                            acagccAAAAAAATGTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt       <  1:327942/233‑1 (MQ=255)
                                                                                                                                                                                                                                              aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGATCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATCATCAACTCCAGAACCAATGTTCTCATCAATCCCGACCAGATGCAGATAATGACGTTGTATCACGTCCAGAGCGGCGGGCAGATCGgtgtccg  >  1:510731/1‑236 (MQ=255)
                                                                                                                                                                                                                                              aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCACCGCCAGATGATAACGTAGAATACCATCCAGTGCGGCGGGCACATCgggggaca  >  2:104129/1‑233 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          
TCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACCTGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACAGACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACC  >  NC_000913/2977932‑2978400

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