Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,316 G→A 6.0% intergenic (‑185/+267) 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,3160GA 6.0% 322.9 / 6.6 101intergenic (‑185/+267)iscR/trmJDNA‑binding transcriptional dual regulator IscR/tRNA Cm32/Um32 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (41/53);  new base A (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 = 9.80e-01

Only 100 of 106 total aligned reads displayed.

TTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCA  >  NC_000913/2662167‑2662465
                                                                                                                                                     |                                                                                                                                                     
ttttttATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATg                                                                                                                                                       >  1:827929/1‑150 (MQ=255)
   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                                                                                                                                  <  1:491235/150‑1 (MQ=255)
                     gTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGc                                                                                                                                  >  2:173735/1‑150 (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:1019259/1‑150 (MQ=255)
                          gCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTgg                                                                                                                             <  2:362136/150‑1 (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                                                                                                                     >  1:608026/1‑150 (MQ=255)
                                  gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                     <  2:508253/150‑1 (MQ=255)
                                  gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                     <  1:648778/150‑1 (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                                                                                                      <  1:1089246/150‑1 (MQ=255)
                                                 ttCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGa                                                                                                      <  2:886621/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                                                                                             >  1:349983/1‑150 (MQ=255)
                                                          aaaTTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCa                                                                                             >  2:787971/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: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)
                                                                                                                            gTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATACGACCTACGGTTCTg                           >  1:613761/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                     
TTTTTTATTGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCA  >  NC_000913/2662167‑2662465

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: