Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,630,857 C→A 5.9% intergenic (‑255/+111) yhiI ← / ← yhiJ putative membrane fusion protein YhiI/DUF4049 domain‑containing protein YhiJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,630,8570CA 5.9% 434.2 / 10.5 129intergenic (‑255/+111)yhiI/yhiJputative membrane fusion protein YhiI/DUF4049 domain‑containing protein YhiJ
Reads supporting (aligned to +/- strand):  ref base C (58/63);  new base A (2/6);  total (60/69)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 140 total aligned reads displayed.

CGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAT  >  NC_000913/3630708‑3631006
                                                                                                                                                     |                                                                                                                                                     
cGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCAATCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTatccatcc                                                                                                                                                       <  2:96880/150‑1 (MQ=255)
     aGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTg                                                                                                                                                  >  1:93520/1‑150 (MQ=255)
      gcgcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTgg                                                                                                                                                 >  1:669665/1‑150 (MQ=255)
        gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGGAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCATGTTGGAt                                                                                                                                               >  1:737761/1‑150 (MQ=255)
        gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGAt                                                                                                                                               >  1:1423340/1‑150 (MQ=255)
        gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGAt                                                                                                                                               <  2:871401/150‑1 (MQ=255)
         cgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATa                                                                                                                                              >  1:356225/1‑150 (MQ=255)
         cgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATa                                                                                                                                              >  1:283083/1‑150 (MQ=255)
           cAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGc                                                                                                                                            >  1:1531646/1‑150 (MQ=255)
             gagCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGt                                                                                                                                          <  2:606680/150‑1 (MQ=255)
             gagCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGt                                                                                                                                          <  2:611680/150‑1 (MQ=255)
               gCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGttt                                                                                                                                        >  1:938790/1‑150 (MQ=255)
               gCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttt                                                                                                                                        <  1:619338/150‑1 (MQ=255)
                  ttGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCAc                                                                                                                                     >  1:1419494/1‑150 (MQ=255)
                   tGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACg                                                                                                                                    <  2:1136902/150‑1 (MQ=255)
                     aGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGcc                                                                                                                                  >  2:253013/1‑150 (MQ=255)
                       cgcCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTg                                                                                                                                <  2:1428672/150‑1 (MQ=255)
                       cgcCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTg                                                                                                                                <  1:553949/150‑1 (MQ=255)
                       cgcCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTg                                                                                                                                <  2:849469/150‑1 (MQ=255)
                        gcCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGt                                                                                                                               >  1:1061595/1‑150 (MQ=255)
                         ccAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGtt                                                                                                                              >  1:1283567/1‑150 (MQ=255)
                          cAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTc                                                                                                                             <  2:110175/150‑1 (MQ=255)
                           aaCGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCGGGTCa                                                                                                                            >  2:1394882/1‑150 (MQ=255)
                            aCGGCGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAATCTCGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAg                                                                                                                           <  1:1137605/150‑1 (MQ=255)
                               gggAAGGCTTAGTGAAAACGGGTAGTAATATTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTg                                                                                                                        >  1:342018/1‑150 (MQ=255)
                                      cTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAAc                                                                                                                 <  1:99972/150‑1 (MQ=255)
                                      cTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAAc                                                                                                                 >  1:947592/1‑150 (MQ=255)
                                        tAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACag                                                                                                               >  2:1001542/1‑150 (MQ=255)
                                        tAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACag                                                                                                               >  2:932259/1‑150 (MQ=255)
                                          gTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACagag                                                                                                             <  1:1337780/150‑1 (MQ=255)
                                                  gggTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtt                                                                                                     <  2:1041365/150‑1 (MQ=255)
                                                   ggTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCttt                                                                                                    <  1:388969/150‑1 (MQ=255)
                                                    gtagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtttt                                                                                                   <  2:1062016/150‑1 (MQ=255)
                                                    gtagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtttt                                                                                                   <  2:162058/150‑1 (MQ=255)
                                                     tagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACTGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGATCATTGCttttt                                                                                                  <  2:96907/150‑1 (MQ=255)
                                                     tagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCttttt                                                                                                  <  2:1507382/150‑1 (MQ=255)
                                                      agtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTa                                                                                                 <  1:845771/150‑1 (MQ=255)
                                                              tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                         >  1:781125/1‑150 (MQ=255)
                                                              tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                         <  2:938790/150‑1 (MQ=255)
                                                              tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                         <  2:782745/150‑1 (MQ=255)
                                                               gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCGTTCCTTTTTACGTTTTATc                                                                                        <  2:1309784/150‑1 (MQ=255)
                                                               gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATc                                                                                        >  1:1345865/1‑150 (MQ=255)
                                                                tcatcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCa                                                                                       <  2:1419494/150‑1 (MQ=255)
                                                                tcatcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCa                                                                                       <  2:734890/150‑1 (MQ=255)
                                                                    catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                   >  2:584142/1‑150 (MQ=255)
                                                                    catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                   <  1:185533/150‑1 (MQ=255)
                                                                      tcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtc                                                                                 >  1:740637/1‑150 (MQ=255)
                                                                       caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtct                                                                                >  2:732204/1‑150 (MQ=255)
                                                                        aGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtctc                                                                               >  1:531181/1‑150 (MQ=255)
                                                                        aGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtctc                                                                               <  2:1423340/150‑1 (MQ=255)
                                                                              ccTGGCCGTAAATGAAACACTCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAg                                                                         <  2:305454/150‑1 (MQ=255)
                                                                              ccTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAg                                                                         <  2:1009444/150‑1 (MQ=255)
                                                                              ccTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAg                                                                         >  1:896116/1‑150 (MQ=255)
                                                                                 ggCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATa                                                                      >  2:340656/1‑150 (MQ=255)
                                                                                          tGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTa                                                             <  1:932259/150‑1 (MQ=255)
                                                                                           gAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTaa                                                            <  2:836046/150‑1 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                          >  1:262731/1‑150 (MQ=255)
                                                                                               cacGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATAtt                                                        >  1:792359/1‑150 (MQ=255)
                                                                                                 cGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGTAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                                      >  2:926221/1‑150 (MQ=255)
                                                                                                 cGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                                      <  1:967019/150‑1 (MQ=255)
                                                                                                   cATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCAtg                                                    <  2:1061595/150‑1 (MQ=255)
                                                                                                     tACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGt                                                                    <  2:629510/132‑1 (MQ=255)
                                                                                                        tCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTc                                               <  2:669665/150‑1 (MQ=255)
                                                                                                         cATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCt                                              >  1:1105555/1‑150 (MQ=255)
                                                                                                              gtgtGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGt                                         <  2:152768/150‑1 (MQ=255)
                                                                                                                 tgATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCa                                      >  1:354742/1‑150 (MQ=255)
                                                                                                                  gATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCag                                     >  2:1042303/1‑150 (MQ=255)
                                                                                                                    tCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCagag                                   >  1:670525/1‑150 (MQ=255)
                                                                                                                    tCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCagag                                   <  2:155050/150‑1 (MQ=255)
                                                                                                                      tataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTAGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAAACAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCt                                 >  2:1399961/1‑150 (MQ=255)
                                                                                                                       ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                <  2:356225/150‑1 (MQ=255)
                                                                                                                       ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                >  2:479186/1‑150 (MQ=255)
                                                                                                                       ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                <  2:1236621/150‑1 (MQ=255)
                                                                                                                       ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                >  2:583387/1‑150 (MQ=255)
                                                                                                                         aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTa                              <  1:185050/150‑1 (MQ=255)
                                                                                                                          aaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTaa                             >  1:47109/1‑150 (MQ=255)
                                                                                                                          aaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTaa                             >  2:1210478/1‑150 (MQ=255)
                                                                                                                             cTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACAt                          <  1:1538672/150‑1 (MQ=255)
                                                                                                                              tGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACAtt                         >  1:538828/1‑150 (MQ=255)
                                                                                                                               ggCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtg                        >  1:1012741/1‑150 (MQ=255)
                                                                                                                                 cTGACAAAACGGTATCCATCCTGTTGGATATCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtgtg                      <  2:913850/150‑1 (MQ=255)
                                                                                                                                 cTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtgtg                      <  2:403056/150‑1 (MQ=255)
                                                                                                                                   gACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgcg                    >  1:1047365/1‑150 (MQ=255)
                                                                                                                                     cAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGcc                  >  2:189439/1‑150 (MQ=255)
                                                                                                                                       aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTc                                               >  1:1188978/1‑119 (MQ=255)
                                                                                                                                       aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCtt                <  2:670525/150‑1 (MQ=255)
                                                                                                                                       aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCtt                >  2:167761/1‑150 (MQ=255)
                                                                                                                                       aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCtt                <  1:1001542/150‑1 (MQ=255)
                                                                                                                                        aaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCtta               >  1:661805/1‑149 (MQ=255)
                                                                                                                                             tATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTCCAGAGCTGTAACATTGTGCGCCTTCTAATc          <  2:379423/150‑1 (MQ=255)
                                                                                                                                              atccatccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCGCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCa         <  1:1435049/150‑1 (MQ=255)
                                                                                                                                                 catccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAg      >  1:603603/1‑150 (MQ=255)
                                                                                                                                                 catccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAg      <  2:47109/150‑1 (MQ=255)
                                                                                                                                                 catccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAg      <  1:748501/150‑1 (MQ=255)
                                                                                                                                                   tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGc    <  2:1164438/150‑1 (MQ=255)
                                                                                                                                                   tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGc    <  1:189439/150‑1 (MQ=255)
                                                                                                                                                    ccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCa   >  2:917624/1‑150 (MQ=255)
                                                                                                                                                    ccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCa   >  2:746545/1‑150 (MQ=255)
                                                                                                                                                     cTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAt  <  2:781125/150‑1 (MQ=255)
                                                                                                                                                     cTGTTGGATAGCGTTTCACGCCTATTCAGTTGCAGGAACAGAGCTTTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAt  >  2:524515/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAT  >  NC_000913/3630708‑3631006

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: