Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,907,574 T→C 7.0% D2G (GAC→GGC)  phoU ← negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,907,5740TC 7.0% 2805.7 / 135.3 875D2G (GAC→GGC) phoUnegative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base T (404/409);  new base C (38/23);  total (442/433)
Fisher's exact test for biased strand distribution p-value = 6.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

Only 100 of 892 total aligned reads displayed.

TCGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAG  >  NC_000913/3907425‑3907723
                                                                                                                                                     |                                                                                                                                                     
 cGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTc                                                                                                                                                      <  2:9537753/150‑1 (MQ=255)
 cGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTc                                                                                                                                                      <  2:2586104/150‑1 (MQ=255)
  gCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGccc                                                                                                                                                     <  2:3877566/150‑1 (MQ=255)
  gCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGccc                                                                                                                                                     <  1:4230538/150‑1 (MQ=255)
  gCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTcc                                                                                                                                                     <  2:674799/150‑1 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                    >  1:326195/1‑150 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                    >  1:2694646/1‑150 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCa                                                                                                                                                    <  2:11761526/150‑1 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGATCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTCTCCa                                                                                                                                                    >  1:7375118/1‑150 (MQ=255)
    cAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCAt                                                                                                                                                   <  1:989782/150‑1 (MQ=255)
            tGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACt                                                                                                                                           <  1:9375740/150‑1 (MQ=255)
                  ggTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGaa                                                                                                                                     <  2:4812781/150‑1 (MQ=255)
                   gTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAAt                                                                                                                                    >  1:4432488/1‑150 (MQ=255)
                      aTGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCaa                                                                                                                                 <  1:6645180/150‑1 (MQ=255)
                        gCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAAcc                                                                                                                               <  2:8271260/150‑1 (MQ=255)
                        gCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAAcc                                                                                                                               <  2:5210939/150‑1 (MQ=255)
                          aTCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGt                                                                                                                             <  2:10735062/150‑1 (MQ=255)
                          aTCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGt                                                                                                                             >  1:3482653/1‑150 (MQ=255)
                            cgcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTaa                                                                                                                           <  2:6459015/150‑1 (MQ=255)
                             gcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                          >  1:910542/1‑150 (MQ=255)
                             gcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                          >  2:7453174/1‑150 (MQ=255)
                                  gATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCTTAACGACCg                                                                                                                     >  2:1061432/1‑150 (MQ=255)
                                  gATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCg                                                                                                                     <  2:11254678/150‑1 (MQ=255)
                                    ttGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTCTTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGt                                                                                                                   >  2:6915534/1‑150 (MQ=255)
                                        aTCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATg                                                                                                               <  2:4609715/150‑1 (MQ=255)
                                         tCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGt                                                                                                              >  2:5554567/1‑150 (MQ=255)
                                         tCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGt                                                                                                              >  1:7910639/1‑150 (MQ=255)
                                          cAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACAGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTa                                                                                                             <  2:2140783/150‑1 (MQ=255)
                                            gAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTgcggtcgg                                                                                                           >  1:10000138/1‑143 (MQ=255)
                                               aGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGtctt                                                                                                        >  1:7670327/1‑150 (MQ=255)
                                                   gctgctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgt                                                                                                    >  1:8859181/1‑150 (MQ=255)
                                                     tgcaCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAATATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgttt                                                                                                  <  1:11629656/150‑1 (MQ=255)
                                                        tCCACCATGCCGCCCTTGGTCATCACCTGCGTGCGGATACTTTCCTGTTCGGCTTTGAACTGGCCGGAAATTTGTTTATTAAGATTGAGACTGTCCATACCGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgtt                                                                                               <  1:4786137/150‑1 (MQ=255)
                                                        tCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgtt                                                                                               >  2:3558150/1‑150 (MQ=255)
                                                        tCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgtt                                                                                               >  2:3275428/1‑150 (MQ=255)
                                                        tCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgtt                                                                                               >  2:1981869/1‑150 (MQ=255)
                                                         cccccaGGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGAAACGACCGGTGATGTAGTCTTCTgtttgttt                                                                                              <  1:8459090/147‑1 (MQ=255)
                                                         ccaccaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTGGAACTGGCCGGAAATATGTTTATTAACATTGAGACTGCCCATACCGCCCTCCTGAAACAACCGTAACGACCGGTGATGTAGTCTTCTgtttgttt                                                                                              <  1:4672217/150‑1 (MQ=255)
                                                           accaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttct                                                                                            >  1:7772485/1‑150 (MQ=255)
                                                           accaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttct                                                                                            >  2:11616116/1‑150 (MQ=255)
                                                           accaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttct                                                                                            >  1:1440334/1‑150 (MQ=255)
                                                           accaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttct                                                                                            >  1:1011514/1‑150 (MQ=255)
                                                                 ccgccCATGGTCATCACCTGGGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGTAAATATTTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTg                                                                                      <  1:11281357/150‑1 (MQ=255)
                                                                 ccgccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTg                                                                                      >  1:7160261/1‑150 (MQ=255)
                                                                 ccgacCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGAGAGGTAGTCTTCTGTTTGTTTCTTCGCTg                                                                                      >  2:9175283/1‑150 (MQ=255)
                                                                  cgccCTTGGTCATCACCTGCGTGCGGATACTTTCCATTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTgg                                                                                     <  1:11805793/150‑1 (MQ=255)
                                                                  cgccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTgg                                                                                     >  1:528185/1‑150 (MQ=255)
                                                                     ccATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCtt                                                                                  <  2:5797558/150‑1 (MQ=255)
                                                                            caccaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAAc                                                                           >  2:1768776/4‑150 (MQ=255)
                                                                               caCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGa                                                                        >  2:9392168/1‑150 (MQ=255)
                                                                                   tgcgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGt                                                                    <  1:3687069/150‑1 (MQ=255)
                                                                                     cgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTcc                                                                  <  1:6302431/150‑1 (MQ=255)
                                                                                      gtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCg                                                                 >  2:77620/1‑150 (MQ=255)
                                                                                      gtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCg                                                                 >  1:11057567/1‑150 (MQ=255)
                                                                                        gcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCgtg                                                               >  2:10410206/1‑150 (MQ=255)
                                                                                          ggATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGtt                                                             >  1:8117838/1‑150 (MQ=255)
                                                                                           gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTg                                                            <  2:5983677/150‑1 (MQ=255)
                                                                                           gATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTg                                                            >  2:7216181/1‑150 (MQ=255)
                                                                                              aCTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGCTGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTg                                                         <  1:9518684/150‑1 (MQ=255)
                                                                                                tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTCAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                       <  2:2975347/150‑1 (MQ=255)
                                                                                                tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                       >  1:5553233/1‑150 (MQ=255)
                                                                                                tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                       >  2:2799096/1‑150 (MQ=255)
                                                                                                tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTCGCTTGGTGAACAGATCGTCCGTGATGCTGaa                                                       >  2:12159129/1‑150 (MQ=255)
                                                                                                 ggccAGTTCAACGCCGAACTGGCCGGAAATATGTTTATTAAGATTGAGACAGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAAc                                                      >  2:10263116/3‑150 (MQ=255)
                                                                                                    cAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtca                                                   >  1:12209180/1‑150 (MQ=255)
                                                                                                     aGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaa                                                  >  2:9736090/1‑150 (MQ=255)
                                                                                                     aGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaa                                                  >  1:2082251/1‑150 (MQ=255)
                                                                                                      gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaat                                                 <  1:1805331/150‑1 (MQ=255)
                                                                                                      gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaat                                                 <  1:2405252/150‑1 (MQ=255)
                                                                                                       ttCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaatc                                                >  2:8648407/1‑150 (MQ=255)
                                                                                                        tCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaatca                                               <  1:5805484/150‑1 (MQ=255)
                                                                                                         cGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaatcaa                                              >  2:10919788/1‑150 (MQ=255)
                                                                                                             gTTGAACTGGCCGTAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCg                                          <  2:2314470/150‑1 (MQ=255)
                                                                                                               tGAACTGGCCGTAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGcc                                        <  1:3397360/150‑1 (MQ=255)
                                                                                                               tGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGcc                                        <  1:10508813/150‑1 (MQ=255)
                                                                                                                gAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGccc                                       >  1:2832876/1‑150 (MQ=255)
                                                                                                                 aaCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCa                                      <  1:11968469/150‑1 (MQ=255)
                                                                                                                 aaCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCa                                      <  2:4252715/150‑1 (MQ=255)
                                                                                                                 aaCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCa                                      <  1:6485881/150‑1 (MQ=255)
                                                                                                                  aCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAg                                     <  2:89660/150‑1 (MQ=255)
                                                                                                                  aCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAg                                     <  2:1367242/150‑1 (MQ=255)
                                                                                                                      gCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTAc                                 <  2:12255197/150‑1 (MQ=255)
                                                                                                                      gCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTAc                                 >  2:8868687/1‑150 (MQ=255)
                                                                                                                       cctgAAATATTTTTATTAAGATTGATAATTTCCATAACGCACTCCTGTATCTACCGTAACGTCCGGTGATTTAGTCTTCTGTTTGTTTCTTCTCTTGCTTGTTGAACAGATCGTCAGTGTTGATGAACTCAATCAATTCTACCAGGTACa                                <  2:5052438/147‑1 (MQ=255)
                                                                                                                       ccGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACa                                <  1:9681458/150‑1 (MQ=255)
                                                                                                                            aaTATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGt                                                                               >  1:10675843/1‑98 (MQ=255)
                                                                                                                              tatGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGc                         >  1:4975372/1‑150 (MQ=255)
                                                                                                                                   ttattaATATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTTTTGTTGTCTTCTTTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTAAATCAATTCGCCCAGGTACATAAACGCCgtgt                    <  1:1028255/150‑1 (MQ=255)
                                                                                                                                      ttaAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGt                 >  1:10882151/1‑150 (MQ=255)
                                                                                                                                       taAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTc                <  2:11908034/150‑1 (MQ=255)
                                                                                                                                        aaGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCg               >  2:8506869/1‑150 (MQ=255)
                                                                                                                                          gATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGa             >  1:2983612/1‑150 (MQ=255)
                                                                                                                                          gATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGa             >  2:6182701/1‑150 (MQ=255)
                                                                                                                                           aTTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaa            >  2:1492309/1‑150 (MQ=255)
                                                                                                                                           aTTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaa            >  1:3926014/1‑150 (MQ=255)
                                                                                                                                           aTTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaa            <  1:8610845/150‑1 (MQ=255)
                                                                                                                                              gagaCTGTCCATAACGCACTCCTGAATCAACCTTAACGACCGGTGATGTAGTCTTATGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGTTACATAAACGCCGTGTGGTCGGaacaa         <  2:4049292/150‑1 (MQ=255)
                                                                                                                                               agaCTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaacaac        >  2:8025731/1‑150 (MQ=255)
                                                                                                                                                  cTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAAcgcg     >  2:4455924/1‑150 (MQ=255)
                                                                                                                                                   tGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACgcgc    >  1:6835804/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAG  >  NC_000913/3907425‑3907723

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