Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 429,554 A→T 7.2% 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 7.2% 384.8 / 12.9 137Y17F (TAT→TTT) ‡secFSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base A (60/66);  new base T (6/4);  total (66/71)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-01

Only 100 of 137 total aligned reads displayed.

CGACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCT  >  NC_000913/429406‑429702
                                                                                                                                                    |                                                                                                                                                    
cGACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTAt                                                                                                                                                     <  2:349621/150‑1 (MQ=255)
 gACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCCCGGCCGTAAAGTCTATg                                                                                                                                                    >  2:632285/1‑150 (MQ=255)
 gACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATg                                                                                                                                                    >  2:1481037/1‑150 (MQ=255)
  aCGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGa                                                                                                                                                   <  2:352359/150‑1 (MQ=255)
   cGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCAGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATCTACTGTTGCACAACTAAACCACGGCCGTAAAGTCTATGAc                                                                                                                                                  >  1:350794/1‑150 (MQ=255)
     tCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACtt                                                                                                                                                <  2:1347696/150‑1 (MQ=255)
     tCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACtt                                                                                                                                                >  2:484830/1‑150 (MQ=255)
      cGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACttt                                                                                                                                               >  1:1086797/1‑150 (MQ=255)
       gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTa                                                                                                                                              <  1:486320/150‑1 (MQ=255)
        aTGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTTCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTAt                                                                                                                                             <  2:1711515/150‑1 (MQ=255)
         tGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATg                                                                                                                                            <  2:386121/150‑1 (MQ=255)
          gTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCATAGACTTTacgg                                                                                                                                           <  1:1384153/150‑4 (MQ=255)
           tttACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATgcg                                                                                                                                          <  1:341236/150‑1 (MQ=255)
            ttACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATgcgc                                                                                                                                         <  1:1356659/150‑1 (MQ=255)
            ttACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATgcgc                                                                                                                                         <  2:796448/150‑1 (MQ=255)
              aCCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCATAAAGTCATAGACTTTacggccgt                                                                                                                                       <  1:1690670/150‑8 (MQ=255)
                cgcgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTggg                                                                                                                                     >  2:1630705/1‑150 (MQ=255)
                  cgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGac                                                                                                                                   <  1:508055/150‑1 (MQ=255)
                        tCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTgg                                                                                                                             <  2:805734/150‑1 (MQ=255)
                         cGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTggg                                                                                                                            >  1:702155/1‑150 (MQ=255)
                          ggTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGc                                                                                                                           <  2:1277497/150‑1 (MQ=255)
                          ggTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGc                                                                                                                           >  1:859547/1‑150 (MQ=255)
                          ggTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCGGGCAATCTGAGGAGGGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGc                                                                                                                           <  2:1284281/150‑1 (MQ=255)
                              cgcgTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTc                                                                                                                       >  1:909116/1‑150 (MQ=255)
                               gcgTGCCATCGTAAACCTGCTATATGGCAGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATTCTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCg                                                                                                                      >  1:235468/1‑150 (MQ=255)
                                cgTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTTTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCgg                                                                                                                     <  1:1429540/150‑1 (MQ=255)
                                  tGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACCGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCa                                                                                                                   <  1:1449498/150‑1 (MQ=255)
                                             aCCTGCTATAAGGCGGCAAGCGCGTCAAGAAGCTGGCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctg                                                                                                        <  1:119120/150‑1 (MQ=255)
                                                tGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCAAGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCAGCTCTGGTctgctg                                                                                                     >  2:1073524/1‑150 (MQ=255)
                                                 gCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGt                                                                                                    >  2:1009554/1‑150 (MQ=255)
                                                 gCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGt                                                                                                    <  2:489147/150‑1 (MQ=255)
                                                  cTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                   <  2:169369/150‑1 (MQ=255)
                                                  cTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                   <  1:676277/150‑1 (MQ=255)
                                                   tataTGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTa                                                                                                  >  2:1674243/1‑150 (MQ=255)
                                                          cggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctg                                                                                           <  1:1309317/150‑1 (MQ=255)
                                                          cggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctg                                                                                           >  1:649915/1‑150 (MQ=255)
                                                           ggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                          <  1:828528/150‑1 (MQ=255)
                                                           ggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                          <  2:205157/150‑1 (MQ=255)
                                                           ggcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACCAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                          <  2:194021/150‑1 (MQ=255)
                                                            gcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgct                                                                                         >  1:43590/1‑150 (MQ=255)
                                                            gcAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgct                                                                                         >  2:1266057/1‑150 (MQ=255)
                                                            gcAAGCGAGTCAAGAAGATGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgct                                                                                         <  2:946955/150‑1 (MQ=255)
                                                                 cgcgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGt                                                                                    <  1:732608/150‑1 (MQ=255)
                                                                  gcgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGtt                                                                                   <  1:1204101/150‑1 (MQ=255)
                                                                   cgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGtga                                                                                  >  1:42653/1‑148 (MQ=255)
                                                                   cgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGATGTTAATCGCTGCTATCGtta                                                                                  >  1:229192/1‑150 (MQ=255)
                                                                      cAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCTCAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTGCGGCATCTCTGGGCTGCTGTTAAGCGCTGCTATCGttattt                                                                               >  1:99073/1‑149 (MQ=255)
                                                                      cAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGttatta                                                                               >  1:1121685/1‑150 (MQ=255)
                                                                      cAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGttatta                                                                               <  1:964653/150‑1 (MQ=255)
                                                                      cAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGGTAATCGCTGCTATCGttatta                                                                               >  2:1397771/1‑150 (MQ=255)
                                                                       aagaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTtattat                                                                              >  2:1705085/1‑150 (MQ=255)
                                                                        agaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATg                                                                             >  2:363718/1‑150 (MQ=255)
                                                                            gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTTGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCg                                                                         >  2:1418250/1‑150 (MQ=255)
                                                                            gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCg                                                                         >  1:1134969/1‑150 (MQ=255)
                                                                             cTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGt                                                                        >  2:1108919/1‑150 (MQ=255)
                                                                               gTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTgc                                                                      >  1:466988/1‑150 (MQ=255)
                                                                               gTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTgc                                                                      <  1:1481037/150‑1 (MQ=255)
                                                                                 cAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTgcgc                                                                    >  1:551488/1‑150 (MQ=255)
                                                                                      tGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCtt                                                               <  1:420418/150‑1 (MQ=255)
                                                                                       taggatggCGAGGTGGCACAGGAATATAATGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTGCGGCATCTCTGGTCTGCGGTTATTCGCTGCTATCGTTATTATGGGCGTGCGCGGCttt                                                              <  2:1251113/143‑1 (MQ=255)
                                                                                            gTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTg                                                         >  2:199598/1‑150 (MQ=255)
                                                                                              gCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTggg                                                       <  2:1085050/150‑1 (MQ=255)
                                                                                              gCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGGCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTggg                                                       <  1:940571/150‑1 (MQ=255)
                                                                                                  tgtgGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCt                                                   <  2:551488/150‑1 (MQ=255)
                                                                                                   ggggCACAGGAATATACTGTTGAACAACTACACCACGGCCGTAAAGTCACTGACTTTATGCGCTGGGACTAGTGTGCTTTCGGCTTCTCTGGTCTGCTGTTTCTCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTg                                                  <  1:249140/148‑1 (MQ=255)
                                                                                                    tgGCACAGGAATATACTGTTGAACAACTAAACCTCGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                 <  2:1455253/150‑1 (MQ=255)
                                                                                                    tgGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                 >  1:652622/1‑150 (MQ=255)
                                                                                                     ggCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGa                                                >  2:1115859/1‑150 (MQ=255)
                                                                                                     ggCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGa                                                <  1:1115859/150‑1 (MQ=255)
                                                                                                         caGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTc                                            <  2:1134969/150‑1 (MQ=255)
                                                                                                         caGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTc                                            <  1:1566931/150‑1 (MQ=255)
                                                                                                           gtaaTATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGTCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATgggg                                                                          <  1:822512/116‑2 (MQ=255)
                                                                                                           ggAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGc                                                                          >  2:822512/1‑118 (MQ=255)
                                                                                                            gAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCAcc                                         <  1:544126/150‑1 (MQ=255)
                                                                                                                 taCTGTTGAACAACTAAACCCCGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTAGTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                    >  2:821292/1‑150 (MQ=255)
                                                                                                                 taCTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                    <  2:51601/150‑1 (MQ=255)
                                                                                                                           caacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAtt                          <  1:1231813/150‑1 (MQ=255)
                                                                                                                            aacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTg                         >  1:960658/1‑150 (MQ=255)
                                                                                                                            aacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTg                         <  2:43590/150‑1 (MQ=255)
                                                                                                                             acTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGa                        <  1:574752/150‑1 (MQ=255)
                                                                                                                              cTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGaa                       <  1:1009554/150‑1 (MQ=255)
                                                                                                                                aaaCCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGTTGTTAGGGGTATTGAAAt                     >  2:306453/1‑150 (MQ=255)
                                                                                                                                 aaCCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAAtt                    <  1:199598/150‑1 (MQ=255)
                                                                                                                                 aaCCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCAGCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAAtt                    <  2:702155/150‑1 (MQ=255)
                                                                                                                                   ccACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTAc                  <  1:1705085/150‑1 (MQ=255)
                                                                                                                                    cACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACg                 <  2:42055/150‑1 (MQ=255)
                                                                                                                                     aCGGGCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGGCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGc                <  1:1104125/150‑1 (MQ=255)
                                                                                                                                      cGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCt               >  1:207266/1‑150 (MQ=255)
                                                                                                                                       gcgcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGAAGGATTTCACCGGTGGTACGGTTATTGAAATTACtctc              >  1:278982/6‑150 (MQ=255)
                                                                                                                                        cgcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCTTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCg             >  2:12927/5‑150 (MQ=255)
                                                                                                                                        cgcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCg             >  2:661633/5‑150 (MQ=255)
                                                                                                                                         ccGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGATGGTTTACACCGGTGGTACGGTATTTAAAATTCCGcacaa            >  2:1216073/1‑146 (MQ=255)
                                                                                                                                          catAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaa           >  2:1188775/3‑150 (MQ=255)
                                                                                                                                          catAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTCGGGGCTGGATTTCACCGGTGGTACCGTTATTGAAATTACGCTCGaa           >  2:776757/3‑150 (MQ=255)
                                                                                                                                          cGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTTTTGAAATTACGCTCGaa           >  1:1512752/1‑150 (MQ=255)
                                                                                                                                            tccaggcTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTTATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGTGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaaa         <  2:601753/144‑1 (MQ=255)
                                                                                                                                            tAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaaa         <  2:909116/150‑1 (MQ=255)
                                                                                                                                                  cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGc   >  2:1182577/1‑150 (MQ=255)
                                                                                                                                                  cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGc   >  2:311847/1‑150 (MQ=255)
                                                                                                                                                   tATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCt  <  2:743052/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGCT  >  NC_000913/429406‑429702

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