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→GAA) gspDgeneral 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

ACCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCG  >  NC_000913/3456854‑3457295
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 
aCCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCGGTGGCCAGCCCTCCTCTCCGCGAAGATTGCGGCGAAAAAGCGaa                                                                                                                                                                                                                            >  2:237341/1‑227 (MQ=255)
aCCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAc                                                                                                                                                                                                                           <  2:46113/228‑1 (MQ=255)
aCCCCCCAACGTTCTTATTCTGACCGGTCGTGCCTCCCCCATTAAAAAAATGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAc                                                                                                                                                                                                                           <  1:110222/228‑1 (MQ=255)
 cccTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAAAGAACTGATTTAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGGCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAcc                                                                                                                                                                                                                          <  1:506018/228‑1 (MQ=255)
   cccccaagTTCGTATTCTGACCGGTGGTGCCTCCCCCATTAATAAACTGATGGAAGTTATAAGGCGCGTTGATGTCATCGGCACAGGG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCGCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCGCCTCTCCGCGAAGATTGTGGCGGTTAAGCGAACCaa                                                                                                                                                                                                                        <  1:442902/221‑1 (MQ=255)
   cTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAACGaa                                                                                                                                                                                                                            >  2:110222/1‑224 (MQ=255)
    ttcAACGGTTTTATTCTGGCCGGTCGTGCCTCCACCATTAATAAACTTATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGGCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCaa                                                                                                                                                                                                                        <  2:126971/225‑1 (MQ=255)
     ccAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTTTCCGCGAAGATTGTTGCGGAAAAACTAACCAAca                                                                                                                                                                                                                      >  2:700789/1‑227 (MQ=255)
            tCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAGACGAACaaaacatctta                                                                                                                                                                                                                >  2:383394/1‑219 (MQ=255)
                aTTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCTTCAGCGGAAGATATCGCCGGGATTCTTAATCAATTAATAAGCGAAACCCACGGAAAAAGCCAAATGCCAACCCCCCTATCCGCGAAGATTGTGGAGGAGaacgagcacaaaccccttcccaca                                                                                                                                                                                                           >  1:662859/1‑206 (MQ=255)
                  tCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAg                                                                                                                                                                                                        <  2:56925/228‑1 (MQ=255)
                      aCCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGAAAAAGCCAGATGCCAGACCTCCTCTCCGCGAAGATTGTCGCGGAAAAGCcaaccaacaccatcatcaccacgcgg                                                                                                                                                                                                     >  1:35269/1‑210 (MQ=255)
                      aCCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCGGGAATATGCGTCAGCGGAAGACCTCGCCGAGATTCTTAATCAATTACACGACGCAAGACACGGTACAGCCCAGATGCACGCCCACCTCTCCCCCAAGATGGTAGCGAAAAAGCAAAcaacacccctaccaacccgggga                                                                                                                                                                                                     >  1:339794/1‑206 (MQ=255)
                       cctgttgggCCTCCACCCGTAATAAAAGGATTGAAATCAAAAAGCGGGTTCATTGCATCGGAACAGAG‑AAGCAGGAAATTATA‑CAGCTTGAATATGCGTAAGCGCAAGATCTCGCCGAGATTCTTAATCCATTAATCAGCGAAAGCCGCGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTAGCA‑TCAGTGGaa                                                                                                                                                                                                   <  1:121443/220‑2 (MQ=255)
                              tGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCAAGCCCTCATCTCCGCGAGAATTGTGGCGGAAAAGCGAACCAACACCCTTCtca‑tcacgcggccgggaaa                                                                                                                                                                                            >  1:573236/1‑215 (MQ=255)
                               gCCTCCACCATTAAAAAACTGAATGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGa                                                                                                                                                                                              <  2:573236/224‑1 (MQ=255)
                                ccaccaccaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTACGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAg                                                                                                                                                                                         <  2:300059/225‑1 (MQ=255)
                                     accaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTGATCGGCACAGAG‑AAGCACCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATAAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGCGGCGGCAAACCGAACCCACTCCTTTATCA‑TCCATGGC‑CCGGACAAAGCag                                                                                                                                                                                      >  2:703292/1‑225 (MQ=255)
                                     accaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAGCGAACCAACCCCCTTCTCT‑CCCCTGGA‑CCCGaaaaacacac                                                                                                                                                                                     >  1:225497/1‑222 (MQ=255)
                                     accaTTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCCCGAAGATTGTGGCGGAAAAACGAAACACCACACTTATCA‑CCCCTGGACCCGGaaaaaaaac                                                                                                                                                                                      >  1:389541/1‑223 (MQ=255)
                                         ttAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTACAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTATGCCGGAGAAGCGTACACACTCTCTTATCA‑TCATGGGACCCGGAAAAAGCACGcagc                                                                                                                                                                                 >  2:455359/1‑225 (MQ=255)
                                         ttAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAg                                                                                                                                                                                 <  2:389541/227‑1 (MQ=255)
                                         ttAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAGCGAACCAAATCCCTTTTCACCCAGCGAC‑CCGGaaaaaaaaccccac                                                                                                                                                                                 >  1:82236/1‑219 (MQ=255)
                                             taaACTGATTGAAGTCATAAATCGCAGTGATGTCAGCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTTGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCAt                                                                                                                                                                            <  2:339794/228‑1 (MQ=255)
                                             taaACTAATTGAAGTCATAAAGCGCGTTGCGTCCATCGGCCCAGAG‑AAGCAACAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCAt                                                                                                                                                                            <  2:662859/228‑1 (MQ=255)
                                             taaAATGATTGAAATTAGAAAGGGCGGTGGTGCTTACGGCACAGAGAAAGCAG‑AAATTATTCCACCGGGAATAGGC‑TCAGCGGAAGAATCCGCCGGAATTTTTAACCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTAGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCAt                                                                                                                                                                            <  1:53979/228‑1 (MQ=255)
                                                aCTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGAGGCGGATAAGCGAAACCACCCTCTTATCA‑TCATGGGC‑CCGGAAAAAGAACGCCAgagactcac                                                                                                                                                                         >  2:95279/1‑220 (MQ=255)
                                                 cTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCAAACCAAAACCACTATCA‑TCAAGGGA‑CCGGGAAAAGCACACCAggcgagtcat                                                                                                                                                                        >  2:369813/1‑219 (MQ=255)
                                                 cTGATGGAAGTCAGAAAGCTCGTTGATTTAATCGGCACAGAG‑AAGCAGCAAATTATC‑CATCTGGAATATGCGTCAGCGGAAGATCCCGCCGAGATTCTTAAGCAATTAACCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCGCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCa                                                                                                                                                                          <  2:606868/226‑1 (MQ=255)
                                                    aTTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAAACAACCCCCTTATCC‑TCATTGGA‑CCCGAAAAAGCACGCCAGCGCATACCTTc                                                                                                                                                                      >  2:128881/1‑227 (MQ=255)
                                                                    gtgTTGAGTTCATCGGCCAAGAG‑AAACTAGAAATTATT‑CATCTGGAATATGCGTAAGCGGAAGATCTCGCCAAGATTCTTAATCAGTTAACCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCAGCCCTTCATTACTGAAAAGCCTTg                                                                                                                                                     <  1:299640/226‑1 (MQ=255)
                                                                     cgTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAAACCTCTTATCA‑TCAGTGGA‑CCGGGAAAAGCACGCCAGCGCATCACTTCATTACTGTAAAACCttaa                                                                                                                                                    >  2:553205/1‑226 (MQ=255)
                                                                                tCGGCACAGAG‑AACCAGCAAATTAAT‑CACCCGGAATACGCGCCAGCGGAAGATCGCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAagaga                                                                                                                                         <  1:383394/228‑1 (MQ=255)
                                                                                   gCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATAGGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCg                                                                                                                                      <  2:82236/228‑1 (MQ=255)
                                                                                     acaGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGa                                                                                                                                                                                                    <  2:565002/165‑1 (MQ=255)
                                                                                       agag‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACCCTCTTATCA‑TCCGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACCGCAAAACATTGATGTCGAAGAAAACGagga                                                                                                                                  >  2:205563/1‑228 (MQ=255)
                                                                                            aaGCAGCAAATTATT‑CAACTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGaggaagga                                                                                                                              <  2:225497/228‑1 (MQ=255)
                                                                                               cagcaAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAACGAACCa                                                                                                                                                                                                                         >  1:539585/1‑136 (MQ=255)
                                                                                                agcaAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCaa                                                                                                                                                                                                                        <  2:539585/136‑1 (MQ=255)
                                                                                                            caacaggAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAAACAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGATCCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTTCTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTttat                                                                                                               <  2:592733/223‑1 (MQ=255)
                                                                                                                              gcgtaacactctcccgcaaTTTTTAACCACTTAATCAGCAAAAGTCACGGAAAAAGACATCTGCCCGCCCTCCGATCCGCGAAGAGTGGGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGC‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGTTTTATTACCTGAAATATGATaa                                                                                              <  1:703292/209‑1 (MQ=255)
                                                                                                                                            cGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTg                                                                                                                                                     <  2:158607/158‑1 (MQ=255)
                                                                                                                                                                     gAAAGCCACGTTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGTCGGATAAGCGAACCCACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCCCGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATTCCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGaa                                                     <  2:289460/228‑1 (MQ=255)
                                                                                                                                                                      acaGCCACGGTCAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAACGACCCACCCCCCTTTTCC‑TCCCCGGG‑ACCGACAAAGCACCCCCGCGCCTCCCCTCATTTCTTGAAAACCCTGCTGTCCAAGAAGGCGGGGACGGACATTCCCGGCGTTATTACCTGCAAACTGCTAAAGCCGCGCACTTGGGGGCAGTACTATACCGTGGTTTCCgca                                                     >  1:486887/3‑226 (MQ=255)
                                                                                                                                                                             cGGGAAAAGCCAGATGCCATCCCTCCTCACCGTCAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGGGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGaa                                             <  1:553205/228‑1 (MQ=255)
                                                                                                                                                                                      ccAGATTCCCGCCCTCCTCACCGCGAAGATTGTGGAGGAAAAGCCAACCAACTCTCTTATCA‑CAAGTGGA‑CCGGAAAAAACACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAGATGaaa                                     <  1:700789/227‑1 (MQ=255)
                                                                                                                                                                                          aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGGTTATTACCTGAAAAATGCTAAAGCCACGAATCTAGGGGAAGGGCTA‑ACCGGGGTTTCCCAAAAGATGAAAAGTGAAAAAggg                                >  2:418122/1‑228 (MQ=255)
                                                                                                                                                                                          aTGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGAAAAGCGAACCAACTCTCTTATCA‑CCCGGGGG‑CCGGAAAAAAGAACCCAGCGCCACACTTCATTTCTTAAAAGACTTGATGGAGCAAGGAACGGGGAAGGAAAAACCCGGGGTTATTACCTCGAATATGTTAAGACCACGAAACAGGTGGGAGGGTCA‑AACGGGGTTTCCCAAAAgactcaaaaaaaaaaagag                                >  1:572560/1‑210 (MQ=255)
                                                                                                                                                                                           ggCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAg                                           <  2:572560/215‑1 (MQ=255)
                                                                                                                                                                                             ccaagccatctcaCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑ACAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAAt                             <  2:165578/219‑1 (MQ=255)
                                                                                                                                                                                               aGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAAATCTCTTATCA‑TCCGTGGA‑CCGGAAAAAGCACGCCAGCGCATCCCTTCATTACTGAAAAACCTTGGTGTCGAAGAGAGAGAGGAAGGGAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGGATCTAGGTGAAGTGCTA‑ACCGGGGTTTTCGAAAAGCTGAAAGATGAAAAAGGGAATg                            >  2:213694/1‑227 (MQ=255)
                                                                                                                                                                                                      ctctCCGCGAAGAGTGTGGCGGAAAAGCGAAACAACTCTCTTAGCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACCTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑GCCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAg                    <  2:129231/228‑1 (MQ=255)
                                                                                                                                                                                                           cgcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACCTCATTACTGAAAAGCCCTGATGTCGAAGAGAGCGAGGAAGGAAAAACCCCGGGTTATTTCCTGAAATAAGCTAAAGCCACGAATATGGGGGAAGGGCTA‑AACGGGGTTTCCCAAAAGCTGCAAAAGAAAAAAgagcgggtccttataaccac               >  1:247378/1‑209 (MQ=255)
                                                                                                                                                                                                            gcgAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAg                                           <  2:247378/199‑1 (MQ=255)
                                                                                                                                                                                                             cgAAGATTGTGGCGGATAAGCGAACCAACTCCCTTATCC‑TCCGTGGA‑CCGGAAAAAGCACCCCCAGGACTCCCTTCTTTACTGCAAAGCCCTGGTGTCGAAAAAAACACGGAAAGAAATACCCCGGTTTATTACCTGAAAAAATGTAAAGACAAGACAATGGGGGAAGTGCTA‑ACCGCGGTTTTCTGAAAGACCAAAAATGAAACAGGGAAAGGGCGTGAGccccccc             >  1:405873/1‑227 (MQ=255)
                                                                                                                                                                                                                      tgGCGGATAAGCGACCCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGAGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGgcg    <  2:35269/228‑1 (MQ=255)
                                                                                                                                                                                                                        gCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGTA‑CCGGAAAAAGCACGCCAGAGCATCACTACATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGcgcg  <  2:405873/228‑1 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                 
ACCCTCCAACGTTCTTATTCTGACCGGTCGTGCCTCCACCATTAATAAACTGATTGAAGTCATAAAGCGCGTTGATGTCATCGGCACAGAG‑AAGCAGCAAATTATT‑CATCTGGAATATGCGTCAGCGGAAGATCTCGCCGAGATTCTTAATCAATTAATCAGCGAAAGCCACGGTAAAAGCCAGATGCCAGCCCTCCTCTCCGCGAAGATTGTGGCGGATAAGCGAACCAACTCTCTTATCA‑TCAGTGGA‑CCGGAAAAAGCACGCCAGCGCATCACTTCATTACTGAAAAGCCTTGATGTCGAAGAGAGCGAGGAAGGAAATACCCGGGTTTATTACCTGAAATATGCTAAAGCCACGAATCTGGTGGAAGTGCTA‑ACCGGTGTTTCCGAAAAGCTGAAAGATGAAAAAGGGAATGCGCGTAAGCCCTCCTCTTCTGGCGCG  >  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: ‑