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,8610CT98.8% 289.0 / ‑4.2 81R488C (CGT→TGT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/50);  total (31/50)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTG  >  NC_000913/3823724‑3824007
                                                                                                                                         |                                                                                                                                                  
ccGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAAcc                                                                                                                                        >  2:1190335/1‑150 (MQ=255)
  gTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATg                                                                                                                                     <  2:827368/151‑1 (MQ=255)
  gTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATg                                                                                                                                     >  1:254942/1‑151 (MQ=255)
    tGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGc                                                                                                                                    <  2:1026527/150‑1 (MQ=255)
    tGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGc                                                                                                                                    <  1:717190/150‑1 (MQ=255)
      aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCttt                                                                                                                                 <  2:373336/151‑1 (MQ=255)
      aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCttt                                                                                                                                 <  2:456883/151‑1 (MQ=255)
        aTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTgg                                                                                                                               >  1:766264/1‑151 (MQ=255)
          cATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTggg                                                                                                                              <  1:948888/150‑1 (MQ=255)
              aCCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGggtggt                                                                                                                          <  1:336978/150‑1 (MQ=255)
               ccGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGggtggt                                                                                                                          >  1:968734/1‑149 (MQ=255)
                 gCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGcc                                                                                                                      >  1:937703/1‑151 (MQ=255)
                   tCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                                                                    >  2:1136330/1‑151 (MQ=255)
                    ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                                                                    <  1:593854/150‑1 (MQ=255)
                    ccGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGt                                                                                                                    >  2:747668/1‑150 (MQ=255)
                      gggCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAg                                                                                                                 <  1:867097/151‑1 (MQ=255)
                        gcgcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCt                                                                                                               <  1:963659/151‑1 (MQ=255)
                          gcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTATGCtg                                                                                                              <  1:943080/150‑1 (MQ=255)
                          gcTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtg                                                                                                              >  2:390769/1‑150 (MQ=255)
                             cGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaatg                                                                                                          >  2:204327/1‑151 (MQ=255)
                             cGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaat                                                                                                           <  2:1239383/150‑1 (MQ=255)
                             cGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCtgaa                                                                                                            <  2:324161/149‑1 (MQ=255)
                                  gAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATc                                                                                                     >  2:600457/1‑151 (MQ=255)
                                    aTGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATccc                                                                                                   >  2:373227/1‑151 (MQ=255)
                                      gCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCg                                                                                                  <  1:837469/150‑1 (MQ=255)
                                        cGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGc                                                                                                 >  1:2992/1‑149 (MQ=255)
                                             ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaaa                                                                                          <  2:27854/151‑1 (MQ=255)
                                             ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaa                                                                                           >  2:1155183/1‑150 (MQ=255)
                                             ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaa                                                                                           <  1:1003142/150‑1 (MQ=255)
                                             ggCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGaaa                                                                                           >  1:597890/1‑150 (MQ=255)
                                                 gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAatat                                                                                       <  1:117856/150‑1 (MQ=255)
                                                 gAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAata                                                                                        <  2:364547/149‑1 (MQ=255)
                                                    cTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTc                                                                                     >  2:763699/1‑149 (MQ=255)
                                                      ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcg                                                                                 >  1:922733/1‑151 (MQ=255)
                                                      ttGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAg                                                                                   <  2:326481/149‑1 (MQ=255)
                                                        gTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcgc                                                                                <  2:188337/150‑1 (MQ=255)
                                                        gTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcgc                                                                                <  1:916854/150‑1 (MQ=255)
                                                        gTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGcgcg                                                                               <  2:136416/151‑1 (MQ=255)
                                                        gTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATGCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAgcgc                                                                                <  2:697742/150‑1 (MQ=255)
                                                               gCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGa                                                                        <  2:107068/151‑1 (MQ=255)
                                                               gCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGa                                                                        <  2:1078591/151‑1 (MQ=255)
                                                                  cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                                                                      <  1:923311/150‑1 (MQ=255)
                                                                  cGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATc                                                                      <  1:1256591/150‑1 (MQ=255)
                                                                   gAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGc                                                                    >  2:269867/1‑151 (MQ=255)
                                                                   gAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGc                                                                    <  2:965800/151‑1 (MQ=255)
                                                                        cgcgcgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGaa                                                               <  1:566906/151‑1 (MQ=255)
                                                                            cgcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTg                                                           >  2:219838/1‑151 (MQ=255)
                                                                             gcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTgg                                                          <  1:385363/151‑1 (MQ=255)
                                                                             gcCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTgg                                                          <  2:242279/151‑1 (MQ=255)
                                                                              ccAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGc                                                         >  2:491520/1‑151 (MQ=255)
                                                                                 aaaTTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAAc                                                      >  1:993968/1‑151 (MQ=255)
                                                                                  aaTTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACg                                                     >  2:197940/1‑151 (MQ=255)
                                                                                  aaTTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACg                                                     >  2:280335/1‑151 (MQ=255)
                                                                                  aaTTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACg                                                     >  2:1060473/1‑151 (MQ=255)
                                                                                  aaTTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACg                                                     >  2:550097/1‑151 (MQ=255)
                                                                                    ttCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCt                                                   >  2:1228153/1‑151 (MQ=255)
                                                                                     tCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCt                                                   <  2:1206361/150‑1 (MQ=255)
                                                                                     tCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCt                                                   <  2:301159/150‑1 (MQ=255)
                                                                                      cGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTg                                                 <  2:124606/151‑1 (MQ=255)
                                                                                         cAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGAc                                               <  2:488311/150‑1 (MQ=255)
                                                                                              gCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATct                                          <  1:1258702/150‑1 (MQ=255)
                                                                                               cTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATctg                                         >  1:1032399/1‑150 (MQ=255)
                                                                                                     aaCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAg                                  <  2:1010250/151‑1 (MQ=255)
                                                                                                      aCCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGa                                 <  2:728229/151‑1 (MQ=255)
                                                                                                          cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGaaa                               <  2:997408/149‑1 (MQ=255)
                                                                                                          cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAAt                              <  2:601523/150‑1 (MQ=255)
                                                                                                          cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATc                             <  2:1243688/151‑1 (MQ=255)
                                                                                                          cAAGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATc                             <  1:1050166/151‑1 (MQ=255)
                                                                                                           aaGCGTGATGATTCTGTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAAGCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCg                            <  2:972333/151‑1 (MQ=255)
                                                                                                                          gTAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCa              <  1:1051234/150‑1 (MQ=255)
                                                                                                                           tAAGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAAt            >  1:824398/1‑151 (MQ=255)
                                                                                                                            aaGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAAt            <  1:619104/150‑1 (MQ=255)
                                                                                                                            aaGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATg           <  2:356347/151‑1 (MQ=255)
                                                                                                                             aGCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGa          <  2:688909/151‑1 (MQ=255)
                                                                                                                              gCCTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATg           >  1:724052/1‑149 (MQ=255)
                                                                                                                               ccTGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGc        >  2:198043/1‑151 (MQ=255)
                                                                                                                                 tGGGCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCg       >  1:969700/1‑150 (MQ=255)
                                                                                                                                  gggCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCgtg     <  1:23025/151‑1 (MQ=255)
                                                                                                                                    gCCGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCgtggt   >  2:402829/1‑151 (MQ=255)
                                                                                                                                     ccGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCgtggtg  <  1:679656/151‑1 (MQ=255)
                                                                                                                                     ccGTTGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCgtggtg  <  2:398060/151‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
CCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTG  >  NC_000913/3823724‑3824007

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