Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,457,073 T→A 22.0% D233E (GAT→GAA gspD → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,457,0730TA22.0% 119.2 / 5.3 57D233E (GAT→GAAgspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (11/30);  new base A (12/2);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 2.71e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCG  >  NC_000913/3456854‑3457295
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 
aCCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCGGTGGCCAGCCCTCCTCTCCGCGAAGATTGCGGCGAAAAAGCGaa                                                                                                                                                                                                                            >  2:237341/1‑227 (MQ=255)
aCCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAc                                                                                                                                                                                                                           <  2:46113/228‑1 (MQ=255)
aCCCCCCAACGTTCTTATTCTGACCGGTCGTGCCTCCCCCATTAAAAAAATGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAc                                                                                                                                                                                                                           <  1:110222/228‑1 (MQ=255)
 cccTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAAAGAACTGATTTAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGGCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAcc                                                                                                                                                                                                                          <  1:506018/228‑1 (MQ=255)
   cccccaagTTCGTATTCTGACCGGTGGTGCCTCCCCCATTAATAAACTGATGGAAGTTATAAGGCGCGTTGATGTCATCGGCACAGGGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCGCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCGCCTCTCCGCGAAGATTGTGGCGGTTAAGCGAACCaa                                                                                                                                                                                                                        <  1:442902/221‑1 (MQ=255)
   cTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAACGaa                                                                                                                                                                                                                            >  2:110222/1‑224 (MQ=255)
    ttcAACGGTTTTATTCTGGCCGGTCGTGCCTCCACCATTAATAAACTTATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGGCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCaa                                                                                                                                                                                                                        <  2:126971/225‑1 (MQ=255)
     ccAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTTTCCGCGAAGATTGTTGCGGAAAAACTAACCAAca                                                                                                                                                                                                                      >  2:700789/1‑227 (MQ=255)
            tCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAGACGAACaaaacatctta                                                                                                                                                                                                                >  2:383394/1‑219 (MQ=255)
                aTTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCTTCAGCGGAAGATATCGCCGGGATTCTTAATCAATTAATAAGCGAAACCCACGGAAAAAGCCAAATGCCAACCCCCCTATCCGCGAAGATTGTGGAGGAGaacgagcacaaaccccttcccaca                                                                                                                                                                                                           >  1:662859/1‑206 (MQ=255)
                  tCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAg                                                                                                                                                                                                        <  2:56925/228‑1 (MQ=255)
                      aCCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGAAAAAGCCAGATGCCAGACCTCCTCTCCGCGAAGATTGTCGCGGAAAAGCcaaccaacaccatcatcaccacgcgg                                                                                                                                                                                                     >  1:35269/1‑210 (MQ=255)
                      aCCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCGGGAATATGCGTCAGCGGAAGACCTCGCCGAGATTCTTAATCAATTACACGACGCAAGACACGGTACAGCCCAGATGCACGCCCACCTCTCCCCCAAGATGGTAGCGAAAAAGCAAAcaacacccctaccaacccgggga                                                                                                                                                                                                     >  1:339794/1‑206 (MQ=255)
                       cctgttgggCCTCCACCCGTAATAAAAGGATTGAAATCAAAAAGCGGGTTCATTGCATCGGAACAGAGAAGCAGGAAATTATACAGCTTGAATATGCGTAAGCGCAAGATCTCGCCGAGATTCTTAATCCATTAATCAGCGAAAGCCGCGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTAGCATCAGTGGaa                                                                                                                                                                                                   <  1:121443/220‑2 (MQ=255)
                              tGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCAAGCCCTCATCTCCGCGAGAATTGTGGCGGAAAAGCGAACCAACACCCTTCtcatcacgcggccgggaaa                                                                                                                                                                                            >  1:573236/1‑215 (MQ=255)
                               gCCTCCACCATTAAAAAACTGAATGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGa                                                                                                                                                                                              <  2:573236/224‑1 (MQ=255)
                                ccaccaccaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTACGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAg                                                                                                                                                                                         <  2:300059/225‑1 (MQ=255)
                                     accaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTGATCGGCACAGAGAAGCACCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATAAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGCGGCGGCAAACCGAACCCACTCCTTTATCATCCATGGCCCGGACAAAGCag                                                                                                                                                                                      >  2:703292/1‑225 (MQ=255)
                                     accaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAGCGAACCAACCCCCTTCTCTCCCCTGGACCCGaaaaacacac                                                                                                                                                                                     >  1:225497/1‑222 (MQ=255)
                                     accaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCCCGAAGATTGTGGCGGAAAAACGAAACACCACACTTATCACCCCTGGACCCGGaaaaaaaac                                                                                                                                                                                      >  1:389541/1‑223 (MQ=255)
                                         ttAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTACAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTATGCCGGAGAAGCGTACACACTCTCTTATCATCATGGGACCCGGAAAAAGCACGcagc                                                                                                                                                                                 >  2:455359/1‑225 (MQ=255)
                                         ttAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAg                                                                                                                                                                                 <  2:389541/227‑1 (MQ=255)
                                         ttAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAGCGAACCAAATCCCTTTTCACCCAGCGACCCGGaaaaaaaaccccac                                                                                                                                                                                 >  1:82236/1‑219 (MQ=255)
                                             taaACTGATTGAAGTCATAAATCGCAGTGATGTCAGCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTTGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                                                                                            <  2:339794/228‑1 (MQ=255)
                                             taaACTAATTGAAGTCATAAAGCGCGTTGCGTCCATCGGCCCAGAGAAGCAACAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                                                                                            <  2:662859/228‑1 (MQ=255)
                                             taaAATGATTGAAATTAGAAAGGGCGGTGGTGCTTACGGCACAGAGAAAGCAGAAATTATTCCACCGGGAATAGGCTCAGCGGAAGAATCCGCCGGAATTTTTAACCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTAGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCAt                                                                                                                                                                            <  1:53979/228‑1 (MQ=255)
                                                aCTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGAGGCGGATAAGCGAAACCACCCTCTTATCATCATGGGCCCGGAAAAAGAACGCCAgagactcac                                                                                                                                                                         >  2:95279/1‑220 (MQ=255)
                                                 cTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCAAACCAAAACCACTATCATCAAGGGACCGGGAAAAGCACACCAggcgagtcat                                                                                                                                                                        >  2:369813/1‑219 (MQ=255)
                                                 cTGATGGAAGTCAGAAAGCTCGTTGATTTAATCGGCACAGAGAAGCAGCAAATTATCCATCTGGAATATGCGTCAGCGGAAGATCCCGCCGAGATTCTTAAGCAATTAACCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCGCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCa                                                                                                                                                                          <  2:606868/226‑1 (MQ=255)
                                                    aTTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAACAACCCCCTTATCCTCATTGGACCCGAAAAAGCACGCCAGCGCATACCTTc                                                                                                                                                                      >  2:128881/1‑227 (MQ=255)
                                                                    gtgTTGAGTTCATCGGCCAAGAGAAACTAGAAATTATTCATCTGGAATATGCGTAAGCGGAAGATCTCGCCAAGATTCTTAATCAGTTAACCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCAGCCCTTCATTACTGAAAAGCCTTg                                                                                                                                                     <  1:299640/226‑1 (MQ=255)
                                                                     cgTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAAACCTCTTATCATCAGTGGACCGGGAAAAGCACGCCAGCGCATCACTTCATTACTGTAAAACCttaa                                                                                                                                                    >  2:553205/1‑226 (MQ=255)
                                                                                tCGGCACAGAGAACCAGCAAATTAATCACCCGGAATACGCGCCAGCGGAAGATCGCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAagaga                                                                                                                                         <  1:383394/228‑1 (MQ=255)
                                                                                   gCACAGAGAAGCAGCAAATTATTCATCTGGAATAGGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                                                                                      <  2:82236/228‑1 (MQ=255)
                                                                                     acaGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGa                                                                                                                                                                                                    <  2:565002/165‑1 (MQ=255)
                                                                                       agagAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACCCTCTTATCATCCGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACCGCAAAACATTGATGTCGAAGAAAACGagga                                                                                                                                  >  2:205563/1‑228 (MQ=255)
                                                                                            aaGCAGCAAATTATTCAACTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGaggaagga                                                                                                                              <  2:225497/228‑1 (MQ=255)
                                                                                               cagcaAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAACGAACCa                                                                                                                                                                                                                         >  1:539585/1‑136 (MQ=255)
                                                                                                agcaAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCaa                                                                                                                                                                                                                        <  2:539585/136‑1 (MQ=255)
                                                                                                            caacaggAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAAACAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGATCCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTTCTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTttat                                                                                                               <  2:592733/223‑1 (MQ=255)
                                                                                                                              gcgtaacactctcccgcaaTTTTTAACCACTTAATCAGCAAAAGTCACGGAAAAAGACATCTGCCCGCCCTCCGATCCGCGAAGAGTGGGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGCCCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGTTTTATTACCTGAAATATGATaa                                                                                              <  1:703292/209‑1 (MQ=255)
                                                                                                                                            cGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                                                                                                     <  2:158607/158‑1 (MQ=255)
                                                                                                                                                                     gAAAGCCACGTTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGTCGGATAAGCGAACCCACTCTCTTATCATCAGTGGACCGGAAAAAGCCCGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATTCCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGaa                                                     <  2:289460/228‑1 (MQ=255)
                                                                                                                                                                      acaGCCACGGTCAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAACGACCCACCCCCCTTTTCCTCCCCGGGACCGACAAAGCACCCCCGCGCCTCCCCTCATTTCTTGAAAACCCTGCTGTCCAAGAAGGCGGGGACGGACATTCCCGGCGTTATTACCTGCAAACTGCTAAAGCCGCGCACTTGGGGGCAGTACTATACCGTGGTTTCCgca                                                     >  1:486887/3‑226 (MQ=255)
                                                                                                                                                                             cGGGAAAAGCCAGATGCCATCCCTCCTCACCGTCAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGGGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGaa                                             <  1:553205/228‑1 (MQ=255)
                                                                                                                                                                                      ccAGATTCCCGCCCTCCTCACCGCGAAGATTGTGGAGGAAAAGCCAACCAACTCTCTTATCACAAGTGGACCGGAAAAAACACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGaaa                                     <  1:700789/227‑1 (MQ=255)
                                                                                                                                                                                          aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGGTTATTACCTGAAAAATGCTAAAGCCACGAATCTAGGGGAAGGGCTAACCGGGGTTTCCCAAAAGATGAAAAGTGAAAAAggg                                >  2:418122/1‑228 (MQ=255)
                                                                                                                                                                                          aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAGCGAACCAACTCTCTTATCACCCGGGGGCCGGAAAAAAGAACCCAGCGCCACACTTCATTTCTTAAAAGACTTGATGGAGCAAGGAACGGGGAAGGAAAAACCCGGGGTTATTACCTCGAATATGTTAAGACCACGAAACAGGTGGGAGGGTCAAACGGGGTTTCCCAAAAgactcaaaaaaaaaaagag                                >  1:572560/1‑210 (MQ=255)
                                                                                                                                                                                           ggCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAg                                           <  2:572560/215‑1 (MQ=255)
                                                                                                                                                                                             ccaagccatctcaCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCAACAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAAt                             <  2:165578/219‑1 (MQ=255)
                                                                                                                                                                                               aGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAAATCTCTTATCATCCGTGGACCGGAAAAAGCACGCCAGCGCATCCCTTCATTACTGAAAAACCTTGGTGTCGAAGAGAGAGAGGAAGGGAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGGATCTAGGTGAAGTGCTAACCGGGGTTTTCGAAAAGCTGAAAGATGAAAAAGGGAATg                            >  2:213694/1‑227 (MQ=255)
                                                                                                                                                                                                      ctctCCGCGAAGAGTGTGGCGGAAAAGCGAAACAACTCTCTTAGCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACCTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAGCCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAg                    <  2:129231/228‑1 (MQ=255)
                                                                                                                                                                                                           cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACCTCATTACTGAAAAGCCCTGATGTCGAAGAGAGCGAGGAAGGAAAAACCCCGGGTTATTTCCTGAAATAAGCTAAAGCCACGAATATGGGGGAAGGGCTAAACGGGGTTTCCCAAAAGCTGCAAAAGAAAAAAgagcgggtccttataaccac               >  1:247378/1‑209 (MQ=255)
                                                                                                                                                                                                            gcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAg                                           <  2:247378/199‑1 (MQ=255)
                                                                                                                                                                                                             cgAAGATTGTGGCGGATAAGCGAACCAACTCCCTTATCCTCCGTGGACCGGAAAAAGCACCCCCAGGACTCCCTTCTTTACTGCAAAGCCCTGGTGTCGAAAAAAACACGGAAAGAAATACCCCGGTTTATTACCTGAAAAAATGTAAAGACAAGACAATGGGGGAAGTGCTAACCGCGGTTTTCTGAAAGACCAAAAATGAAACAGGGAAAGGGCGTGAGccccccc             >  1:405873/1‑227 (MQ=255)
                                                                                                                                                                                                                      tgGCGGATAAGCGACCCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGAGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGgcg    <  2:35269/228‑1 (MQ=255)
                                                                                                                                                                                                                        gCGGATAAGCGAACCAACTCTCTTATCATCAGTGTACCGGAAAAAGCACGCCAGAGCATCACTACATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGcgcg  <  2:405873/228‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 
ACCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAGAAGCAGCAAATTATTCATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCATCAGTGGACCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTAACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCG  >  NC_000913/3456854‑3457295

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: