Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,317 C→T 5.9% intergenic (‑186/+266) iscR ← / ← trmJ DNA‑binding transcriptional dual regulator IscR/tRNA Cm32/Um32 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,662,3170CT 5.9% 317.3 / 6.3 101intergenic (‑186/+266)iscR/trmJDNA‑binding transcriptional dual regulator IscR/tRNA Cm32/Um32 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/53);  new base T (3/3);  total (44/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 106 total aligned reads displayed.

TTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAT  >  NC_000913/2662170‑2662466
                                                                                                                                                   |                                                                                                                                                     
tttATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATgcgg                                                                                                                                                     >  1:795714/1‑150 (MQ=255)
 ttATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATgcggc                                                                                                                                                    <  2:350335/150‑1 (MQ=255)
  tATTGTAAAAGTGGTGGTTTTGCATAGGAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATgcggcg                                                                                                                                                   <  2:305495/150‑1 (MQ=255)
   aTTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGt                                                                                                                                                  >  1:1179289/1‑150 (MQ=255)
    ttGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTg                                                                                                                                                 <  1:568516/150‑1 (MQ=255)
     tGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGa                                                                                                                                                <  2:1405918/150‑1 (MQ=255)
     tGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGa                                                                                                                                                <  1:1460340/150‑1 (MQ=255)
      gTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCATCTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGaa                                                                                                                                               >  1:1356983/1‑150 (MQ=255)
      gTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTTGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGaa                                                                                                                                               >  1:773782/1‑150 (MQ=255)
      gTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGttca                                                                                                                                               <  2:382129/150‑4 (MQ=255)
        aaaaGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACg                                                                                                                                             <  2:1402160/150‑1 (MQ=255)
         aaaGTGGGGGTATTGCATAGCAGGGTCAAGGCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGc                                                                                                                                            <  1:695230/150‑1 (MQ=255)
          aaGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGcc                                                                                                                                           <  2:678862/150‑1 (MQ=255)
            gTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCat                                                                                                                                         <  1:720383/150‑1 (MQ=255)
             tGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCata                                                                                                                                        <  1:934908/150‑1 (MQ=255)
             tGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCata                                                                                                                                        <  1:755894/150‑1 (MQ=255)
                gggTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATcc                                                                                                                                     >  1:373237/1‑150 (MQ=255)
                  gTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGc                                                                                                                                   >  2:173735/1‑150 (MQ=255)
                  gTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGc                                                                                                                                   <  1:491235/150‑1 (MQ=255)
                   tATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATATGGCGTTCACGCCGCATCCgac                                                                                                                                  <  2:1044889/150‑3 (MQ=255)
                    aTTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGccc                                                                                                                                 >  1:882767/1‑150 (MQ=255)
                      tGCATAGCAGGGGCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATATGGCGTTCACGCCGCATCCgacaac                                                                                                                               <  1:1359269/150‑6 (MQ=255)
                       gCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTgg                                                                                                                              <  2:362136/150‑1 (MQ=255)
                       gCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTgg                                                                                                                              >  2:1019259/1‑150 (MQ=255)
                       gCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTgg                                                                                                                              >  2:591050/1‑150 (MQ=255)
                         aTAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTgggg                                                                                                                            <  1:1243883/150‑1 (MQ=255)
                          tAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGt                                                                                                                           >  1:252626/1‑150 (MQ=255)
                          tAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGt                                                                                                                           <  1:173735/150‑1 (MQ=255)
                             cAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCt                                                                                                                        <  2:517683/150‑1 (MQ=255)
                               gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                      <  2:508253/150‑1 (MQ=255)
                               gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                      <  1:648778/150‑1 (MQ=255)
                               gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                      >  1:608026/1‑150 (MQ=255)
                                ggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTa                                                                                                                     >  2:1361956/1‑150 (MQ=255)
                                 gTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTac                                                                                                                    >  1:274034/1‑150 (MQ=255)
                                   cAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTacac                                                                                                                  >  2:95075/1‑150 (MQ=255)
                                    aaGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACt                                                                                                                 >  1:281453/1‑150 (MQ=255)
                                      gTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGt                                                                                                               <  2:1179289/150‑1 (MQ=255)
                                        cTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAg                                                                                                             <  1:194556/150‑1 (MQ=255)
                                         tGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAgg                                                                                                            >  1:186270/1‑150 (MQ=255)
                                          gACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGc                                                                                                           >  2:1036597/1‑150 (MQ=255)
                                           aCATTCCCGAGTAAATTGGGCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGcc                                                                                                          <  2:794465/150‑1 (MQ=255)
                                            cATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCt                                                                                                         <  1:1078627/150‑1 (MQ=255)
                                              ttCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGa                                                                                                       <  2:886621/150‑1 (MQ=255)
                                              ttCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGa                                                                                                       <  1:1089246/150‑1 (MQ=255)
                                               tCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGAt                                                                                                      >  2:651647/1‑150 (MQ=255)
                                                cccGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATa                                                                                                     >  1:809937/1‑150 (MQ=255)
                                                  cGACTAAATCAGTCAAGTAAATAGTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAg                                                                                                   >  2:584158/1‑150 (MQ=38)
                                                       aaaTTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCa                                                                                              >  2:787971/1‑150 (MQ=255)
                                                       aaaTTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCa                                                                                              >  1:349983/1‑150 (MQ=255)
                                                       aaaTTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCa                                                                                              <  1:538185/150‑1 (MQ=255)
                                                         aTTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCAtt                                                                                            <  2:281453/150‑1 (MQ=255)
                                                          ttGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTa                                                                                           >  2:837823/1‑150 (MQ=255)
                                                            ggTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTAcg                                                                                         >  2:150655/1‑150 (MQ=255)
                                                             gTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTAcgc                                                                                        <  1:574357/150‑1 (MQ=255)
                                                               cAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATaa                                                                                                    >  2:1260103/1‑136 (MQ=255)
                                                               cAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATaa                                                                                                    <  1:1260103/136‑1 (MQ=255)
                                                               cAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGt                                                                                      >  2:417358/1‑150 (MQ=255)
                                                                aaCTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                     <  2:373564/150‑1 (MQ=255)
                                                                     tttACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                                                                <  2:1221946/150‑1 (MQ=255)
                                                                       tACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                              <  2:1356983/150‑1 (MQ=255)
                                                                       tACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                              >  1:632157/1‑150 (MQ=255)
                                                                       tACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                              <  1:1091744/150‑1 (MQ=255)
                                                                            gACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAAc                                                                         >  1:823826/1‑150 (MQ=255)
                                                                              cTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACgg                                                                       >  1:818029/1‑150 (MQ=255)
                                                                                    tAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCg                                                                 <  2:274034/150‑1 (MQ=255)
                                                                                      ttCGGGTATTTAACCTTCAGTGGCATTTTTTTATCGTGGCGTTTTTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGa                                                               <  2:602494/149‑1 (MQ=255)
                                                                                        cGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATg                                                             <  2:364089/150‑1 (MQ=255)
                                                                                         gggTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                            >  1:1313663/1‑150 (MQ=255)
                                                                                          ggTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcg                                                           <  1:1057588/150‑1 (MQ=255)
                                                                                            tATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATgcggc                                                         <  2:1282412/150‑1 (MQ=255)
                                                                                                tAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACAGCTGTCGGATGCGGCGTGa                                                     <  2:236344/150‑1 (MQ=255)
                                                                                                 aaCCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGGGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGaa                                                    <  2:1326671/150‑1 (MQ=255)
                                                                                                     ttCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGcc                                                <  1:480027/150‑1 (MQ=255)
                                                                                                      tCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTTAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCt                                               <  1:182889/150‑1 (MQ=255)
                                                                                                      tCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCt                                               <  2:823826/150‑1 (MQ=255)
                                                                                                       cAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCtt                                              <  2:818029/150‑1 (MQ=255)
                                                                                                           gCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATcc                                          <  1:1036597/150‑1 (MQ=255)
                                                                                                           gCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATcc                                          >  1:994983/1‑150 (MQ=255)
                                                                                                              aTTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGAc                                       >  2:1511587/1‑150 (MQ=255)
                                                                                                                ttttttATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCt                                     <  2:795714/150‑1 (MQ=255)
                                                                                                                ttttttATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCt                                     >  1:1364756/1‑150 (MQ=255)
                                                                                                                     tATCGTGGCGTTTGGAGCTGTTGTCGGATGAGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGGCTGATATGATGCATTACGCGTCGCATCAGGCTACGGCTGTCGGTTGCGGCGTGAACGCCTTATGCGACCTTCGGt                                <  2:827929/150‑1 (MQ=17)
                                                                                                                       tCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTc                              <  2:349983/150‑1 (MQ=255)
                                                                                                                            gCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                         <  2:1350577/150‑1 (MQ=255)
                                                                                                                                   ccccTGTTGTCGCCGCCGGCGGGCACGCCATAGCCGCCCTGGGGTTCCGTACCCTGCAGCCCTGCTCAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGCTGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAg                  <  2:299712/148‑1 (MQ=14)
                                                                                                                                    aCCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAgg                 >  1:981099/1‑150 (MQ=255)
                                                                                                                                    aCCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAgg                 <  2:215472/150‑1 (MQ=255)
                                                                                                                                      cTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGcc               >  1:1414249/1‑150 (MQ=255)
                                                                                                                                      cTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGcc               <  1:591050/150‑1 (MQ=255)
                                                                                                                                       ccagggcGGATATGGCGTTCACGCCGCATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCt              >  1:653285/7‑150 (MQ=14)
                                                                                                                                        gttgtCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTg             <  2:186270/150‑1 (MQ=255)
                                                                                                                                        cagggcGGATATGGCGTTCACGCCGCATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTg             >  1:122100/6‑150 (MQ=14)
                                                                                                                                         agggcGGATATGGCGTTCACGCCGCATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGa            >  1:665044/5‑150 (MQ=14)
                                                                                                                                          tgtCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGAt           >  1:895237/1‑150 (MQ=255)
                                                                                                                                             cGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAg        <  2:809937/150‑1 (MQ=255)
                                                                                                                                              ggATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGa       >  1:1110966/1‑150 (MQ=255)
                                                                                                                                               gATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGAc      >  2:600110/1‑150 (MQ=255)
                                                                                                                                               gATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGAc      >  1:729487/1‑150 (MQ=255)
                                                                                                                                                  gcggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCa   <  1:1361956/150‑1 (MQ=255)
                                                                                                                                                   cggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAt  <  1:584158/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAT  >  NC_000913/2662170‑2662466

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