Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,811 G→A R471H (CGT→CAT)  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,8110GA99.1% 400.6 / NA 105R471H (CGT→CAT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (51/53);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCT  >  NC_000913/3823663‑3823958
                                                                                                                                                    |                                                                                                                                                   
gTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATc                                                                                                                                                   >  2:1425000/1‑151 (MQ=255)
gTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATc                                                                                                                                                   >  2:310817/1‑151 (MQ=255)
gTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATc                                                                                                                                                   <  1:131362/151‑1 (MQ=255)
 tGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATc                                                                                                                                                   >  1:401757/1‑150 (MQ=255)
  gggCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATtcatca                                                                                                                                                  <  1:444336/150‑1 (MQ=255)
  gggCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATc                                                                                                                                                   >  1:1510601/1‑149 (MQ=255)
  gggCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAg                                                                                                                                                 <  1:95144/151‑1 (MQ=255)
    gcgcACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGtt                                                                                                                                               >  1:225806/1‑151 (MQ=255)
         cgcgTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGa                                                                                                                                          <  1:157045/151‑1 (MQ=255)
          gcgTTGACCACCAGCCTTCCCCGCTGTCGCAGCCGCTTACCAGCGCTCAAACCGTTGAAATCATTACCGCTCAGGGCGCTCGCCCGAATGCCGCTTGGCTGAGCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTg                                                                                                                                           <  2:303682/149‑1 (MQ=255)
            gTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGaa                                                                                                                                         <  1:1009139/149‑1 (MQ=255)
            gTTGACCACCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCCCCGGGCGCTCGCCCGAATGCCCCTTCGCTGAACTTTGTCGTTAGCTCGAAAGCGCACACCAAAATTCATCAGATGCTGaaa                                                                                                                                        >  2:1165760/1‑150 (MQ=255)
              tGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAAc                                                                                                                                     >  1:188951/1‑151 (MQ=255)
                     cAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGc                                                                                                                              >  2:55398/1‑151 (MQ=255)
                        ccTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGtg                                                                                                                           >  2:854294/1‑151 (MQ=255)
                          ttACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGtgat                                                                                                                         >  2:1466528/1‑151 (MQ=255)
                            aCCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGtgatg                                                                                                                        <  2:103300/150‑1 (MQ=255)
                              ccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGAtt                                                                                                                     >  1:712542/1‑151 (MQ=255)
                               cGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTc                                                                                                                    >  2:275192/1‑151 (MQ=255)
                                   gTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                                                                                                                  <  1:93646/149‑1 (MQ=255)
                                     cGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTa                                                                                                                >  2:913792/1‑149 (MQ=255)
                                      gCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGc                                                                                                             <  2:558370/151‑1 (MQ=255)
                                        aGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCt                                                                                                           <  1:766730/151‑1 (MQ=255)
                                              ttACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGcc                                                                                                      >  1:643364/1‑150 (MQ=255)
                                                 ccAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgtc                                                                                                   <  2:373139/150‑1 (MQ=255)
                                                  cAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgtcgt                                                                                                 >  2:600990/1‑151 (MQ=255)
                                                  cAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgtcgt                                                                                                 <  1:387307/151‑1 (MQ=255)
                                                   aGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCgtcgtc                                                                                                <  2:1508765/151‑1 (MQ=255)
                                                   aGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCgtcgtc                                                                                                <  1:734825/151‑1 (MQ=255)
                                                    gCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCt                                                                                               <  1:753255/151‑1 (MQ=255)
                                                    gCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCt                                                                                               >  2:925331/1‑151 (MQ=255)
                                                          aaaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCaa                                                                                         <  2:1073995/151‑1 (MQ=255)
                                                          aaaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCa                                                                                          >  1:957853/1‑150 (MQ=255)
                                                             ccGTTGAAATCATTACCGCTCCGGGTGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAAcc                                                                                       >  1:550443/1‑150 (MQ=255)
                                                             ccGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAAcc                                                                                       >  1:420581/1‑150 (MQ=255)
                                                                 tGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCt                                                                                  >  2:971604/1‑151 (MQ=255)
                                                                 tGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCt                                                                                  >  1:1100445/1‑151 (MQ=255)
                                                                  gAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCt                                                                                  <  1:1276063/150‑1 (MQ=255)
                                                                   aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCttt                                                                                >  1:166040/1‑151 (MQ=255)
                                                                    aaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTg                                                                               <  1:958626/151‑1 (MQ=255)
                                                                     aTCATTACCGCTCCGGGCGCTCGCCCGAGTGCCGCTTGGCGGAACTTTGTCGTTAGCTCGAAAGCCCGCGCGGAAAGTCGGCGGTTGCTGAAAACCCTCAAGCGTGTTGATTCTGTAAGCCGGGGCCGTCGTCTGCTCAGCCATGCttaa                                                                               <  2:919636/150‑3 (MQ=255)
                                                                     aTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTgg                                                                              >  2:965466/1‑151 (MQ=255)
                                                                       cATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTggg                                                                             <  2:377268/150‑1 (MQ=255)
                                                                       cATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTggg                                                                             <  1:140681/150‑1 (MQ=255)
                                                                        aTTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtg                                                                           >  2:253979/1‑151 (MQ=255)
                                                                        aTTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtg                                                                           >  1:1085749/1‑151 (MQ=255)
                                                                        aTTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGggtg                                                                           >  1:1260877/1‑151 (MQ=255)
                                                                           aCCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTa                                                                        <  2:212980/151‑1 (MQ=255)
                                                                            ccGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAg                                                                       <  2:702758/151‑1 (MQ=255)
                                                                                 ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                   >  1:1064439/1‑150 (MQ=255)
                                                                                  cGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTaa                                                                 <  1:212511/151‑1 (MQ=255)
                                                                                   gggCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTa                                                                  >  1:620332/1‑149 (MQ=255)
                                                                                   gggCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAg                                                                <  2:604590/151‑1 (MQ=255)
                                                                                    ggCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGc                                                               <  1:770571/151‑1 (MQ=255)
                                                                                            cccGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatgaa                                                       <  1:259315/151‑1 (MQ=255)
                                                                                              cGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAAt                                                     <  1:173116/151‑1 (MQ=255)
                                                                                                aaTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATcc                                                   >  1:699934/1‑151 (MQ=255)
                                                                                                   gCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                                 <  1:1269706/150‑1 (MQ=255)
                                                                                                   gCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGc                                                <  1:1179586/151‑1 (MQ=255)
                                                                                                    ccGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGc                                                <  1:671898/150‑1 (MQ=255)
                                                                                                    ccGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCa                                               <  2:835035/151‑1 (MQ=255)
                                                                                                    ccGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCa                                               <  1:272125/151‑1 (MQ=255)
                                                                                                     cGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGc                                                >  1:1407083/1‑149 (MQ=255)
                                                                                                       cTTGTCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGCGGTAGCCGTAAGCTGAATGAAATCCCGCAg                                              >  2:136865/1‑149 (MQ=255)
                                                                                                       cTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAgg                                             >  1:833648/1‑150 (MQ=255)
                                                                                                         tGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaa                                           >  2:636060/1‑150 (MQ=255)
                                                                                                         tGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaa                                           <  2:1338038/150‑1 (MQ=255)
                                                                                                           gCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAat                                        <  2:1270912/151‑1 (MQ=255)
                                                                                                              gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAata                                       <  1:897984/149‑1 (MQ=255)
                                                                                                              gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATAtt                                     <  1:1146858/151‑1 (MQ=255)
                                                                                                              gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATAtt                                     >  1:1474981/1‑151 (MQ=255)
                                                                                                               aaCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                                    <  2:1084894/151‑1 (MQ=255)
                                                                                                                   ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcg                                >  2:1246336/1‑151 (MQ=255)
                                                                                                                   ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                  <  1:459523/149‑1 (MQ=255)
                                                                                                                   ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                  <  2:77913/149‑1 (MQ=255)
                                                                                                                    tGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcgc                               >  2:157417/1‑151 (MQ=255)
                                                                                                                       cGTTAGCTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGa                             <  2:53357/150‑1 (MQ=255)
                                                                                                                             cTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGAt                      >  1:1237834/1‑151 (MQ=255)
                                                                                                                             cTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGAt                      >  1:504924/1‑151 (MQ=255)
                                                                                                                             cTCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGAt                      >  2:1084379/1‑151 (MQ=255)
                                                                                                                              tCGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                     >  2:1115927/1‑151 (MQ=255)
                                                                                                                               cGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAGCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                     <  2:1509188/150‑1 (MQ=255)
                                                                                                                               cGAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCg                    <  2:1269790/151‑1 (MQ=255)
                                                                                                                                gAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGc                   <  2:293996/151‑1 (MQ=255)
                                                                                                                                gAAAGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGc                   <  2:225797/151‑1 (MQ=255)
                                                                                                                                  aaGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCAt                 <  1:764077/151‑1 (MQ=255)
                                                                                                                                  aaGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCAt                 <  2:1169334/151‑1 (MQ=255)
                                                                                                                                   aGCGCGCGCCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATg                >  2:42/1‑151 (MQ=255)
                                                                                                                                     cgcgcgcCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGaa              <  1:191157/151‑1 (MQ=255)
                                                                                                                                     cgcgcgcCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGa               >  1:188572/1‑150 (MQ=255)
                                                                                                                                        gcgcCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCt           >  1:367753/1‑151 (MQ=255)
                                                                                                                                         cgcCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTg          <  1:275204/151‑1 (MQ=255)
                                                                                                                                         cgcCAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTg          <  2:484188/151‑1 (MQ=255)
                                                                                                                                           ccAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGc        >  2:778421/1‑151 (MQ=255)
                                                                                                                                            cAAAATTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCa       <  2:127473/151‑1 (MQ=255)
                                                                                                                                             aaaaTTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCaa      <  1:1351616/151‑1 (MQ=255)
                                                                                                                                             aaaaTTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCa       >  1:889589/1‑150 (MQ=255)
                                                                                                                                               aaTTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACg    >  2:1248260/1‑151 (MQ=255)
                                                                                                                                                aTTCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGc   >  2:1459146/1‑151 (MQ=255)
                                                                                                                                                 ttCATCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCt  >  1:1265697/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCT  >  NC_000913/3823663‑3823958

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