Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 429,554 A→T 5.6% 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.6% 260.3 / 2.4 94Y17F (TAT→TTT) secFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base A (48/39);  new base T (2/4);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCT  >  NC_000913/429409‑429702
                                                                                                                                                 |                                                                                                                                                    
cctCGATGTGAACCGCTAGTGTCGGTACGCGTGCCAGCGAACCCCTGCTATCTGGCGGCACGCGCGTCAAGCCGCGGTCAATCTGAGGCGTGCGAGGTGGCACAGGAATATACTGTTGACCAACTAAACCACGGCCGTAAAGTCTATGa                                                                                                                                                   <  2:1523356/147‑1 (MQ=255)
  tCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGAGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACt                                                                                                                                                 <  2:1169291/149‑1 (MQ=255)
   cGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTATAGTCTATGACtt                                                                                                                                                >  2:1480147/1‑149 (MQ=255)
    gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGGGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACttt                                                                                                                                               <  2:1495637/149‑1 (MQ=255)
    gATGTTTACCGCGATTGGCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGGCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACttt                                                                                                                                               <  2:1465079/149‑1 (MQ=255)
      tGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGGCTTTAt                                                                                                                                             >  2:314354/1‑149 (MQ=255)
      tGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTAt                                                                                                                                             <  2:521600/149‑1 (MQ=255)
       gTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATg                                                                                                                                            >  2:780643/1‑149 (MQ=255)
        tttACCTCGATTTGCTGTAGGCGAGCCATCGTAAAGGTGCTAAATGGCGTCAAGCGCGTCAAGAAGGTGTCTATCTGAGGAGTGCGATGTGGCTCAGGAATATACTGGTGAACAACTAAACCACGGCCATAAAGTCATAGACTTTacgg                                                                                                                                           <  2:968632/149‑4 (MQ=255)
        tttAACGAGATTTTCGTTACGCGTGCCATCGTAAACCTGCTATATTTCGGCAAGATTTTCAATAAGCTGTAAATATGATGATTTCGATGTTGATCAGGAATATACTGTTTAAAAAATAAACCACGGAAATAAAGTCATAGACTTTacgg                                                                                                                                           <  2:974317/149‑4 (MQ=255)
        ggtACCGCGCGTGTCGGTTCGCGTGACATCCTAAACCGGCAATATGGCGGCAAGCGCGTCAAGCAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCATAGACTTTacgg                                                                                                                                           <  2:844961/147‑4 (MQ=255)
           aCCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCt                                                                                                                                        >  1:646724/1‑149 (MQ=255)
             cgcgATTGTCGGGACGCGTGCCATCGTAAACCTCCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTgg                                                                                                                                      <  1:791807/149‑1 (MQ=255)
                     tCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTg                                                                                                                              >  2:1375659/1‑149 (MQ=255)
                        gTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTgggg                                                                                                                           >  2:129743/1‑148 (MQ=255)
                        gTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGc                                                                                                                           >  2:490379/1‑149 (MQ=255)
                             cgTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCg                                                                                                                      <  1:1221350/149‑1 (MQ=255)
                                 ccATCGTAAACCTGCTATATGGCGGCAAGCTCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGTAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTTTGACTTTATGCGCTGGGACTACTGGGCTTTCGTCAt                                                                                                                  >  2:1048739/1‑149 (MQ=255)
                                    tCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCAtctc                                                                                                               <  2:454641/149‑1 (MQ=255)
                                    tCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGATGAACAACTAAACCACGGCCGTAAAGTCTATGACATTATGCGCTGGGACTACTGGTCTTTCGGAAtctc                                                                                                               >  1:1525868/1‑149 (MQ=255)
                                     cGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATctct                                                                                                              <  1:1147442/149‑1 (MQ=255)
                                      gTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTg                                                                                                             >  1:941826/1‑149 (MQ=255)
                                      gTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTg                                                                                                             >  1:761474/1‑149 (MQ=255)
                                        aaaCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGt                                                                                                           >  2:739421/1‑149 (MQ=255)
                                        aaaCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGt                                                                                                           >  2:260287/1‑149 (MQ=255)
                                        aaaCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGGCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGt                                                                                                           <  1:1121961/149‑1 (MQ=255)
                                           ccTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctg                                                                                                        >  1:391418/1‑149 (MQ=255)
                                           ccTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctg                                                                                                        >  2:1227887/1‑149 (MQ=255)
                                           ccTGCTATATGGCGGCAAGCGCCTCAAGAAGCTGTCAATCTGAGCAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACAACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctg                                                                                                        >  1:397046/1‑149 (MQ=255)
                                            cTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctgc                                                                                                       <  2:1485878/149‑1 (MQ=255)
                                                    tgcctgcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCGGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATc                                                                                               <  1:792714/144‑1 (MQ=255)
                                                    tGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATc                                                                                               >  2:1428432/1‑149 (MQ=255)
                                                     ggcggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCg                                                                                              >  1:1264487/1‑149 (MQ=255)
                                                       ctgcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgct                                                                                            <  2:397046/147‑1 (MQ=255)
                                                       cggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgct                                                                                            <  1:354739/149‑1 (MQ=255)
                                                       cggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgct                                                                                            <  2:391418/149‑1 (MQ=255)
                                                       cggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgct                                                                                            <  1:1375659/149‑1 (MQ=255)
                                                       cggcAAGCGCGTCAAGAAGATGTCAATCTGAGGAGTGCGATGTGGCTCAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgct                                                                                            <  2:1067845/149‑1 (MQ=255)
                                                        ggcAAGCGCATCTAGAATCTGTCAATCTGATGAGTGCGATGTGGCACAGGAATTTACTGTTGAACAAATAAACCACGGCCGTTAAGTCTTTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctg                                                                                           <  1:193933/149‑1 (MQ=255)
                                                         gcAAGCGCGTCAATTAGCTGTAAATCTGATGAGTGCGATGTTGCACAGGAATATACTGTTGAACAACTAAACCACTGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                          <  2:865795/149‑1 (MQ=255)
                                                         gcAAGCGCGTAACGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAAAATACTGTTGAACCACTAAACCACGGCCGTACAGTCTCTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                          <  1:1003644/149‑1 (MQ=255)
                                                          cAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgct                                                                                         <  1:490379/149‑1 (MQ=255)
                                                           aaGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTa                                                                                        <  1:584270/149‑1 (MQ=255)
                                                              cgcgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCg                                                                                     <  2:1471019/149‑1 (MQ=255)
                                                                 gTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGtta                                                                                  <  1:260287/149‑1 (MQ=255)
                                                                 gTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGCCCGTAAAGTCTATGACTTTATGCGATGGGACTACTGCGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGtta                                                                                  >  1:146940/1‑149 (MQ=255)
                                                                  tCAAGAAGCTGTCAATCTGAGGAGTGCGATGTCGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATCCGCTGCGACTACTGCGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGtaat                                                                                 >  2:1186428/1‑149 (MQ=255)
                                                                    aagaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCGGTTAATCGCTGCTATCGttatta                                                                               >  2:1111805/1‑149 (MQ=255)
                                                                    aagaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTGTCGGCCTCTCTGGTCTGCTGTTAATCGCTGCTATCGttagta                                                                               >  1:1169418/1‑149 (MQ=255)
                                                                     agaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTtattat                                                                              >  1:816931/1‑149 (MQ=255)
                                                                      gaagCTGTCAATCTGAGGATTGCGATGTGGCACAGGATTATACTTTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTTGTCTTCTGTTAATCGCTTCTATCGTTATTATg                                                                             <  1:787253/149‑1 (MQ=255)
                                                                       aagCTTGCAATCTGAGGAGTGCGATGTGGCACAGGAAGATACTGTTGAACAACTACACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACAACTTGGCTATCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATgg                                                                            <  2:896035/149‑1 (MQ=255)
                                                                       aagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAAAATACTGTTGAACAACTAAACCACGGCCTTAAAGTCTATGACTTTATGCGCGGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATgg                                                                            <  2:388437/149‑1 (MQ=255)
                                                                         gCTTTCAATCTGACGAGGGCGTTTTGGCCAAGGAATATACTGTTGAACTAATAAAACACGGCCGTAAAGTATATTATTTTATGCGATGTGACTACTTGGCTTTCGGCATCTATGGTATGCTGTTAATAGCTGCTATCGTTATTATGGGc                                                                          <  2:1027312/149‑1 (MQ=255)
                                                                         gCTGTCAATCTGAGGAGTGCGATGTGGCACAGTAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGAATTTATGCGCTGGGTCTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATggga                                                                          >  1:223372/1‑148 (MQ=255)
                                                                         gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGc                                                                          >  1:938424/1‑149 (MQ=255)
                                                                                aTCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTAAGGGCGTGcgcg                                                                   >  2:321846/1‑149 (MQ=255)
                                                                                    gaggagTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCtt                                                               <  1:946435/149‑1 (MQ=255)
                                                                                       gagTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTGTATGCGCTGGGACTACTGGGCTTTCGGCATCGCTGGTCTGCTGTTAATCGCTGCGATCGTTATGATGGGCGTGCGCGGCTTTac                                                            >  2:1088219/1‑148 (MQ=255)
                                                                                          tGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTACTCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTg                                                         >  1:1419781/1‑149 (MQ=255)
                                                                                           gCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTGATTAGGGGCGTGCGCGGCTTTAACTgg                                                        >  2:908523/1‑149 (MQ=255)
                                                                                             gATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTCGTCTGCTGTTAATCGCTCCTATCGTTATTATGGGCGTGCTCGGCTTTAACTgggg                                                      >  1:215484/1‑149 (MQ=255)
                                                                                              aTGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTggggg                                                     >  2:990327/1‑149 (MQ=255)
                                                                                                  ggCACAGGAATATACTGTTGAACAACTAAACGACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTGTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                 <  2:1422613/149‑1 (MQ=255)
                                                                                                  ggCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                 >  1:1151598/1‑149 (MQ=255)
                                                                                                   gCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCGGGa                                                >  1:1212291/1‑149 (MQ=255)
                                                                                                    cacaGGAATCTACTGTTGAACAACTAAACCCCGGCCGTAAAGGCTATGACTTTATGCGCTGGGAGTACTGGGCTTTCGGCATCTCTTGGCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGAt                                               <  1:1381580/149‑1 (MQ=255)
                                                                                                     acaGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGAtt                                              >  2:48743/1‑149 (MQ=255)
                                                                                                        tgAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGTCTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCa                                           <  2:20142/148‑1 (MQ=255)
                                                                                                          aaTATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCAcc                                         <  2:898405/149‑1 (MQ=255)
                                                                                                               aCTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCTCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                    >  1:1475598/1‑149 (MQ=255)
                                                                                                               aCTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                    <  1:321846/149‑1 (MQ=255)
                                                                                                                 tGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGcgcg                                                                   <  2:15309/116‑1 (MQ=255)
                                                                                                                 tGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGcgcg                                                                   >  1:15309/1‑116 (MQ=255)
                                                                                                                 tGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCACTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTa                                  >  2:1460102/1‑149 (MQ=255)
                                                                                                                 tGTTGAACAACTAAACAACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTa                                  >  2:1354897/1‑149 (MQ=255)
                                                                                                                    tGAACAAATAAACCACGGCCGTAAAGTCTATGACTTTTTGCGCTGGGACTAATTGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATAGTTATTATGGGAGTGAGCGGGTTTAACTGGGTGCTTGATTTCACCGTTGGTACgg                               >  1:1299243/1‑149 (MQ=255)
                                                                                                                     gAACAACTAAACCACGGCCGTAAAGTCTATGACGTTAAGCGCTGGGCCTACTGGGCTTTCGGCATCTCTGGTCTTCTGTTAATCGCTGCTATCGGTATTATCGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCAGCGGTCGTACGGt                              >  1:1454018/1‑149 (MQ=255)
                                                                                                                      aacaacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGtt                             >  1:636786/1‑149 (MQ=255)
                                                                                                                       acaacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCTCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTa                            >  2:119488/1‑149 (MQ=255)
                                                                                                                        caacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAt                           <  2:474022/149‑1 (MQ=255)
                                                                                                                            tAAACCACGGCCGTAAAGTATATGAGTTTATGCGCTGGGAATACTGGGATTTAGGCATATATGGACTGCTGTTAATCCCTGCTATCTTTATTATGGACGTGCGAGGCTTTAACTGGGGGCTGGAATTCACCGGTGGTACGGTTATTGaa                       >  2:1103493/1‑149 (MQ=255)
                                                                                                                             aaaCCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTTTTTTGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGaaa                      >  2:986093/1‑149 (MQ=255)
                                                                                                                                  aCGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACg                 >  2:1263540/1‑149 (MQ=255)
                                                                                                                                    ggCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCt               >  1:1471310/1‑149 (MQ=255)
                                                                                                                                       catAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTGTCGGCATCTATGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGTGTGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGa            >  1:1048691/3‑149 (MQ=255)
                                                                                                                                          aaaGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaaa         >  1:615625/1‑149 (MQ=255)
                                                                                                                                           aaGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaaaa        <  2:615625/149‑1 (MQ=255)
                                                                                                                                             gTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAAcc      <  1:739421/149‑1 (MQ=255)
                                                                                                                                              tCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATTGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGATATTACTCTCGAAAAACCg     >  2:1522328/1‑149 (MQ=255)
                                                                                                                                               cTATGACTTTATGCGCTGTGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTGAATCGCTGCTATCGTTATTATGTGCGTGCGAGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGTTTTTTGAAATTACGCTCGAAAAACCgg    >  2:417118/1‑149 (MQ=255)
                                                                                                                                               cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATTGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCgg    <  1:314354/149‑1 (MQ=255)
                                                                                                                                               cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCgg    <  1:59906/149‑1 (MQ=255)
                                                                                                                                               cTATGACTTTATGCGCGGGGACTACTGGGCTTTCGGCATCTCTGGTCTTCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCgg    <  1:1182477/149‑1 (MQ=255)
                                                                                                                                                 aTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCt  >  1:1497208/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCT  >  NC_000913/429409‑429702

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