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,8610CT97.5% 282.6 / ‑2.0 82R488C (CGT→TGT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (34/45);  total (35/47)
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.70e-01

CAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTGG  >  NC_000913/3823713‑3824008
                                                                                                                                                    |                                                                                                                                                   
cAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgttgt                                                                                                                                                   >  1:813187/1‑151 (MQ=255)
  gCGGTCAAACCGTTGAAATCATTACCGCTCCGGTCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCTTGATGATTCTGTAAGCCTGGGCCGTTGTCt                                                                                                                                                 <  1:406816/151‑1 (MQ=255)
     gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCt                                                                                                                                              >  2:286916/1‑151 (MQ=255)
      tCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTc                                                                                                                                             >  2:389124/1‑151 (MQ=255)
       cAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTc                                                                                                                                             >  2:964866/1‑150 (MQ=255)
        aaaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCa                                                                                                                                            >  2:1034256/1‑150 (MQ=255)
                gAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCtt                                                                                                                                   <  2:268913/151‑1 (MQ=255)
                 aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCttt                                                                                                                                  <  2:1768/151‑1 (MQ=255)
                        tACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGggtgg                                                                                                                            <  1:430604/150‑1 (MQ=255)
                         aCCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGggtggt                                                                                                                           <  1:373363/150‑1 (MQ=255)
                          ccGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAg                                                                                                                         <  1:410011/151‑1 (MQ=255)
                           cGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGc                                                                                                                        <  2:1118974/151‑1 (MQ=255)
                            gCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGcc                                                                                                                       <  2:654468/151‑1 (MQ=255)
                             cTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCg                                                                                                                      >  1:1112863/1‑151 (MQ=255)
                               ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                                                                     <  1:825234/150‑1 (MQ=255)
                               ccGGGCGCTCGCACGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCg                                                                                                                      <  2:33236/149‑1 (MQ=255)
                                 gggCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTaa                                                                                                                   <  1:1103733/150‑1 (MQ=255)
                                   gcgcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCt                                                                                                                <  2:786127/151‑1 (MQ=255)
                                     gcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtg                                                                                                               >  2:151251/1‑150 (MQ=255)
                                     gcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGAGGCAGCCGTAAGCt                                                                                                                >  2:1011568/1‑149 (MQ=255)
                                        cGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatg                                                                                                           >  1:1152495/1‑151 (MQ=255)
                                          cccGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatga                                                                                                          <  2:112747/150‑1 (MQ=255)
                                           ccGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGaaa                                                                                                        <  1:911045/151‑1 (MQ=255)
                                            cGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAAt                                                                                                       <  2:1220535/151‑1 (MQ=255)
                                             gAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATc                                                                                                      <  1:127494/151‑1 (MQ=255)
                                              aaTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATcc                                                                                                     >  1:217815/1‑151 (MQ=255)
                                               aTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATccc                                                                                                    <  2:856679/151‑1 (MQ=255)
                                                tGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                                                                                   >  1:398280/1‑151 (MQ=255)
                                                        ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaaa                                                                                           <  2:1024794/151‑1 (MQ=255)
                                                           tGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAatat                                                                                        >  2:577191/1‑151 (MQ=255)
                                                            gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGAGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATAtt                                                                                       >  1:1087005/1‑151 (MQ=255)
                                                             aaCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                                                                                      <  1:693438/151‑1 (MQ=255)
                                                              aCTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                                                                                      >  1:303743/1‑150 (MQ=255)
                                                                tttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgc                                                                                   >  2:1012919/1‑151 (MQ=255)
                                                                 ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                                                                    <  1:1182199/149‑1 (MQ=255)
                                                                 ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                                                                    <  2:440628/149‑1 (MQ=255)
                                                                 ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                                                                    <  1:648434/149‑1 (MQ=255)
                                                                     cGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAg                                                                              <  2:571099/151‑1 (MQ=255)
                                                                      gTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGc                                                                             <  2:704952/151‑1 (MQ=255)
                                                                         aGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTgg                                                                          >  2:275138/1‑151 (MQ=255)
                                                                          gCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGa                                                                         <  2:534893/151‑1 (MQ=255)
                                                                             cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                                                                       <  1:567324/150‑1 (MQ=255)
                                                                              gAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGc                                                                     >  2:194882/1‑151 (MQ=255)
                                                                                aaGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCAt                                                                   >  1:974296/1‑151 (MQ=255)
                                                                                 aGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATg                                                                  >  2:610620/1‑151 (MQ=255)
                                                                                 aGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATg                                                                  >  2:435702/1‑151 (MQ=255)
                                                                                 aGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGCGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATg                                                                  >  2:283352/1‑151 (MQ=255)
                                                                                  gcgcgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGa                                                                 <  2:105159/151‑1 (MQ=255)
                                                                                   cgcgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGa                                                                 >  1:1031765/1‑150 (MQ=255)
                                                                                    gcgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAg                                                               <  2:990035/151‑1 (MQ=255)
                                                                                    gcgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAg                                                               <  1:600752/151‑1 (MQ=255)
                                                                                     cgcgcCAAAATTCGTCAGTTGCTGAAAAGCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGaa                                                                <  2:87210/149‑1 (MQ=255)
                                                                                     cgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTATCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGaa                                                                <  1:646113/149‑1 (MQ=255)
                                                                                        gcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTgg                                                           <  1:739201/151‑1 (MQ=255)
                                                                                         ccAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGc                                                          >  2:746924/1‑151 (MQ=255)
                                                                                          cAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCa                                                         <  1:539830/151‑1 (MQ=255)
                                                                                          cAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCa                                                         <  2:288875/151‑1 (MQ=255)
                                                                                           aaaaTTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGc                                                          <  2:133797/149‑1 (MQ=255)
                                                                                                    cAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACg                                               <  1:194571/151‑1 (MQ=255)
                                                                                                        tGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATct                                           >  2:677544/1‑151 (MQ=255)
                                                                                                          cTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgc                                         >  1:752518/1‑151 (MQ=255)
                                                                                                           tGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgct                                        >  2:208483/1‑151 (MQ=255)
                                                                                                            gAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgctg                                       >  2:339874/1‑151 (MQ=255)
                                                                                                             aaaaaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgctg                                       <  1:849345/150‑1 (MQ=255)
                                                                                                              aaaaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGc                                     <  1:917841/151‑1 (MQ=255)
                                                                                                               aaaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGc                                     >  1:92385/1‑150 (MQ=255)
                                                                                                                aaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctgctg                                       >  2:674003/1‑147 (MQ=255)
                                                                                                                 aCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGc                                     >  1:1140916/1‑148 (MQ=255)
                                                                                                                    tCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAAt                               >  2:144959/1‑151 (MQ=255)
                                                                                                                       aGCGTGATGATTCTGTAAGCCTGGGCCGTTTTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAAt                               <  1:671346/148‑1 (MQ=255)
                                                                                                                       aGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCgg                            >  1:67918/1‑151 (MQ=255)
                                                                                                                        gCGTGATGATTCTGTAAGCCTGGGCCGTGGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAAGAGTCAGCGCGAGCTGGAGCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGa                           <  1:972073/151‑1 (MQ=255)
                                                                                                                            gatgatTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACtt                        <  2:943090/150‑1 (MQ=255)
                                                                                                                            gatgatTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTg                       <  2:530448/151‑1 (MQ=255)
                                                                                                                             atgatTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTgg                      <  1:893413/151‑1 (MQ=255)
                                                                                                                             atgatTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTgg                      >  2:591436/1‑151 (MQ=255)
                                                                                                                                atTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTaa                   <  1:186902/151‑1 (MQ=255)
                                                                                                                                atTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTa                    >  2:128039/1‑150 (MQ=255)
                                                                                                                                 ttCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTaa                   <  2:398293/150‑1 (MQ=255)
                                                                                                                                     gTAAGCCTGGGCCGTTGTCTGCTCACCCATGCTTTGGGTGGTAGCCGTAAGCTGAATTAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGGAGCTGGCAACGCTTGACGATCTGCTGGGAGAAATCGGACTTGGTAACg                 <  1:1182791/148‑1 (MQ=255)
                                                                                                                                       aaGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATg            >  2:810819/1‑151 (MQ=255)
                                                                                                                                          ccTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGc         >  2:722272/1‑151 (MQ=255)
                                                                                                                                             gggCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCg        <  1:765304/149‑1 (MQ=255)
                                                                                                                                                 cgttgtCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCgtggt    >  2:706603/1‑149 (MQ=255)
                                                                                                                                                 cgttgtCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGtggtgg  <  1:964808/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTGG  >  NC_000913/3823713‑3824008

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