Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,662,324 A→C 6.7% intergenic (‑193/+259) 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,3240AC 6.7% 322.5 / 10.9 102intergenic (‑193/+259)iscR/trmJDNA‑binding transcriptional dual regulator IscR/tRNA Cm32/Um32 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (40/55);  new base C (3/4);  total (43/59)
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.97e-01

Only 100 of 107 total aligned reads displayed.

TGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGT  >  NC_000913/2662175‑2662473
                                                                                                                                                     |                                                                                                                                                     
tGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGa                                                                                                                                                       <  2:1405918/150‑1 (MQ=255)
tGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGa                                                                                                                                                       <  1:1460340/150‑1 (MQ=255)
 gTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCATCTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGaa                                                                                                                                                      >  1:1356983/1‑150 (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                                                                                                                                            >  2:574357/1‑150 (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: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:173735/150‑1 (MQ=255)
                     tAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGt                                                                                                                                  >  1:252626/1‑150 (MQ=255)
                        cAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCt                                                                                                                               <  2:517683/150‑1 (MQ=255)
                          gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                             <  1:648778/150‑1 (MQ=255)
                          gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                             >  1:608026/1‑150 (MQ=255)
                          gggTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGt                                                                                                                             <  2:508253/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                                                                                                              <  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)
                                                    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                                                                                                           <  1:1260103/136‑1 (MQ=255)
                                                          cAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGt                                                                                             >  2:417358/1‑150 (MQ=255)
                                                           aaCTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                            <  2:373564/150‑1 (MQ=255)
                                                            aCTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCg                                                                                           <  2:1572060/150‑1 (MQ=255)
                                                                tttACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                                                                       <  2:1221946/150‑1 (MQ=255)
                                                                  tACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                     <  1:1091744/150‑1 (MQ=255)
                                                                  tACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                     <  2:1356983/150‑1 (MQ=255)
                                                                  tACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAg                                                                                     >  1:632157/1‑150 (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)
                                                                                                            tttttATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                           >  1:1518671/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)
                                                                                                                                 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             >  1:729487/1‑150 (MQ=255)
                                                                                                                                            tGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGc           <  1:95075/150‑1 (MQ=255)
                                                                                                                                             gcggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCa          <  1:1361956/150‑1 (MQ=255)
                                                                                                                                              cggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAt         <  1:584158/150‑1 (MQ=255)
                                                                                                                                               ggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt        >  1:1470616/1‑150 (MQ=255)
                                                                                                                                               ggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt        <  2:1313663/150‑1 (MQ=255)
                                                                                                                                               ggcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCAtt        <  1:1123061/150‑1 (MQ=255)
                                                                                                                                                gcgTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa       <  1:651647/150‑1 (MQ=255)
                                                                                                                                                   tGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgc    <  1:787971/150‑1 (MQ=255)
                                                                                                                                                    gAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTAcgcg   >  2:1465288/1‑150 (MQ=255)
                                                                                                                                                     aaCGCAATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCCACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGt  <  2:882767/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTAAAAGTGGGGGTATTGCATAGCAGGGTCAAGTCTGACATTCCCGAGTAAATTGGTCAACTATTTACTTGACTGATTTAGTCGGGTATTTAACCTTCAGTGCCATTTTTTTATCGTGGCGTTTGTACCTGTTGTCGGATGCGGCGTGAACGCCATATCCGCCCTGGGGTTCTGTACACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGT  >  NC_000913/2662175‑2662473

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: