Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,599,756 A→G 5.5% intergenic (‑98/+173) 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,7560AG 5.5% 269.8 / 5.0 91intergenic (‑98/+173)livJ/rpoHbranched‑chain amino acid ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor
Reads supporting (aligned to +/- strand):  ref base A (49/36);  new base G (3/2);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TGAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATC  >  NC_000913/3599607‑3599902
                                                                                                                                                     |                                                                                                                                                  
tGAATGCCAGCGCGATACATCCTGCCAGTAACGATTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGGTGTTTTTCAGGACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAgcagca                                                                                                                                                    <  2:81932/150‑1 (MQ=255)
 gAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCAt                                                                                                                                                   >  1:213775/1‑150 (MQ=255)
     gCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActc                                                                                                                                               <  1:1615283/150‑1 (MQ=255)
      ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                              <  2:1664408/150‑1 (MQ=255)
      ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                              <  2:356410/150‑1 (MQ=255)
      ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                              >  1:650870/1‑150 (MQ=255)
          cgcgATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTa                                                                                                                                          <  1:1696480/150‑1 (MQ=255)
             gATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGCGACTCCCCATTCTTCTGGGTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTaaag                                                                                                                                       <  2:994240/150‑2 (MQ=255)
             gATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTaaaa                                                                                                                                       >  1:385979/1‑150 (MQ=255)
             gATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTaaaa                                                                                                                                       <  1:551937/150‑1 (MQ=255)
              aTACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAAc                                                                                                                                      >  2:860377/1‑150 (MQ=255)
              aTACAACCTGCCAGTAACGCAATACCCTTTCTGTTCATCCTGAGAATCCCCATTCTTCTGGATATTACGTGTGTTCTGATGTTGTTTTTCAGCACTTTATTTCGATTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAAc                                                                                                                                      <  1:361144/150‑1 (MQ=255)
                 cATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATa                                                                                                                                   >  2:393574/1‑150 (MQ=255)
                 cATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATa                                                                                                                                   >  1:1532301/1‑150 (MQ=255)
                  aTCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTAAAGCacgg                                                                                                                                  <  1:1212303/150‑4 (MQ=255)
                  aTCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGGTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAc                                                                                                                                  >  2:138897/1‑150 (MQ=255)
                      tGCCAGAACCGCTTTACCCTTTATGTTCATCCTCAGAATCCCCATTCATCTGGTTATTACGTGTGTTGTGATGTTGTTTATCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCAGCTAAAACATACCCGa                                                                                                                              >  2:1736360/1‑150 (MQ=255)
                         cAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAttt                                                                                                                           >  2:891431/1‑150 (MQ=255)
                            tAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTa                                                                                                                        <  2:1285258/150‑1 (MQ=255)
                             aaCGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTAt                                                                                                                       >  1:823211/1‑150 (MQ=255)
                                    tACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTg                                                                                                                <  1:1187970/150‑1 (MQ=255)
                                      cccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                              <  2:1211883/150‑1 (MQ=255)
                                           tATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGc                                                                                                         >  1:403741/1‑150 (MQ=255)
                                           tATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGc                                                                                                         >  2:1345061/1‑150 (MQ=255)
                                             tGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTa                                                                                                       >  1:1541715/1‑150 (MQ=255)
                                              gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                      >  1:636544/1‑150 (MQ=255)
                                              gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                      >  2:526382/1‑150 (MQ=255)
                                              gTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAt                                                                                                      >  1:516279/1‑150 (MQ=255)
                                                   tacTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGa                                                                                                 >  1:11387/3‑150 (MQ=255)
                                                    ccTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGAc                                                                                                <  2:1532301/150‑1 (MQ=255)
                                                    ccTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGGTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGAc                                                                                                <  1:109420/150‑1 (MQ=255)
                                                     cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                               >  2:158818/1‑150 (MQ=255)
                                                     cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                               >  1:500467/1‑150 (MQ=255)
                                                     cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACACGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                               >  1:759282/1‑150 (MQ=255)
                                                      tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAg                                                                                              >  2:134648/1‑150 (MQ=255)
                                                       gagaATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGt                                                                                             >  2:1112972/1‑150 (MQ=255)
                                                         gaATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGttt                                                                                           >  1:196399/1‑150 (MQ=255)
                                                          aaTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTtta                                                                                          >  2:1774134/1‑150 (MQ=255)
                                                          aaTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTtta                                                                                          <  2:679417/150‑1 (MQ=255)
                                                           aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                         >  2:1642051/1‑150 (MQ=255)
                                                            tCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatt                                                                                        <  2:1187097/150‑1 (MQ=255)
                                                             ccccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatta                                                                                       >  1:1052920/1‑150 (MQ=255)
                                                             ccccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatta                                                                                       <  2:1513795/150‑1 (MQ=255)
                                                             ccccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatta                                                                                       <  1:1032581/150‑1 (MQ=255)
                                                               ccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAAc                                                                                     <  2:213775/150‑1 (MQ=255)
                                                                      tctGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGc                                                                              >  1:626752/1‑150 (MQ=255)
                                                                          gTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTgg                                                                          <  2:650870/150‑1 (MQ=255)
                                                                           ttattaCGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTggg                                                                         <  2:1464932/150‑1 (MQ=255)
                                                                               taCGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGAtt                                                                     >  2:558438/1‑150 (MQ=255)
                                                                                aCGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTg                                                                    <  2:195285/150‑1 (MQ=255)
                                                                                   tgtgtTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCg                                                                 <  2:295634/150‑1 (MQ=255)
                                                                                      gttgtGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGTCGttt                                                              >  2:710355/1‑150 (MQ=255)
                                                                                      gttgtGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGttt                                                              >  2:157409/1‑150 (MQ=255)
                                                                                      gttgtGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGttt                                                              <  2:850182/150‑1 (MQ=255)
                                                                                        tgtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTg                                                            <  2:636414/150‑1 (MQ=255)
                                                                                         gtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGc                                                           >  1:85831/1‑150 (MQ=255)
                                                                                          tgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGcc                                                          <  1:575281/150‑1 (MQ=255)
                                                                                             tgttgtTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGa                                                       >  1:1360544/1‑150 (MQ=255)
                                                                                             tgttgtTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGa                                                       >  1:1247238/1‑150 (MQ=255)
                                                                                              gttgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAg                                                      >  1:1143876/1‑150 (MQ=255)
                                                                                              gttgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAg                                                      >  2:606146/1‑150 (MQ=255)
                                                                                                tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                    <  1:1573698/150‑1 (MQ=255)
                                                                                                tgttTTACAGCACTTTATTTCGTTTTATGCATGACTACCCGGGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                    <  2:226942/150‑1 (MQ=255)
                                                                                                 gttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCCTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGGGGGGATTGGCGTTTTGCCGGAgggg                                                   <  2:558622/150‑1 (MQ=255)
                                                                                                      tCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAAtt                                              >  1:862084/1‑150 (MQ=255)
                                                                                                      tCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGACGGAGGGGAAAtt                                              >  1:1640527/1‑150 (MQ=255)
                                                                                                         ggaCTTTATTTAGTTTTATGCATGACTACCCGTTCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGTATTGGCGTTTTGCCGGAGGGGAAATTGat                                           <  1:577616/148‑1 (MQ=255)
                                                                                                          cACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatt                                          <  1:466975/150‑1 (MQ=255)
                                                                                                               tATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTAc                                     <  1:1716417/150‑1 (MQ=255)
                                                                                                                  ttCGTTTTATGCATGTCTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACaga                                  >  1:1598870/1‑150 (MQ=255)
                                                                                                                  ttCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACaga                                  >  1:1647352/1‑150 (MQ=255)
                                                                                                                  ttCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACaga                                  >  2:1198342/1‑150 (MQ=255)
                                                                                                                   tCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTCTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACagag                                 <  1:891431/150‑1 (MQ=255)
                                                                                                                     gTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGc                               >  2:849996/1‑150 (MQ=255)
                                                                                                                     gTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTCATAAAACATACCCGATTTTTATCATATTGGAATAGCTATTTTGACAGTTTATTAACAATCAGCGTGGGGATTGGCGTTTTGCCGGAGTGGATAATGATTATAACACAGGc                               >  2:139837/1‑150 (MQ=255)
                                                                                                                      ttttAGGCATGACTCACCGGGCTTTAGCAGCATACGCTGCTAAAACATACCCGATTTTGATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTTCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGcc                              <  2:316108/150‑1 (MQ=255)
                                                                                                                       tttATGCCTGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGTAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGccc                             <  2:1799721/150‑1 (MQ=255)
                                                                                                                        ttATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCa                            <  2:673667/150‑1 (MQ=255)
                                                                                                                           tGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAATCATACCCGATATTTATGATATTGCAATAGCTATTTTGACAGTATATTAACAATCTGCGTGGGGATTGGAGTTTTGACTGAGGGGATATACATTTATACACAGTCCCaaaa                         >  1:287114/1‑150 (MQ=255)
                                                                                                                             cATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCaaaaaa                       >  1:812814/1‑150 (MQ=255)
                                                                                                                               tGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaaca                     <  2:186175/150‑1 (MQ=255)
                                                                                                                                    aCCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaaaac                >  1:1667972/1‑150 (MQ=255)
                                                                                                                                     cccGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAAcc               >  2:568639/1‑150 (MQ=255)
                                                                                                                                     cccGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGAATATTACAGAGGCACAAAAAACAAAAAcc               >  2:289490/1‑150 (MQ=255)
                                                                                                                                         tgttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCgg           >  1:1370590/4‑150 (MQ=255)
                                                                                                                                           ttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc         >  2:985026/2‑150 (MQ=255)
                                                                                                                                           ttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc         >  2:978654/2‑150 (MQ=255)
                                                                                                                                           cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAATACAAAAACCCCGGAc         >  1:117543/1‑150 (MQ=255)
                                                                                                                                           cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGGGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc         <  2:830975/150‑1 (MQ=255)
                                                                                                                                             ttAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGActc       >  1:564159/1‑150 (MQ=255)
                                                                                                                                             ttAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGActc       <  1:134648/150‑1 (MQ=255)
                                                                                                                                              ttgcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGActct      >  2:748109/3‑150 (MQ=255)
                                                                                                                                              tAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGActct      <  1:885047/150‑1 (MQ=255)
                                                                                                                                                gcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCa    >  1:1166413/1‑150 (MQ=255)
                                                                                                                                                gcagcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGGTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCa    >  2:943956/1‑150 (MQ=255)
                                                                                                                                                  agcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATc  >  1:541314/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATC  >  NC_000913/3599607‑3599902

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: