Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,861 C→T R488C (CGT→TGT)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,8610CT99.1% 410.6 / ‑5.7 111R488C (CGT→TGT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (48/62);  total (48/63)
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.96e-01

Only 100 of 111 total aligned reads displayed.

CGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGT  >  NC_000913/3823716‑3824003
                                                                                                                                                 |                                                                                                                                              
  gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCt                                                                                                                                         <  2:438467/151‑1 (MQ=255)
  gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCt                                                                                                                                         <  2:1129146/151‑1 (MQ=255)
  gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCt                                                                                                                                         <  1:421961/151‑1 (MQ=255)
    cAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCa                                                                                                                                       <  2:352902/151‑1 (MQ=255)
     aaaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCa                                                                                                                                       >  1:135470/1‑150 (MQ=255)
        ccGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCa                                                                                                                                   <  1:725297/151‑1 (MQ=255)
         cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCAt                                                                                                                                  >  2:458234/1‑151 (MQ=255)
         cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCAt                                                                                                                                  <  2:232887/151‑1 (MQ=255)
         cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCACCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCAt                                                                                                                                  >  2:256658/1‑151 (MQ=255)
           ttGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATg                                                                                                                                 <  1:622098/150‑1 (MQ=255)
            tGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGc                                                                                                                                >  2:954416/1‑150 (MQ=255)
             gAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCtt                                                                                                                              >  1:950286/1‑151 (MQ=255)
              aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCttt                                                                                                                             >  2:39973/1‑151 (MQ=255)
                aTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTgg                                                                                                                           <  2:121775/151‑1 (MQ=255)
                    ttATCGCTCCGGGCGATCGCCCGACTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGATAGCGCGCGCCCAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGGTTTTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGggtgg                                                                                                                       >  1:564524/1‑151 (MQ=255)
                     tACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGggtgg                                                                                                                       <  1:1139368/150‑1 (MQ=255)
                        cGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGCGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGc                                                                                                                   >  1:448704/1‑151 (MQ=255)
                         gCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGcc                                                                                                                  <  1:71864/151‑1 (MQ=255)
                         gCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGcc                                                                                                                  >  1:1372979/1‑151 (MQ=255)
                         gCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGcc                                                                                                                  >  1:551629/1‑151 (MQ=255)
                            ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                                                                >  2:344698/1‑150 (MQ=255)
                            ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTa                                                                                                               <  1:1376423/151‑1 (MQ=255)
                             cGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                                                                <  2:936298/149‑1 (MQ=255)
                             cGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTaa                                                                                                              <  2:867850/151‑1 (MQ=255)
                             cGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTa                                                                                                               >  1:22883/1‑150 (MQ=255)
                                     cGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaat                                                                                                       <  2:345489/150‑1 (MQ=255)
                                      gCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatga                                                                                                     <  2:175710/151‑1 (MQ=255)
                                       cccGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatgaa                                                                                                    <  2:1059135/151‑1 (MQ=255)
                                        ccGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGaaa                                                                                                   <  2:214738/151‑1 (MQ=255)
                                           aaTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATcc                                                                                                <  1:1099467/151‑1 (MQ=255)
                                            aTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATccc                                                                                               <  1:1292860/151‑1 (MQ=255)
                                              gCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                                                                              <  1:1171445/150‑1 (MQ=255)
                                              gCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                                                                              <  2:1345082/150‑1 (MQ=255)
                                               ccGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCa                                                                                            <  2:814976/151‑1 (MQ=255)
                                                cGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCACGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGc                                                                                             >  1:1013181/1‑149 (MQ=255)
                                                cGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGc                                                                                             >  2:1290055/1‑149 (MQ=255)
                                                 gCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAgg                                                                                          <  2:233574/151‑1 (MQ=255)
                                                  cTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAgg                                                                                          >  2:946708/1‑150 (MQ=255)
                                                    tGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGCTGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaa                                                                                        >  1:33925/1‑150 (MQ=255)
                                                     ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaa                                                                                       <  2:1369014/150‑1 (MQ=255)
                                                     ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaa                                                                                       <  2:256744/150‑1 (MQ=255)
                                                     ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATGCTGTAAGCATGGGCCGTTGTCTGCTCAACCATTCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaaa                                                                                      <  1:1140447/151‑1 (MQ=255)
                                                         gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATAtt                                                                                  >  1:1365181/1‑151 (MQ=255)
                                                          aaCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                                                                                 <  1:1082117/151‑1 (MQ=255)
                                                          aaCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGACGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                                                                                 >  2:355668/1‑151 (MQ=255)
                                                           aCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAata                                                                                    >  2:1218608/1‑147 (MQ=255)
                                                           aCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCa                                                                                <  2:709422/151‑1 (MQ=255)
                                                              ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                                                               <  1:537925/149‑1 (MQ=255)
                                                               tGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcg                                                                             <  1:1159289/150‑1 (MQ=255)
                                                                  cGTTAGCTCGAAAGCGCGCGCCAAAATTCGTGAGTTGCTGAAACACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGa                                                                          <  2:41088/150‑1 (MQ=255)
                                                                    ttAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTAAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCt                                                                       >  1:727832/1‑151 (MQ=255)
                                                                     tAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTg                                                                      >  1:762288/1‑151 (MQ=255)
                                                                      aGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTgg                                                                     <  1:1388600/151‑1 (MQ=255)
                                                                        cTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTg                                                                      <  1:70492/148‑1 (MQ=255)
                                                                         tCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                                                                  >  1:82924/1‑151 (MQ=255)
                                                                          ctaaaGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCg                                                                 >  2:1366473/3‑151 (MQ=255)
                                                                          cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGAt                                                                   <  1:885956/149‑1 (MQ=255)
                                                                          cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGAt                                                                   <  1:27027/149‑1 (MQ=255)
                                                                          cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                                                                  <  2:204503/150‑1 (MQ=255)
                                                                          cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                                                                  <  1:1114244/150‑1 (MQ=255)
                                                                           gAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGc                                                                >  1:634706/1‑151 (MQ=255)
                                                                                cgcgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGa                                                            >  1:40797/1‑150 (MQ=255)
                                                                                   gcgcCAAAATTCGTCAGTTGCTGAAAAACCACAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCt                                                        <  2:778964/151‑1 (MQ=255)
                                                                                     gcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTgg                                                      >  1:581073/1‑151 (MQ=255)
                                                                                       cAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGc                                                     >  2:266267/1‑150 (MQ=255)
                                                                                       cAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCa                                                    <  2:375973/151‑1 (MQ=255)
                                                                                               gTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGa                                            <  2:140576/151‑1 (MQ=255)
                                                                                                tCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGAc                                           >  1:1300008/1‑151 (MQ=255)
                                                                                                   gTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGAt                                        <  2:953578/151‑1 (MQ=255)
                                                                                                   gTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGAt                                        <  2:846905/151‑1 (MQ=255)
                                                                                                   gTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGAt                                        >  1:284618/1‑151 (MQ=255)
                                                                                                   gTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCCCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGAt                                        <  2:783367/151‑1 (MQ=255)
                                                                                                    ttGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATc                                       >  2:596496/1‑151 (MQ=255)
                                                                                                     tGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATc                                       <  2:388784/150‑1 (MQ=255)
                                                                                                      gCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATct                                      <  2:953201/150‑1 (MQ=255)
                                                                                                        tGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgct                                   >  1:1006643/1‑151 (MQ=255)
                                                                                                        tGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgct                                   >  1:1354332/1‑151 (MQ=255)
                                                                                                         gAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgctg                                  >  1:537080/1‑151 (MQ=255)
                                                                                                         gAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgc                                    <  1:247200/149‑1 (MQ=255)
                                                                                                          aaaaaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGTCAACGCTTGACGATCTGCTgg                                 >  1:335416/1‑151 (MQ=255)
                                                                                                            aaaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGc                                >  2:503578/1‑150 (MQ=255)
                                                                                                                 tCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAAt                          >  2:1059501/1‑151 (MQ=255)
                                                                                                                  cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGaaa                           <  2:22885/149‑1 (MQ=255)
                                                                                                                  cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGaaa                           <  2:134731/149‑1 (MQ=255)
                                                                                                                  cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATc                         >  2:1269811/1‑151 (MQ=255)
                                                                                                                   aaGCGTTATGATTCTGTAAGCCTGGGCCGTTGTCTGCGCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATc                         <  1:372443/150‑1 (MQ=255)
                                                                                                                   aaGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCg                        >  2:847184/1‑151 (MQ=255)
                                                                                                                    aGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCgg                       <  2:262085/151‑1 (MQ=255)
                                                                                                                       gTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACt                    <  2:40799/151‑1 (MQ=255)
                                                                                                                           tgatTCTGTAAGCCTGGGCCGTTGTCGGCTCAACCATGCTTTGGGTGGTAGCAGTAAGCTGAAAGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTgg                 <  1:373600/150‑1 (MQ=255)
                                                                                                                            gatTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTa               <  1:106125/151‑1 (MQ=255)
                                                                                                                             atTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTaa              <  2:1420141/151‑1 (MQ=255)
                                                                                                                             atTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTa               >  1:996807/1‑150 (MQ=255)
                                                                                                                             atTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTa               >  2:1346440/1‑150 (MQ=255)
                                                                                                                                cTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGc           >  2:1158572/1‑151 (MQ=255)
                                                                                                                                 tGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGc           >  2:867447/1‑150 (MQ=255)
                                                                                                                                    aaGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAAt        <  1:39971/150‑1 (MQ=255)
                                                                                                                                    aaGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAAt        <  2:503722/150‑1 (MQ=255)
                                                                                                                                       ccTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGc    >  1:1262890/1‑151 (MQ=255)
                                                                                                                                          gggCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCgt  <  1:286348/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGT  >  NC_000913/3823716‑3824003

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