Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,907,574 T→C 5.2% 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 5.2% 1083.2 / 28.6 326D2G (GAC→GGC) phoUnegative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base T (147/161);  new base C (6/11);  total (154/172)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

TCGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAG  >  NC_000913/3907425‑3907723
                                                                                                                                                     |                                                                                                                                                     
 cGCCAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTc                                                                                                                                                      <  1:1090118/150‑1 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                    >  1:2606743/1‑150 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                    >  1:3679355/1‑150 (MQ=255)
   ccAGATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCa                                                                                                                                                    >  1:3730912/1‑150 (MQ=255)
      gATCGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATaa                                                                                                                                                 >  1:3593795/1‑150 (MQ=255)
         cGCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGc                                                                                                                                              <  2:4021224/150‑1 (MQ=255)
          gCTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCa                                                                                                                                             <  1:2805796/150‑1 (MQ=255)
           cTGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCAc                                                                                                                                            <  1:809235/150‑1 (MQ=255)
            tGTCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGACCTATGTTTATTAAGATTGAGACTGTCCATATCGCACt                                                                                                                                           >  1:4092595/1‑150 (MQ=255)
              tCCTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGTCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTcc                                                                                                                                         >  1:1271872/1‑150 (MQ=255)
                cTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTg                                                                                                                                       <  1:3392038/150‑1 (MQ=255)
                cTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTg                                                                                                                                       >  2:3124181/1‑150 (MQ=255)
                cTGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTg                                                                                                                                       >  2:3489806/1‑150 (MQ=255)
                 tGGTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGa                                                                                                                                      >  1:2488457/1‑150 (MQ=255)
                  ggTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGaa                                                                                                                                     >  1:4025022/1‑150 (MQ=255)
                   gTTATGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAAt                                                                                                                                    >  2:3416414/1‑150 (MQ=255)
                       tGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAAc                                                                                                                                <  1:781349/150‑1 (MQ=255)
                       tGCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAAc                                                                                                                                <  2:50388/150‑1 (MQ=255)
                        gCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGTGACTGTCCATAACGCACTCCTGAATCAAcc                                                                                                                               >  2:1841211/1‑150 (MQ=255)
                        gCATCGCGGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAAcc                                                                                                                               >  2:2511421/1‑150 (MQ=255)
                             gcgGTGATTGCATCAGAAAGCTGCTGCTCCACTATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                          >  2:2878414/1‑150 (MQ=255)
                             gcgGTGATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAAc                                                                                                                          <  2:3450073/150‑1 (MQ=255)
                                  gATTGCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCg                                                                                                                     <  2:1347931/150‑1 (MQ=255)
                                      gCATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGa                                                                                                                 <  2:3349025/150‑1 (MQ=255)
                                       cATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGAt                                                                                                                >  2:3839860/1‑150 (MQ=255)
                                       cATCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGAt                                                                                                                >  2:194566/1‑150 (MQ=255)
                                        aTCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATg                                                                                                               <  2:2030530/150‑1 (MQ=255)
                                         tCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGt                                                                                                              >  2:4033340/1‑150 (MQ=255)
                                         tCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGt                                                                                                              >  1:3962428/1‑150 (MQ=255)
                                         tCAGAAAGCTGCTGCTCCACCATGCCGCCCATGGACATCACCTGCGAGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTGTATTAAGATTGGGACTGTCCATAACGCACACCTGAATAAACCGTAACGACCGGTGATGt                                                                                                              >  1:3566477/1‑150 (MQ=255)
                                                   gctgctCCACCATGCCGCCCTTGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCAGGAAATATGTTTATTAAGATTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAAAGACCGGTGATGTAGTCTTCtgt                                                                                                    <  2:1776871/150‑1 (MQ=255)
                                                      gctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttg                                                                                                 >  2:4159519/1‑150 (MQ=255)
                                                       ctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgt                                                                                                <  2:3994023/150‑1 (MQ=255)
                                                       ctCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgt                                                                                                <  2:117662/150‑1 (MQ=255)
                                                        tCCACCATGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCtgtttgtt                                                                                               >  2:1634588/1‑150 (MQ=255)
                                                         ccaccaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCTGTTCGGCTTTGAACTGGCCGGAATTATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTTAATCAACCGTAACGACCGGTGATGTAGTCTTCTgtttgttt                                                                                              <  1:1410860/150‑1 (MQ=255)
                                                         ccaccaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCAAAACGCACTCCTTAATCAACCGTAACGACCGGTGATGTAGTCTTCTgtttgttt                                                                                              <  1:544357/150‑1 (MQ=255)
                                                          caccaTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttc                                                                                             >  1:1975461/1‑150 (MQ=255)
                                                             caTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttcttc                                                                                          >  1:1002386/1‑150 (MQ=255)
                                                             caTGCCGCCCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTttcttc                                                                                          <  2:1017302/150‑1 (MQ=255)
                                                                   gccCATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGc                                                                                    >  1:1553410/1‑150 (MQ=255)
                                                                      cATGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTg                                                                                 <  2:285461/150‑1 (MQ=255)
                                                                        tGGTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGt                                                                               <  2:956447/150‑1 (MQ=255)
                                                                         ggTCATCACCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTg                                                                              <  2:537244/150‑1 (MQ=255)
                                                                           tcatcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGaa                                                                            >  1:4107798/1‑150 (MQ=255)
                                                                            catcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAAc                                                                           <  2:3892010/150‑1 (MQ=255)
                                                                            catcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAAc                                                                           <  2:2127116/150‑1 (MQ=255)
                                                                            catcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAAc                                                                           >  1:1382790/1‑150 (MQ=255)
                                                                             atcaCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACa                                                                          <  2:3029593/150‑1 (MQ=255)
                                                                               caCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGa                                                                        >  1:1451112/1‑150 (MQ=255)
                                                                                aCCTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGAt                                                                       >  2:1293951/1‑150 (MQ=255)
                                                                                 ccTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATc                                                                      <  2:608881/150‑1 (MQ=255)
                                                                                  cTGCGTGCGGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCg                                                                     <  1:2340582/150‑1 (MQ=255)
                                                                                   tgcgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGt                                                                    >  1:242265/1‑150 (MQ=255)
                                                                                   tgcgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGt                                                                    <  2:2483975/150‑1 (MQ=255)
                                                                                   tgcgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGAGGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGACCAGATCGt                                                                    >  2:3530887/1‑150 (MQ=255)
                                                                                    gcgtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTc                                                                   <  2:2374412/150‑1 (MQ=255)
                                                                                      gtgcgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCg                                                                 >  2:2899066/1‑150 (MQ=255)
                                                                                         cgGATACTTTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCgtgt                                                              <  1:2682493/150‑1 (MQ=255)
                                                                                           gATACTGTCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTg                                                            <  2:2202600/150‑1 (MQ=255)
                                                                                                tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                       >  1:4198475/1‑150 (MQ=255)
                                                                                                tttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGaa                                                       >  1:4183780/1‑150 (MQ=255)
                                                                                                 ttCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTCACGACCGGTGATGTAGTCGTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAAc                                                      >  2:3598249/1‑150 (MQ=255)
                                                                                                  tCCAGTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACt                                                     <  2:1451112/150‑1 (MQ=255)
                                                                                                      gTTCGGCGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACtcaat                                                 <  1:3332345/150‑1 (MQ=255)
                                                                                                         cGGCGTTTAACTGGCCGTAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCTTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTTGTTAACAGATCTTCCGTGTTGCTTAACtcaatcaa                                              <  2:3289615/150‑1 (MQ=255)
                                                                                                            cGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTc                                           >  2:155193/1‑150 (MQ=255)
                                                                                                            cGTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTc                                           >  1:3157926/1‑150 (MQ=255)
                                                                                                             gTTGAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCg                                          >  2:4250350/1‑150 (MQ=255)
                                                                                                                gAACTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGccc                                       <  2:3429107/150‑1 (MQ=255)
                                                                                                                 aaCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCa                                      <  1:3017931/150‑1 (MQ=255)
                                                                                                                 aaCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCa                                      <  2:1570753/150‑1 (MQ=255)
                                                                                                                  aCTGGCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAg                                     <  2:1484909/150‑1 (MQ=255)
                                                                                                                      gCCGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTAc                                 >  2:4121026/1‑150 (MQ=255)
                                                                                                                       ccGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACa                                >  2:3785989/1‑150 (MQ=255)
                                                                                                                        cGGAAATATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACAt                               >  1:861532/1‑150 (MQ=255)
                                                                                                                            aaTATGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAAc                           <  2:2595901/150‑1 (MQ=255)
                                                                                                                               atGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGcc                        >  1:622343/1‑150 (MQ=255)
                                                                                                                                tGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCg                       >  2:271352/1‑150 (MQ=255)
                                                                                                                                tGTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCg                       <  2:3100735/150‑1 (MQ=255)
                                                                                                                                 gTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgt                      <  2:118928/150‑1 (MQ=255)
                                                                                                                                 gTTTATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgt                      >  1:1668449/1‑150 (MQ=255)
                                                                                                                                  tttATTAAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgtg                     >  1:1205334/1‑150 (MQ=255)
                                                                                                                                   ttattaAGATTTAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACGGGTGTTGTATTCTTCTGTTTGTTTCTTCGCTGGTTTGGTGAACAGATCGTCGGTTTTGCTGAACTCAATCAATTCGCCCTGGTACATACACGCCgtgt                    <  1:1634588/150‑1 (MQ=255)
                                                                                                                                   ttattaAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCgtgt                    <  2:80042/150‑1 (MQ=255)
                                                                                                                                    tattaAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGtgtg                   <  2:92427/150‑1 (MQ=255)
                                                                                                                                      ttaAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGt                 >  1:918321/1‑150 (MQ=255)
                                                                                                                                       taAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGTTGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAAGAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTc                >  1:2896701/1‑150 (MQ=255)
                                                                                                                                       taAGATTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTc                <  2:3692138/150‑1 (MQ=255)
                                                                                                                                          gtttGAGACTGTCCATAACGCACTCCTGAATCAACCTTAACGACCGGTGTTGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGTTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGa             <  1:3493447/148‑1 (MQ=255)
                                                                                                                                           ttttagaCTGTCAATAACGCACTCCTGAATCAACCGTAACGACCGGTTATTTTTTCTTATGTTTGTTTCTTCGGTGGCTTGGTTAACAGATCGTCGGTGTTGCTTAACTCAATAAATTCGCCCAGTTAAATAAACGCCGTGTGGTCGGaa            <  1:427918/146‑1 (MQ=255)
                                                                                                                                           aTTGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaa            <  2:263023/150‑1 (MQ=255)
                                                                                                                                           aTTGAGACTGCCCATACCGCACTCCTTAAACAACCGTAACGACCGGTTATGTAGTCTTCTGTTTGTTTCTGCGCTGGGTTGGTGAACTGATCGTCCGTGTTGCTGAACTCAATCATTTCGCCCTTTTACATAAACGCCGTGTGGTCGGaa            <  1:1930156/150‑1 (MQ=255)
                                                                                                                                           aTTGAGACTGCCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaa            >  2:1253201/1‑150 (MQ=255)
                                                                                                                                           aTTGAGACCCCCCATACCGCACTCCTCAATAAACCGAAACGACCGGTGATGTAGTCTTCTTTTTGTTTCTTCGCTGGCTTGGTGAACATATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCTGGTACATAAACGCCGTGTGGTCGGaa            <  1:1829525/150‑1 (MQ=255)
                                                                                                                                            ttGAGACTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGaac           >  1:2955320/1‑150 (MQ=255)
                                                                                                                                                gaCTGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAAcg       >  2:231697/1‑150 (MQ=255)
                                                                                                                                                   tGTCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACgcgc    >  1:2027970/1‑150 (MQ=255)
                                                                                                                                                     tCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAg  >  1:3809482/1‑150 (MQ=255)
                                                                                                                                                     tCCATAACGCACTCCTGAATCAACCGTAACGACCGGTGATGTAGTCTTCTGTTTGTTTCTTCGCTGGCTTGGTGAACAGATCGTCCGTGTTGCTGAACTCAATCAATTCGCCCAGGTACATAAACGCCGTGTGGTCGGAACAACGCGCAg  >  1:2265383/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: