Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 429,555 T→A 7.4% Y17* (TAT→TAA) ‡ secF → Sec translocon accessory complex subunit SecF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913429,5550TA 7.4% 289.1 / 9.0 94Y17* (TAT→TAA) ‡secFSec translocon accessory complex subunit SecF
Reads supporting (aligned to +/- strand):  ref base T (45/42);  new base A (3/4);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGC  >  NC_000913/429407‑429701
                                                                                                                                                    |                                                                                                                                                  
gACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATg                                                                                                                                                   <  2:1056962/150‑1 (MQ=255)
gACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATg                                                                                                                                                   <  1:337988/150‑1 (MQ=255)
  cGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCATGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGAc                                                                                                                                                 >  2:723594/1‑150 (MQ=255)
  cGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGAc                                                                                                                                                 >  1:1001164/1‑150 (MQ=255)
  cGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGAc                                                                                                                                                 >  2:1174284/1‑150 (MQ=255)
      gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTa                                                                                                                                             >  1:644966/1‑150 (MQ=255)
      gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTa                                                                                                                                             >  1:744602/1‑150 (MQ=255)
      gATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTa                                                                                                                                             <  2:453655/150‑1 (MQ=255)
        tGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATg                                                                                                                                           >  2:1356523/1‑150 (MQ=255)
         gTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATgc                                                                                                                                          <  2:644966/150‑1 (MQ=255)
         gTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCATAGACTTTacgg                                                                                                                                          <  1:1173942/150‑4 (MQ=255)
          tttACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCATAGACTTTacggc                                                                                                                                         <  1:174638/150‑5 (MQ=255)
            tACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCt                                                                                                                                       >  2:1320253/1‑150 (MQ=255)
            tACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCATAGACTTTacggccg                                                                                                                                       <  2:964376/150‑7 (MQ=255)
             aCCGCGATTGGCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGTCCGTAAAGTCATAGACTTTacggccgt                                                                                                                                      <  1:950047/150‑8 (MQ=255)
              ccGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTgg                                                                                                                                     <  1:1380112/150‑1 (MQ=255)
               cgcgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTggg                                                                                                                                    <  1:974863/150‑1 (MQ=255)
               cgcgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTggg                                                                                                                                    <  1:819280/150‑1 (MQ=255)
               cgcgATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTggg                                                                                                                                    <  1:1231813/150‑1 (MQ=255)
                  gATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGact                                                                                                                                 <  2:776343/150‑1 (MQ=255)
                     tGTCGGTACGAGTGCCATCGTAATCCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGGCAATCTGAGGAGTGCGATGTGGAACAGGAATTTACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGactact                                                                                                                              <  2:1337316/150‑1 (MQ=255)
                      gTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTg                                                                                                                             <  2:908659/150‑1 (MQ=255)
                        cGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTggg                                                                                                                           >  1:1287553/1‑150 (MQ=255)
                        cGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTggg                                                                                                                           <  2:778577/150‑1 (MQ=255)
                         ggTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGc                                                                                                                          >  1:995626/1‑150 (MQ=255)
                              gcgTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCg                                                                                                                     >  2:1401968/1‑150 (MQ=255)
                                    cATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCAtct                                                                                                               <  1:1134300/150‑1 (MQ=255)
                                    cATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCAtct                                                                                                               >  1:1618741/1‑150 (MQ=255)
                                     aTCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCAtctc                                                                                                              <  1:555925/150‑1 (MQ=255)
                                         tAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGt                                                                                                          <  2:384706/150‑1 (MQ=255)
                                             ccTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctgc                                                                                                      >  1:233095/1‑150 (MQ=255)
                                              cTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctgct                                                                                                     >  1:1462650/1‑150 (MQ=255)
                                              cTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctgct                                                                                                     >  1:82194/1‑150 (MQ=255)
                                               tGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTctgctg                                                                                                    <  2:1001164/150‑1 (MQ=255)
                                                 cTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                  <  1:571806/150‑1 (MQ=255)
                                                   agatGGCGGGAAGCGCGGCTATAAGCTGTAGCTGTGAGGAGGGGGATGTGGCTCAGGTATTTACTGTTCAACAACTTAACCAGGGCCGTAAACTCTATGCGTTTATGCGCGGCGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTaa                                                                                                <  2:494526/148‑1 (MQ=255)
                                                    tatGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAAt                                                                                               <  1:1174518/150‑1 (MQ=255)
                                                    tatGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAAt                                                                                               >  2:290300/1‑150 (MQ=255)
                                                      tGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                  <  1:1600519/145‑1 (MQ=255)
                                                      tGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGtt                                                                                                  >  2:1600519/1‑145 (MQ=255)
                                                          ggcATGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCgctgc                                                                                         <  1:723594/150‑1 (MQ=255)
                                                          ggcAAGCGCGGCAAGAAGCTGTCAATCTGAGGAGGGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGGCTGCTGTTAATCgctgc                                                                                         <  2:1397606/150‑1 (MQ=255)
                                                                cgcgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGt                                                                                   <  1:870631/150‑1 (MQ=255)
                                                                  cgTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGtta                                                                                 >  1:1565000/1‑150 (MQ=255)
                                                                     cAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGttatta                                                                              <  2:333799/150‑1 (MQ=255)
                                                                     cAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGttatta                                                                              <  2:395577/150‑1 (MQ=255)
                                                                       agaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATg                                                                            >  2:274643/1‑150 (MQ=255)
                                                                       agaagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATg                                                                            <  2:1462650/150‑1 (MQ=255)
                                                                         aagCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATggg                                                                          <  2:33515/150‑1 (MQ=255)
                                                                           gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCg                                                                        >  1:1148044/1‑150 (MQ=255)
                                                                           gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCg                                                                        <  2:269008/150‑1 (MQ=255)
                                                                           gCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACAACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTATGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCg                                                                        >  2:86765/1‑150 (MQ=255)
                                                                             tGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTg                                                                      >  2:131866/1‑150 (MQ=255)
                                                                              gTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTgc                                                                     <  2:1287553/150‑1 (MQ=255)
                                                                                 aaTCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGcgcg                                                                  <  1:1356523/150‑1 (MQ=255)
                                                                                   tCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGc                                                                <  1:290300/150‑1 (MQ=255)
                                                                                     ttaggacTGCGATTTGGCACAGGAATATACTGTTGAAAAAAAACACAACTGCCGTAATGTATATTAATTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCtt                                                              >  2:1508263/3‑150 (MQ=255)
                                                                                      gaggagTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCttt                                                             <  2:970935/150‑1 (MQ=255)
                                                                                       aggaggGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCGGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGGGCGGCTTTa                                                            <  2:944969/150‑1 (MQ=255)
                                                                                         gagTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTa                                                            >  2:528211/1‑148 (MQ=255)
                                                                                         gagTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATCGGCGTGCGCGGCTTTa                                                            <  1:528211/148‑1 (MQ=255)
                                                                                            tGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTgg                                                       >  2:725420/1‑150 (MQ=255)
                                                                                             gCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTggg                                                      >  1:209675/1‑150 (MQ=255)
                                                                                                   tgGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                >  2:306794/1‑150 (MQ=255)
                                                                                                   tgGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTgg                                                >  2:642969/1‑150 (MQ=255)
                                                                                                    ggCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGa                                               >  2:1563152/1‑150 (MQ=255)
                                                                                                    ggCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGa                                               >  2:975050/1‑150 (MQ=255)
                                                                                                      cacaGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGAtt                                             >  2:723650/1‑150 (MQ=255)
                                                                                                      aacaGGAATATACTGTTGAACAACTATACCACGGCCGTAAATTCTTTGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGTGCGTGCGCGGCTTTAACTGGGGGCTGTAtt                                             <  1:1033743/149‑1 (MQ=255)
                                                                                                          ggAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCAc                                         >  2:1571171/1‑150 (MQ=255)
                                                                                                                taCTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                   >  2:1343819/1‑150 (MQ=255)
                                                                                                                taCTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCggtgg                                   <  2:233095/150‑1 (MQ=255)
                                                                                                                      tGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGt                             >  1:490032/1‑150 (MQ=255)
                                                                                                                      tGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGt                             <  1:1563152/150‑1 (MQ=255)
                                                                                                                      tGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGt                             >  1:743610/1‑150 (MQ=255)
                                                                                                                       gAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGtt                            >  1:958727/1‑150 (MQ=255)
                                                                                                                       gAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGtt                            >  2:145389/1‑150 (MQ=255)
                                                                                                                       gAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGtt                            >  2:811650/1‑150 (MQ=255)
                                                                                                                          caacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAtt                         <  1:1401968/150‑1 (MQ=255)
                                                                                                                          caacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTAtt                         >  2:1035488/1‑150 (MQ=255)
                                                                                                                           aacTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTg                        <  2:958727/150‑1 (MQ=255)
                                                                                                                            acTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGa                       >  2:1244832/1‑150 (MQ=255)
                                                                                                                            acTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGa                       >  1:601191/1‑150 (MQ=255)
                                                                                                                                  ccACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTAc                 >  1:1610917/1‑150 (MQ=255)
                                                                                                                                   cACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGTGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACg                <  2:995626/150‑1 (MQ=255)
                                                                                                                                    aCGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACgg               >  1:323960/1‑149 (MQ=255)
                                                                                                                                    aCGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGc               >  1:96155/1‑150 (MQ=255)
                                                                                                                                        gcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGa           >  2:317150/4‑150 (MQ=255)
                                                                                                                                        gcatAAAGTCATAGACTTTATGCGCTGGGACTACTGGGATTTCGGCATCTCTGGTCTGCTGTTAATAGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGa           >  2:285328/4‑150 (MQ=255)
                                                                                                                                         catAAAGTCATAGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaa          >  1:24601/3‑150 (MQ=255)
                                                                                                                                            aaaGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaaaa       >  2:412290/1‑150 (MQ=255)
                                                                                                                                            aaaGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGaaaaa       >  2:879833/1‑150 (MQ=255)
                                                                                                                                                 cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGc  <  1:1244832/150‑1 (MQ=255)
                                                                                                                                                 cTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGc  <  1:1174284/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACGTCGATGTTTACCGCGATTGTCGGTACGCGTGCCATCGTAAACCTGCTATATGGCGGCAAGCGCGTCAAGAAGCTGTCAATCTGAGGAGTGCGATGTGGCACAGGAATATACTGTTGAACAACTAAACCACGGCCGTAAAGTCTATGACTTTATGCGCTGGGACTACTGGGCTTTCGGCATCTCTGGTCTGCTGTTAATCGCTGCTATCGTTATTATGGGCGTGCGCGGCTTTAACTGGGGGCTGGATTTCACCGGTGGTACGGTTATTGAAATTACGCTCGAAAAACCGGC  >  NC_000913/429407‑429701

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