Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,799 G→A R467H (CGC→CAC)  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,7990GA98.8% 318.3 / ‑4.1 83R467H (CGC→CAC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (37/45);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

GTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATG  >  NC_000913/3823652‑3823944
                                                                                                                                                   |                                                                                                                                                 
gTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcacgc                                                                                                                                                <  1:951907/151‑1 (MQ=255)
gTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcacgc                                                                                                                                                <  1:341572/151‑1 (MQ=255)
 tCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGcc                                                                                                                                               >  1:24220/1‑151 (MQ=255)
         gCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTc                                                                                                                                        <  2:370911/150‑1 (MQ=255)
         gCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTc                                                                                                                                        <  2:541257/150‑1 (MQ=255)
         gCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTc                                                                                                                                        <  2:894303/150‑1 (MQ=255)
          cgggggCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCTCTTACCAGCGGTCAACCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCaaaa                                                                                                                                           <  2:563902/143‑1 (MQ=255)
          cGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGt                                                                                                                                      <  1:451677/151‑1 (MQ=255)
             gggCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTc                                                                                                                                     >  1:369977/1‑149 (MQ=255)
                  cACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCt                                                                                                                              <  2:1184196/151‑1 (MQ=255)
                   aCGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTg                                                                                                                             >  1:445739/1‑151 (MQ=255)
                    cgcgTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTg                                                                                                                             <  2:1018812/150‑1 (MQ=255)
                       gTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGaaa                                                                                                                          <  2:815855/150‑1 (MQ=255)
                        ttGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGaaaa                                                                                                                         >  2:387542/1‑150 (MQ=255)
                             cGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCa                                                                                                                   <  2:939270/151‑1 (MQ=255)
                               ccAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAg                                                                                                                 <  2:872708/151‑1 (MQ=255)
                                cAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGc                                                                                                                >  1:795756/1‑151 (MQ=255)
                                   ccTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtg                                                                                                             <  1:671882/151‑1 (MQ=255)
                                   ccTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtg                                                                                                             >  2:1084531/1‑151 (MQ=255)
                                       aCCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgatga                                                                                                         <  1:766551/151‑1 (MQ=255)
                                       aCCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgatg                                                                                                          <  1:352613/150‑1 (MQ=255)
                                        cccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTgatgat                                                                                                        <  2:72994/151‑1 (MQ=255)
                                        cccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTgatgat                                                                                                        <  2:976443/151‑1 (MQ=255)
                                        cccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATAATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgatg                                                                                                          >  1:691709/1‑149 (MQ=255)
                                            cTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                                                                                                    >  2:626038/1‑151 (MQ=255)
                                              gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                                                                                                    <  2:971987/149‑1 (MQ=255)
                                              gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                                                                                                    <  2:382461/149‑1 (MQ=255)
                                              gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                                                                                                    <  1:686963/149‑1 (MQ=255)
                                              gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                                                                                                   <  2:1163352/150‑1 (MQ=255)
                                               tCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTa                                                                                                  >  1:198983/1‑150 (MQ=255)
                                                cGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTaa                                                                                                 <  1:372185/150‑1 (MQ=255)
                                                 gCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGc                                                                                               <  1:1038366/151‑1 (MQ=255)
                                                 gCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGc                                                                                               >  1:354719/1‑151 (MQ=255)
                                                     ccGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTgg                                                                                           <  2:223878/151‑1 (MQ=255)
                                                      cGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTgg                                                                                           >  1:986959/1‑150 (MQ=255)
                                                         ttACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcg                                                                                       >  1:1034172/1‑151 (MQ=255)
                                                          tACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgt                                                                                      >  2:232927/1‑151 (MQ=255)
                                                              aGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCgtcgtc                                                                                  <  1:770546/151‑1 (MQ=255)
                                                               gCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCt                                                                                 <  2:638929/151‑1 (MQ=255)
                                                               gCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCt                                                                                 <  2:825830/151‑1 (MQ=255)
                                                                 ggTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGc                                                                               >  2:988314/1‑151 (MQ=255)
                                                                    cAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCa                                                                            <  2:416927/151‑1 (MQ=255)
                                                                      aaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAAc                                                                          <  2:1242502/151‑1 (MQ=255)
                                                                      aaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAAc                                                                          >  1:578476/1‑151 (MQ=255)
                                                                        ccGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCa                                                                        <  1:671880/151‑1 (MQ=255)
                                                                         cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCAt                                                                       <  2:233235/151‑1 (MQ=255)
                                                                          gTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATg                                                                      >  2:28514/1‑151 (MQ=255)
                                                                            tGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCt                                                                    >  1:415124/1‑151 (MQ=255)
                                                                              aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCttt                                                                  >  2:77987/1‑151 (MQ=255)
                                                                                 tCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTggg                                                               >  2:448698/1‑151 (MQ=255)
                                                                                     tACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtgg                                                            <  2:308101/150‑1 (MQ=255)
                                                                                      actgCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtggt                                                           >  2:162160/4‑150 (MQ=255)
                                                                                        cGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGc                                                        >  2:724915/1‑151 (MQ=255)
                                                                                        cGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGc                                                        <  2:242736/151‑1 (MQ=255)
                                                                                           tCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                     >  1:121024/1‑151 (MQ=255)
                                                                                           tCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                     >  2:960725/1‑151 (MQ=255)
                                                                                            ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCg                                                      <  1:579867/149‑1 (MQ=255)
                                                                                            ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                     >  2:1188657/1‑150 (MQ=255)
                                                                                            ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                     >  2:126840/1‑150 (MQ=255)
                                                                                            ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTa                                                    <  2:23581/151‑1 (MQ=255)
                                                                                             tgggCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTa                                                    >  1:188318/2‑150 (MQ=255)
                                                                                                gcgcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCt                                                <  1:1050764/151‑1 (MQ=255)
                                                                                                  gcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtg                                               >  2:668414/1‑150 (MQ=255)
                                                                                                    tCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaat                                            >  2:741891/1‑151 (MQ=255)
                                                                                                     cGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatg                                           <  2:9869/151‑1 (MQ=255)
                                                                                                      gCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatga                                          <  1:438701/151‑1 (MQ=255)
                                                                                                      gCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatg                                           >  2:563858/1‑150 (MQ=255)
                                                                                                        ccGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGaaa                                        <  1:683851/151‑1 (MQ=255)
                                                                                                         cGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAAt                                       >  2:985415/1‑151 (MQ=255)
                                                                                                            aTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATcc                                     <  2:883465/150‑1 (MQ=255)
                                                                                                             tGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                   >  1:515756/1‑151 (MQ=255)
                                                                                                             tGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                   >  1:1057855/1‑151 (MQ=255)
                                                                                                                cGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAg                                <  2:52197/151‑1 (MQ=255)
                                                                                                                  cTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAgg                               >  1:168815/1‑150 (MQ=255)
                                                                                                                   ttGGCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGa                              >  1:438950/1‑150 (MQ=255)
                                                                                                                     ggCTGAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaaa                           <  1:754242/151‑1 (MQ=255)
                                                                                                                         gAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAata                         <  1:991054/149‑1 (MQ=255)
                                                                                                                         gAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAata                         <  1:712477/149‑1 (MQ=255)
                                                                                                                         gAACTTTGTCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATAtt                       >  2:119077/1‑151 (MQ=255)
                                                                                                                          aaCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                      >  2:591159/1‑151 (MQ=255)
                                                                                                                                 tCGTTAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcgc                 >  2:861442/1‑149 (MQ=255)
                                                                                                                                     tAGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTg           >  1:1128170/1‑151 (MQ=255)
                                                                                                                                      aGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCta           <  1:972328/150‑2 (MQ=255)
                                                                                                                                      aGCTCGAAAGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTgg          <  2:325867/151‑1 (MQ=255)
                                                                                                                                              aGCGCACGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATg  <  1:60759/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATG  >  NC_000913/3823652‑3823944

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