Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,230 T→A 5.0% intergenic (‑120/+193) sapA ← / ← ymjA antimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2300TA 5.0% 285.5 / 4.4 100intergenic (‑120/+193)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base T (47/48);  new base A (3/2);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

CCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA  >  NC_000913/1357083‑1357374
                                                                                                                                                   |                                                                                                                                                
ccAATAGAGACGCTAATACATGGCGCATAATATTACTGGATTTAATGTTAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGGTCGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                                                                <  1:287590/150‑1 (MQ=255)
ccAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                                                                >  1:1108024/1‑150 (MQ=255)
       agaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                                                                         >  1:1115358/1‑150 (MQ=255)
          cGATAATACCTGGCGCATAATATTATGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTc                                                                                                                                      >  2:1555409/1‑150 (MQ=255)
           gATAATACCTGGCGCCTACATTTAAGGGATTTTATGTAAAGAGCCTATCTTACTAATATTTAATGACATTGGCCATTACCGCTTGTGTTCATGGGTCTTAATGAGGGCCGTGTGGTTGTTCGTTGAACAGTTTTCGTCGCCATCTGTTCa                                                                                                                                     <  2:940096/150‑1 (MQ=255)
             taataCCTGGCGCATAATATTAAGGTATTTTATGTAAAGAGGCTATATTACTAATATTTAATGACATTTGCCATTACCGTTTGTTTTAAGTTGTCGTAATGAGGTCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                                   <  2:103368/150‑1 (MQ=255)
             taataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                                                                                   <  1:1402827/150‑1 (MQ=255)
               ataCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGTCTATCTTACTAATATTTAATGACATTTGCCATTACCGTATGTGTTCAGCGGTCGTAATGATGGCCGTGTGGTTTGTCGTTGAACAGATGGCAACGTCACCTGTTCAAcga                                                                                                                                 <  1:1493134/150‑3 (MQ=255)
                  ccTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                                                                              >  2:1410993/1‑150 (MQ=255)
                   cTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                                                             >  2:472612/1‑150 (MQ=255)
                      gcgcATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAcc                                                                                                                          <  2:606476/150‑1 (MQ=255)
                      gcgcATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAcc                                                                                                                          <  2:660754/150‑1 (MQ=255)
                       cgcATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                                                                         >  2:1383997/1‑150 (MQ=255)
                        gcATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCg                                                                                                                        <  1:776755/150‑1 (MQ=255)
                          ataataTTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTg                                                                                                                      <  2:1333684/150‑1 (MQ=255)
                             atatTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGCCATTTGCCATTACCGTTTGTGTGCAGGGTTCGTAATTAGGGCCGTGTGTTTGGACGTTGAACAGGTGAGGTCGCCAGATGTTCAACATTCGTCCCCGTGAtt                                                                                                                   >  2:1519309/1‑150 (MQ=255)
                             atatTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAtt                                                                                                                   <  1:947444/150‑1 (MQ=255)
                               atTAAGGGATTTTATGTAAAGAGGCCATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtc                                                                                                                 >  2:722738/1‑150 (MQ=255)
                                  aaGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGTTCGTAATGATGGCAGTGTTGTTGTTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                                                                              >  1:907094/1‑150 (MQ=255)
                                  aaGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                                                                              >  1:632890/1‑150 (MQ=255)
                                     ggATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                                                                           >  1:723026/1‑150 (MQ=255)
                                     ggATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                                                                           <  2:327070/150‑1 (MQ=255)
                                      gATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTAGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataa                                                                                                          >  2:779900/1‑150 (MQ=255)
                                       aTTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                                         <  2:1595945/150‑1 (MQ=255)
                                        tgttATGTAAAGAGGCTCTCTTACTAATATTTAATGACATTTGCCATTTCCGTTAGTGTTCAGGGGGCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataact                                                                                                        <  2:257949/148‑1 (MQ=255)
                                         tttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                                                                       <  2:1452564/150‑1 (MQ=255)
                                         tttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                                                                       >  1:1426983/1‑150 (MQ=255)
                                           tATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataactata                                                                                                     <  2:1516495/150‑1 (MQ=255)
                                           tATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataactata                                                                                                     >  1:1227284/1‑150 (MQ=255)
                                           tATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCAATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataactata                                                                                                     >  2:1263880/1‑150 (MQ=255)
                                             tGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCtataactataac                                                                                                   >  2:331837/1‑150 (MQ=255)
                                                aaaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTca                                                                                                <  2:642374/150‑1 (MQ=255)
                                                aaaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGATTGTCTTCAGGGGTCGTAAAGAGAGCCGTGTGGTTGGTCGTTGAACAGCTGACGTCGCCATCTGTTCAACATTCGAACCCGTGATTTCTCTATAACTAAAACTca                                                                                                >  2:127985/1‑150 (MQ=255)
                                                 aaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTTAGGAGGGCCGTGTGGTGGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                                                                               <  2:1001343/150‑1 (MQ=255)
                                                     ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                           <  1:1396327/150‑1 (MQ=255)
                                                     ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                           <  1:381709/150‑1 (MQ=255)
                                                     ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                           >  2:462263/1‑150 (MQ=255)
                                                      gCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaac                                                                                          <  1:8695/150‑1 (MQ=255)
                                                         aTCTTACTAATATTTAATGACATTTGCCATAACCGTTTGTGTTCAGGGGTCGTAATGAGGGTCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                                                                       <  2:717946/150‑1 (MQ=255)
                                                           cTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                                                                     >  2:838786/1‑150 (MQ=255)
                                                           cTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                                                                     <  1:1383997/150‑1 (MQ=255)
                                                               cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCg                                                                                 >  1:391015/1‑150 (MQ=255)
                                                                  atatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                              <  1:722738/150‑1 (MQ=255)
                                                                  atatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                              <  2:1478125/150‑1 (MQ=255)
                                                                   tatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGGCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagga                                                                             <  1:1296580/150‑1 (MQ=255)
                                                                    atTTAATTAAATTTGCCATTATCGTTTGTGTTCTTGGGTCGTAATGAGGGCCGTGTTGTTGGTCGTTGAACAGGTGACTTGTCCATCTGTTCAACATTCGTAGCCGTGATTTCTCTATATCTATAACTCACAGAACAACTTAGCgaggag                                                                            <  2:1338778/150‑1 (MQ=255)
                                                                      ttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCa                                                                          <  1:1320834/150‑1 (MQ=255)
                                                                           gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                                                                     <  2:1108024/150‑1 (MQ=255)
                                                                            aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                                                                    <  1:1203117/150‑1 (MQ=255)
                                                                            aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                                                                    >  1:103936/1‑150 (MQ=255)
                                                                                 tGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGa                                                               >  2:318800/1‑150 (MQ=255)
                                                                                 tGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGa                                                               >  2:536739/1‑150 (MQ=255)
                                                                                  gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTTGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                                                              >  2:81635/1‑150 (MQ=255)
                                                                                  gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATTACTATAACTCACAGAACAACTTAGCGAGGAGCAGCGCGGTACGAt                                                              >  2:99949/1‑150 (MQ=255)
                                                                                  gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                                                              <  2:907094/150‑1 (MQ=255)
                                                                                  gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                                                              >  1:1516331/1‑150 (MQ=255)
                                                                                  gCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                                                              <  2:632890/150‑1 (MQ=255)
                                                                                    cATTACCGTTTGTGTTCAGGGGGGGTCATGAGGTCCGTGTGGTTGGTCGTTGAACAGGGGACGTCGCCATCTGTTCAACATTCGTACCCTTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGCGGAGCAGAGCGGTACGATCg                                                            <  2:1203748/150‑1 (MQ=255)
                                                                                           gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaaa                                                     >  1:1134129/1‑150 (MQ=255)
                                                                                            tttGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAc                                                    <  1:838786/150‑1 (MQ=255)
                                                                                              tgtgTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAccc                                                  >  1:53598/1‑150 (MQ=255)
                                                                                               gtgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCa                                                 <  1:1267251/150‑1 (MQ=255)
                                                                                               gtgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCa                                                 <  1:331837/150‑1 (MQ=255)
                                                                                                tgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCAt                                                >  1:789617/1‑150 (MQ=255)
                                                                                                tgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCAt                                                >  1:491072/1‑150 (MQ=255)
                                                                                                     aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAg                                           <  1:779900/150‑1 (MQ=255)
                                                                                                       gggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                                         <  1:1278469/150‑1 (MQ=255)
                                                                                                       gggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                                         <  1:472612/150‑1 (MQ=255)
                                                                                                       gggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                                         <  1:487806/150‑1 (MQ=255)
                                                                                                          tCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCAAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGc                                      >  1:124023/1‑150 (MQ=255)
                                                                                                           cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCt                                     >  2:1070977/1‑150 (MQ=255)
                                                                                                             tAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGc                                   >  1:1153619/1‑150 (MQ=255)
                                                                                                             tAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGc                                   >  1:532992/1‑150 (MQ=255)
                                                                                                             tAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGc                                   <  1:1519309/150‑1 (MQ=255)
                                                                                                                 gAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTg                               >  2:1445019/1‑150 (MQ=255)
                                                                                                                     gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCgtg                           >  2:720018/1‑150 (MQ=255)
                                                                                                                     gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCgtg                           <  1:943036/150‑1 (MQ=255)
                                                                                                                      ccGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCgtgt                          >  1:1148247/1‑150 (MQ=255)
                                                                                                                        gtgtGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCa                        >  2:629444/1‑150 (MQ=255)
                                                                                                                         tgtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAt                       <  1:687114/150‑1 (MQ=255)
                                                                                                                          gtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCAtt                      >  1:1091532/1‑150 (MQ=255)
                                                                                                                              ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACATCTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTa                  <  1:941610/150‑1 (MQ=255)
                                                                                                                                ggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCt                                     >  1:1570023/1‑129 (MQ=255)
                                                                                                                                ggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCt                                     <  2:1570023/129‑1 (MQ=255)
                                                                                                                                ggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGc                <  1:464147/150‑1 (MQ=255)
                                                                                                                                ggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGc                >  1:807031/1‑150 (MQ=255)
                                                                                                                                 gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGcc               >  1:827523/1‑150 (MQ=255)
                                                                                                                                  atgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCa              >  1:912770/3‑150 (MQ=255)
                                                                                                                                   tgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAg             >  2:1352696/2‑150 (MQ=255)
                                                                                                                                   cGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGAGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAAGATCGTGAGAAAACCATTAAAGCACGCTGCAGTGAGTGTCAGTGTTAGACAg             >  2:1203445/1‑150 (MQ=255)
                                                                                                                                   cGTTGAACAGGTGAAGTCGCCATCTGTTCAACATTCGTACCAGAGATTTCTCTATAACTATAAATAAAAGAACAAATTAGCGAGTAGCAGAGCGGTACTATCGTGAGAAAAACATTAAAGCAAGATGAAGTGAGTGTCATTGTTAGCCAg             >  2:1203448/1‑150 (MQ=255)
                                                                                                                                     ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAAATTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGAt           >  1:699197/1‑150 (MQ=255)
                                                                                                                                      tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATg          <  1:641175/150‑1 (MQ=255)
                                                                                                                                       gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgc         >  1:1320096/1‑150 (MQ=255)
                                                                                                                                         aCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgcgg       <  1:1263880/150‑1 (MQ=255)
                                                                                                                                          cAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATgcggc      <  2:738441/150‑1 (MQ=255)
                                                                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAAGTCACAGAACAAATTATCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGGGTGTCATTGTTAGACAGATGcgtcgt    >  1:71747/1‑150 (MQ=255)
                                                                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGTGGCGt    >  1:725300/1‑150 (MQ=255)
                                                                                                                                             gTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCAAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTg   <  2:124023/150‑1 (MQ=255)
                                                                                                                                              tGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGa  <  1:462263/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGA  >  NC_000913/1357083‑1357374

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: