Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,599,759 C→T 5.5% intergenic (‑101/+170) livJ ← / ← rpoH branched‑chain amino acid ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,599,7590CT 5.5% 232.5 / 2.2 72intergenic (‑101/+170)livJ/rpoHbranched‑chain amino acid ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (2/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

CCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGT  >  NC_000913/3599613‑3599908
                                                                                                                                                  |                                                                                                                                                     
ccAGCGCGATACATCCTTCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                    <  1:772006/150‑1 (MQ=255)
ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCCTTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                    >  1:1725937/1‑150 (MQ=255)
ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                    <  1:1017297/150‑1 (MQ=255)
 cAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTg                                                                                                                                                   >  1:674889/1‑150 (MQ=255)
  ctcgcgATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATGCTTCTGGTTATTACGGGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGc                                                                                                                                                  <  2:1841984/148‑1 (MQ=255)
       gATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTACAGGTTATTACGTGTGTTGTGATGTTGTATTTCAGCACTTTATTTCGATTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTaaaa                                                                                                                                             >  1:586102/1‑150 (MQ=255)
        ttACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAAc                                                                                                                                            <  2:1890773/149‑1 (MQ=255)
         tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                           >  2:1129502/1‑150 (MQ=255)
         tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                           <  2:1780552/150‑1 (MQ=255)
           cATCCTGCCAGTAACGCTTTACCCTTTATGATCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTAAAGCacg                                                                                                                                         <  2:1953080/150‑3 (MQ=255)
            aTCCTGCCAGTAACGCTTTACCCTTTAAGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAc                                                                                                                                        <  1:1102821/150‑1 (MQ=255)
             tCCTTCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAcc                                                                                                                                       <  1:1796575/150‑1 (MQ=255)
              ccTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAccc                                                                                                                                      >  2:480804/1‑150 (MQ=255)
               cTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTAAAGCacgggta                                                                                                                                     <  2:660142/150‑7 (MQ=255)
                 gCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAt                                                                                                                                   >  2:328962/1‑150 (MQ=255)
                  ccAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAtt                                                                                                                                  >  1:174689/1‑150 (MQ=255)
                    aGTAACTCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAtttt                                                                                                                                <  1:1522783/150‑1 (MQ=255)
                    aGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAtttt                                                                                                                                <  1:865298/150‑1 (MQ=255)
                     gTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAttttt                                                                                                                               <  2:927182/150‑1 (MQ=255)
                                cccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTACAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                    >  1:76458/1‑150 (MQ=255)
                                 ccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGaa                                                                                                                   >  1:1642510/1‑150 (MQ=255)
                                 ccTTTATGTTCATCCTGAGAATCCACATTCTTCTGGTTTTTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGaa                                                                                                                   <  2:1233488/150‑1 (MQ=255)
                                  cTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAAt                                                                                                                  <  2:548522/150‑1 (MQ=255)
                                    ttATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTTTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACAACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAg                                                                                                                <  2:83356/150‑1 (MQ=255)
                                    ttATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAg                                                                                                                >  2:193042/1‑150 (MQ=255)
                                      aTGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTATATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATAGGAATATCt                                                                                                              >  1:753394/1‑150 (MQ=255)
                                        gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                            >  1:1766542/1‑150 (MQ=255)
                                        gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                            >  2:557185/1‑150 (MQ=255)
                                              ccTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGAc                                                                                                      <  2:955752/150‑1 (MQ=255)
                                               cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                                     <  2:1746006/150‑1 (MQ=255)
                                                     aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                               <  1:1983816/150‑1 (MQ=255)
                                                     aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                               >  2:529046/1‑150 (MQ=255)
                                                       ccccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatta                                                                                             <  2:29117/150‑1 (MQ=255)
                                                        cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                            >  2:513387/1‑150 (MQ=255)
                                                        cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                            >  1:834761/1‑150 (MQ=255)
                                                                tctGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGc                                                                                    >  1:2025711/1‑150 (MQ=255)
                                                                   ggTTATTACGTGTGTTGTGATGTTGTTATTCAGCACTTTATTTCGTTTTATGAATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATATCTATTTTGACAGTTTATTAACAATCTGCGTg                                                                                 >  1:558171/1‑150 (MQ=255)
                                                                      tattaCGTGTGTTGTGATGGTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTgggg                                                                              <  2:720480/150‑1 (MQ=255)
                                                                           cGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTgg                                                                         >  2:1177117/1‑150 (MQ=255)
                                                                           cGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTCCTTTAGCAGCATACTCTGCTAAAACATACCCGATATTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTgt                                                                         >  2:1605632/1‑149 (MQ=255)
                                                                               tgttgtGATGTTGTTTTTCAGCACTTTATTTCGTGTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGtt                                                                     <  2:1487849/150‑1 (MQ=255)
                                                                               tgttggGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGtt                                                                     <  2:1444140/150‑1 (MQ=255)
                                                                                  tgtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTg                                                                  >  1:1606193/1‑150 (MQ=255)
                                                                                  tgtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTg                                                                  >  2:1164816/1‑150 (MQ=255)
                                                                                      aTGGTGGTTTTCAGCACTTTATTTCGTTATATGCATGCCGCCCCGGGCTTTCGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATAATGGAATAGCTATTTTGACAGTGTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCgg                                                              <  1:967220/150‑1 (MQ=255)
                                                                                       tgttgtTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGa                                                             >  2:356108/1‑150 (MQ=255)
                                                                                          tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                          <  2:1766541/150‑1 (MQ=255)
                                                                                           gttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATTCTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGGTTATTTACAATCTGCGTGGGGATTGGCGTTTTGCAGGAgggg                                                         >  2:1357586/1‑150 (MQ=255)
                                                                                           gttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAgggg                                                         >  2:725630/1‑150 (MQ=255)
                                                                                                 cAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTg                                                   <  2:174689/150‑1 (MQ=255)
                                                                                                   gCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGat                                                 >  1:667507/1‑150 (MQ=255)
                                                                                                   gCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGat                                                 <  2:1916486/150‑1 (MQ=255)
                                                                                                     aCTTTATTTCGTTTTATGCATGCCTAACCGTTCTTTAGCAGCATACTCTGGTAAAACATAGCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatta                                               <  2:393031/150‑1 (MQ=255)
                                                                                                      cTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGattat                                              >  1:678884/1‑150 (MQ=255)
                                                                                                             tCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACagag                                       >  1:223456/1‑150 (MQ=255)
                                                                                                                       cATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCaaaaaa                             >  2:1573591/1‑150 (MQ=255)
                                                                                                                        aTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaac                            <  1:2041370/150‑1 (MQ=255)
                                                                                                                        aTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGGTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaac                            <  2:1666628/150‑1 (MQ=255)
                                                                                                                        aTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAAGAGCTATTTTGACAGTTTATTAACAATCTGCGGGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaac                            <  1:341102/150‑1 (MQ=255)
                                                                                                                         tGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaaca                           >  1:1165419/1‑150 (MQ=255)
                                                                                                                         tGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaaca                           >  2:659207/1‑150 (MQ=255)
                                                                                                                              aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                      >  1:1400093/1‑150 (MQ=255)
                                                                                                                              aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                      >  1:462519/1‑150 (MQ=255)
                                                                                                                              aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCAGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                      >  2:1424464/1‑150 (MQ=255)
                                                                                                                                 cGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAAcccc                   >  1:522183/1‑150 (MQ=255)
                                                                                                                                  gTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCg                  <  1:41763/150‑1 (MQ=255)
                                                                                                                                    gttttAGCAGAGTATGCTGCTAAAACATACCCGTTTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCggc                >  1:198555/3‑149 (MQ=255)
                                                                                                                                    gttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCggc                >  1:1740976/3‑149 (MQ=255)
                                                                                                                                     cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc               >  1:776408/1‑150 (MQ=255)
                                                                                                                                      tttAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGAAATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAct              >  1:1073814/1‑150 (MQ=255)
                                                                                                                                         agcagcATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACtctc           <  2:1606192/150‑1 (MQ=255)
                                                                                                                                          gcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCa          >  2:638764/1‑150 (MQ=255)
                                                                                                                                          gcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCa          >  1:388048/1‑150 (MQ=255)
                                                                                                                                              caTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCa      <  2:674889/150‑1 (MQ=255)
                                                                                                                                                 aCTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGATTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAggg   <  1:328962/150‑1 (MQ=255)
                                                                                                                                                  ctctGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGATGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGt  >  1:1426383/1‑150 (MQ=255)
                                                                                                                                                  ctctGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGt  >  1:734292/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGT  >  NC_000913/3599613‑3599908

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