Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 429,554 A→T 5.4% Y17F (TAT→TTT)  secF → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913429,5540AT 5.4% 404.8 / 9.1 140Y17F (TAT→TTT) secFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base A (78/52);  new base T (3/5);  total (81/59)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 141 total aligned reads displayed.

GCGACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGC  >  NC_000913/429405‑429701
                                                                                                                                                      |                                                                                                                                                   
gCGACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTa                                                                                                                                                     >  1:654085/1‑150 (MQ=255)
    cGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGAc                                                                                                                                                 <  1:1510232/150‑1 (MQ=255)
    cGTAGATGTTTACAGAGATTGTCGGTAAGCGTGCCATCGTAAACCTGCTATATGGAGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGTCACAGGAATATACTGTTG‑AACAACTAAACCCCGGACGTAAAGTCTATGAc                                                                                                                                                 >  2:439814/1‑150 (MQ=255)
     ggcGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACt                                                                                                                                                <  1:821571/148‑1 (MQ=255)
       cGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACttt                                                                                                                                              <  2:1589588/150‑1 (MQ=255)
       cGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTATATGACttt                                                                                                                                              >  1:1807830/1‑150 (MQ=255)
        gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTa                                                                                                                                             >  1:1539951/1‑150 (MQ=255)
        gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTa                                                                                                                                             >  2:49904/1‑150 (MQ=255)
         aTGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTAt                                                                                                                                            >  2:143576/1‑150 (MQ=255)
           gTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCATAGACTTTacgg                                                                                                                                          <  2:1946287/150‑4 (MQ=255)
           gTTTACCGAGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCATAGACTTTacgg                                                                                                                                          <  2:538078/150‑4 (MQ=255)
            tttACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCATAGACTTTacggc                                                                                                                                         <  1:1907648/150‑5 (MQ=255)
             ttACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATgcgc                                                                                                                                        >  2:744351/1‑150 (MQ=255)
              tACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCt                                                                                                                                       >  1:904267/1‑150 (MQ=255)
                ccGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGCAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGAGCTgg                                                                                                                                     >  1:325879/1‑150 (MQ=255)
                ccGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTgg                                                                                                                                     <  2:239494/150‑1 (MQ=255)
                ccGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTgg                                                                                                                                     <  2:214729/150‑1 (MQ=255)
                ccGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTgg                                                                                                                                     >  2:1335933/1‑150 (MQ=255)
                ccGCGAATGTCGGTACGCGTGCCATCGTAAAACTCCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGACGAGTGCGATGTGGCACAGGAAAATACTGTTG‑AACAACTAAAACACGGCCGTAAAGTCTATGACTTTATGAGCTgg                                                                                                                                     >  1:342042/1‑150 (MQ=255)
                   cgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGac                                                                                                                                  <  2:480066/150‑1 (MQ=255)
                   cgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGac                                                                                                                                  <  1:1488597/150‑1 (MQ=255)
                    gATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATTTACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATCACTGGATGCGCTGGGact                                                                                                                                 >  1:1143924/1‑150 (MQ=255)
                        gTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCTAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTg                                                                                                                             >  1:72555/1‑150 (MQ=255)
                        gTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTg                                                                                                                             >  1:1372735/1‑150 (MQ=255)
                        gTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTg                                                                                                                             >  1:1104804/1‑150 (MQ=255)
                           ggTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTCCTGGGc                                                                                                                          >  1:386696/1‑150 (MQ=255)
                            gTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑ATCAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCt                                                                                                                         >  2:426774/1‑150 (MQ=255)
                               cgcgTGCCTTCTTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAATCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTc                                                                                                                      <  2:124532/150‑1 (MQ=255)
                                 cgTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCgg                                                                                                                    >  2:743241/1‑150 (MQ=255)
                                        tCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATctct                                                                                                             >  1:1106453/1‑150 (MQ=255)
                                           tAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGt                                                                                                          >  2:578680/1‑150 (MQ=255)
                                            aaaCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTc                                                                                                         >  2:1670000/1‑150 (MQ=255)
                                              aCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACACCTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctg                                                                                                       >  1:1052457/1‑150 (MQ=255)
                                              aCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctg                                                                                                       >  1:1012408/1‑150 (MQ=255)
                                                 tGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctgctg                                                                                                    <  1:1326680/150‑1 (MQ=255)
                                                 tGCTATATGGCGGCAAGCGAGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGGTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATTCGATGGGTCTACTGGGCTTTCGGCATCTCTGGTctgctg                                                                                                    >  1:279390/1‑150 (MQ=255)
                                                  gCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGt                                                                                                   <  2:562464/150‑1 (MQ=255)
                                                  gCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGt                                                                                                   >  2:1852680/1‑150 (MQ=255)
                                                   cTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCATAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                  <  1:1397445/150‑1 (MQ=255)
                                                   cTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGGGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                  <  1:525262/150‑1 (MQ=255)
                                                       atGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGTAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATc                                                                                              <  1:1598397/150‑1 (MQ=255)
                                                       atGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTTTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATc                                                                                              <  1:957090/150‑1 (MQ=255)
                                                        tGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCg                                                                                             >  2:265674/1‑150 (MQ=255)
                                                         ggcggcACGCGCGTCAAGAATCTGTCAAGCTGAGGAGTGAGATGTGTAACAGGAATATAATGTTT‑AACAAATAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTTCTGGGCTTTCGTCATCTCTGGTCTGCTGTTAATCgc                                                                                            <  2:324669/150‑1 (MQ=255)
                                                         ggcggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgc                                                                                            >  2:63324/1‑150 (MQ=255)
                                                         ggcggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgc                                                                                            <  2:1950672/150‑1 (MQ=255)
                                                         ggcggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgc                                                                                            >  2:69334/1‑150 (MQ=255)
                                                          gcggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgct                                                                                           >  1:1643712/1‑150 (MQ=255)
                                                           cggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAAGATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTTATCgcgg                                                                                          >  1:371807/1‑148 (MQ=255)
                                                           cggcAAGCGCGGCAAGAAGCTGGCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctg                                                                                          <  1:1313287/150‑1 (MQ=255)
                                                            ggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                         <  1:869590/150‑1 (MQ=255)
                                                                     gTCAAGAAGCTGTCAATCTGAGGAGTGCGATGAGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGttat                                                                                >  1:1396634/1‑150 (MQ=255)
                                                                        aagaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTtattat                                                                             <  1:564397/150‑1 (MQ=255)
                                                                         agaagCTTTCAATCTGAGTAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCTGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATg                                                                            <  1:1743510/150‑1 (MQ=255)
                                                                         agaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATg                                                                            <  2:823346/150‑1 (MQ=255)
                                                                         agaagCTGTCAAGCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATg                                                                            <  1:1056001/150‑1 (MQ=255)
                                                                             gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCg                                                                        >  1:150499/1‑150 (MQ=255)
                                                                              cTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGACATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGt                                                                       >  1:279569/1‑150 (MQ=255)
                                                                                  cAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTgcgc                                                                   >  2:896698/1‑150 (MQ=255)
                                                                                  cAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTgcgc                                                                   >  1:1481539/1‑150 (MQ=255)
                                                                                     tCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAATCCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGc                                                                <  1:1459548/150‑1 (MQ=255)
                                                                                        gaggagTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCttt                                                             >  2:75168/1‑150 (MQ=255)
                                                                                        gaggagTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCttt                                                             >  1:1293545/1‑150 (MQ=255)
                                                                                         aggagTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTa                                                            <  2:686075/150‑1 (MQ=255)
                                                                                            tgTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACt                                                         <  2:763616/149‑1 (MQ=255)
                                                                                            agTGCGATGTGGCACAGGAATATACTGTTGAAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAAc                                                          <  1:714250/150‑1 (MQ=255)
                                                                                            agTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTTAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACt                                                         <  1:143576/150‑1 (MQ=255)
                                                                                              tGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTgg                                                       >  2:1862031/1‑150 (MQ=255)
                                                                                                cGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTgggg                                                     >  1:406140/1‑150 (MQ=255)
                                                                                                  aTGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGc                                                   >  2:137056/1‑150 (MQ=255)
                                                                                                     tgGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGGCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                <  2:925424/150‑1 (MQ=255)
                                                                                                       gCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGAt                                              >  2:841008/1‑150 (MQ=255)
                                                                                                        cacaGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGAtt                                             >  2:1803274/1‑150 (MQ=255)
                                                                                                           aGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCa                                          >  1:142018/1‑150 (MQ=255)
                                                                                                            ggAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCAc                                         >  2:1477992/1‑150 (MQ=255)
                                                                                                             gAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTGCGGCATCTCTGGGCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCAcc                                        <  2:1774168/150‑1 (MQ=255)
                                                                                                                 ataCTGTTG‑AACAACTAAACCACGGCCGTACAGGCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGGTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtg                                    <  2:891383/150‑1 (MQ=255)
                                                                                                                 ataCTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtg                                    <  1:557618/150‑1 (MQ=255)
                                                                                                                  taCTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGCTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                   <  2:933913/150‑1 (MQ=255)
                                                                                                                    cTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTa                                 >  1:1911294/1‑150 (MQ=255)
                                                                                                                       ttG‑AACAACTAAACCACGGCCGTAAAGTCTCTGACTTTATGCGCGGGGACTACTGGGCTTTCGGCTTCTCTGGTCTGCTGTTAATCGCTGCTATCGTTTTTATGGGCGTGCGCGGCTTTAACTGGGGGCGGGATTTCACCGGTGGTACgg                              <  1:159101/150‑1 (MQ=255)
                                                                                                                       ttG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACgg                              <  2:1372735/150‑1 (MQ=255)
                                                                                                                       ttG‑AACAACTAAACCACGGCCCTACAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGGATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACgg                              <  2:1490578/150‑1 (MQ=255)
                                                                                                                        tG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGt                             >  2:794314/1‑150 (MQ=255)
                                                                                                                        tG‑AACAACTAAACCACGGCCGTAAAGGCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTAATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGt                             <  2:169402/150‑1 (MQ=255)
                                                                                                                         g‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGtt                            >  2:1640527/1‑150 (MQ=255)
                                                                                                                           caacaatAACCCACGCCCGTCAAGTCTCTCTCTTAATGCGCGTGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAAGCGCTGCTATCGTTTTTATGGGCGGGCGCGGCTTTAACTGGGGGGTTCATTTCCCCGGTGGTACGGTTa                           <  2:329584/144‑1 (MQ=255)
                                                                                                                           aacaacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTa                           >  2:1016090/1‑150 (MQ=255)
                                                                                                                            acaacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAt                          <  2:1552609/150‑1 (MQ=255)
                                                                                                                             caacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAtt                         >  2:1066210/1‑150 (MQ=255)
                                                                                                                             caacTAAACCACGGCCGTAAAGGCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAtt                         <  2:1440232/150‑1 (MQ=255)
                                                                                                                               acTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGa                       >  1:185935/1‑150 (MQ=255)
                                                                                                                                  aaaCCACGGCCGTAACCTCTATCACTTTATGCGCTGGGACTACTGGGCTTGCGGCAGCTCTGGGCTGCTGTGAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAAt                    <  2:325879/150‑1 (MQ=255)
                                                                                                                                    aCCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTa                  <  1:49904/150‑1 (MQ=255)
                                                                                                                                      cACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACg                >  2:1657967/1‑150 (MQ=255)
                                                                                                                                         gcgcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTTACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTc             >  1:1281034/6‑150 (MQ=255)
                                                                                                                                          gCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCg            >  1:1669771/1‑150 (MQ=255)
                                                                                                                                          cgcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTCGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCg            >  2:1271302/5‑150 (MQ=255)
                                                                                                                                           gcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGa           >  1:533800/4‑150 (MQ=255)
                                                                                                                                                   tCTATGACTTTATGAGCTGGGACTACTGGGCTTTAGGCATCTCTGGTCTGCTGTTAATAGCTTCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGTGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCgg   <  1:827813/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCGACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTG‑AACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGC  >  NC_000913/429405‑429701

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