Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,987 G→T G530C (GGT→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,9870GT98.7% 287.2 / ‑4.1 76G530C (GGT→TGT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (34/41);  total (34/42)
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.42e-01

GTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTAT  >  NC_000913/3823846‑3824129
                                                                                                                                             |                                                                                                                                              
gTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGc                                                                                                                                         <  1:364421/149‑1 (MQ=255)
    gCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGATGATGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAAt                                                                                                                                      >  2:330551/1‑148 (MQ=255)
      ttGGGCCGTCGTCTGCTCAACCATGCTTTGGGAGGTAGCCGTAAGCTGAATGAAATCGCACAGGAAAATATTCAGAGTGAGCTGGATTGCATTAAACTGGAAATGCTAGACGATCTGCTGGCATAAATGGGACTTTGTAACGCAATGa                                                                                                                                    >  2:762873/2‑148 (MQ=255)
       tGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAGTCGGACTTTGTAACGCAATGAGCgt                                                                                                                                >  1:1207610/1‑151 (MQ=255)
       tGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCgt                                                                                                                                >  1:861575/1‑151 (MQ=255)
           ccGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCgtggtg                                                                                                                            >  1:944149/1‑151 (MQ=255)
            cgtcgtCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGtggtgg                                                                                                                           <  2:984862/151‑1 (MQ=255)
            cgtcgtCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGtggtgg                                                                                                                           <  1:683407/151‑1 (MQ=255)
            cgtcgtCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGtggtgg                                                                                                                           <  2:1127251/151‑1 (MQ=255)
             gtcgtcTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTggtggt                                                                                                                          <  2:205675/151‑1 (MQ=255)
               cgtcTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTc                                                                                                                         <  2:255250/150‑1 (MQ=255)
                 tcggCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCCGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTAACGATCTGCTGGCAGAAATCGGACTTTGTAACGGAATGAGCGTGGTGGTcg                                                                                                                        <  2:893888/146‑1 (MQ=255)
                    gCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGaa                                                                                                                    <  2:1191876/150‑1 (MQ=255)
                          ccATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTg                                                                                                             <  2:37794/151‑1 (MQ=255)
                          ccATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTg                                                                                                             <  1:312153/151‑1 (MQ=255)
                                    ggtggtAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTAGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATgggg                                                                                                   >  2:758680/1‑151 (MQ=255)
                                     gtggtAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATg                                                                                                      >  2:469351/1‑147 (MQ=255)
                                       ggtAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACg                                                                                                >  1:448142/1‑151 (MQ=255)
                                       ggtAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAAACTGCAACATGGGGAc                                                                                                 <  1:162382/150‑1 (MQ=255)
                                          aGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCt                                                                                             >  2:1188734/1‑151 (MQ=255)
                                           gCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTc                                                                                            <  1:259398/151‑1 (MQ=255)
                                           gCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTc                                                                                            >  2:276553/1‑151 (MQ=255)
                                           gCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTc                                                                                            >  2:235638/1‑151 (MQ=255)
                                            ccGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTcc                                                                                           >  1:360224/1‑151 (MQ=255)
                                                aaGCGGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTc                                                                                       <  2:994168/151‑1 (MQ=255)
                                                 aGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTcc                                                                                      >  1:209087/1‑151 (MQ=255)
                                                 aGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTcc                                                                                      >  2:331854/1‑151 (MQ=255)
                                                     gaatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCg                                                                                  <  1:1056017/151‑1 (MQ=255)
                                                     gaatgaaATACCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGCCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCATCCATTCCACCg                                                                                  >  1:916487/1‑151 (MQ=255)
                                                      aatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCgg                                                                                 <  1:83502/151‑1 (MQ=255)
                                                        tgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCgg                                                                                 >  2:554442/1‑149 (MQ=255)
                                                         gaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCa                                                                               <  2:846421/150‑1 (MQ=255)
                                                          aaaTCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAAc                                                                             >  1:701330/1‑151 (MQ=255)
                                                          aaaTCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAAc                                                                             >  2:12217/1‑151 (MQ=255)
                                                               ccGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCaa                                                                         >  2:550835/1‑150 (MQ=255)
                                                                   aGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCa                                                                    <  1:1018004/151‑1 (MQ=255)
                                                                    ggAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCAc                                                                   >  2:849220/1‑151 (MQ=255)
                                                                     gAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACg                                                                  >  1:451010/1‑151 (MQ=255)
                                                                      aaaaTATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACgg                                                                 <  1:121664/151‑1 (MQ=255)
                                                                      aaaaTATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACgg                                                                 >  2:1121743/1‑151 (MQ=255)
                                                                        aaTATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGAc                                                               <  2:1036427/151‑1 (MQ=255)
                                                                                  gcgAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTCACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGc                                                         <  2:81369/147‑1 (MQ=255)
                                                                                  gcgAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCAt                                                     <  1:509522/151‑1 (MQ=255)
                                                                                  gcgAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCAt                                                     <  1:31104/151‑1 (MQ=255)
                                                                                    gAGCTGGAACGCATGAAGCTGTCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGAGTGGGGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGACACGGACATCTGCCCAtt                                                    >  2:1198301/1‑150 (MQ=255)
                                                                                     aGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTaa                                                  <  2:191015/151‑1 (MQ=255)
                                                                                        tttctCGCCTGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAgg                                               >  1:523864/5‑151 (MQ=255)
                                                                                         ggATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAg                                                <  2:309104/149‑1 (MQ=255)
                                                                                         ggATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAg                                                <  2:1195397/149‑1 (MQ=255)
                                                                                         ggATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGACTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGt                                              <  2:892468/151‑1 (MQ=255)
                                                                                           aTCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGTCCATTAAAGGTGc                                            >  1:314682/1‑151 (MQ=255)
                                                                                               cATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGAt                                        >  1:1075786/1‑151 (MQ=255)
                                                                                               cATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGAt                                        >  2:1005056/1‑151 (MQ=255)
                                                                                                  gAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGc                                     <  2:516675/151‑1 (MQ=255)
                                                                                                  gAAGCTGGCAACGCTTGACGAACTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGc                                     <  1:1211391/151‑1 (MQ=255)
                                                                                                        ggCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTg                               <  2:128314/151‑1 (MQ=255)
                                                                                                         gCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTg                               >  1:776214/1‑150 (MQ=255)
                                                                                                         gCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGa                              <  1:331866/151‑1 (MQ=255)
                                                                                                          cAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGAt                             <  2:1019124/151‑1 (MQ=255)
                                                                                                                ttGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCtt                       >  1:785902/1‑151 (MQ=255)
                                                                                                                 tGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCt                        <  2:1023689/149‑1 (MQ=255)
                                                                                                                   aCGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGc                    <  1:909293/151‑1 (MQ=255)
                                                                                                                    cGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCg                   >  2:259122/1‑151 (MQ=255)
                                                                                                                     gATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCg                   >  2:187762/1‑150 (MQ=255)
                                                                                                                             ggCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAAtgctgc          <  2:445370/151‑1 (MQ=255)
                                                                                                                             ggCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAAtgctgc          >  1:657751/1‑151 (MQ=255)
                                                                                                                             ggCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAAtgctg           <  2:1175950/150‑1 (MQ=255)
                                                                                                                             ggCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAAtgctg           <  1:409888/150‑1 (MQ=255)
                                                                                                                                aGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTgccgc       >  2:751962/1‑151 (MQ=255)
                                                                                                                                aGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTgccgc       <  2:664110/151‑1 (MQ=255)
                                                                                                                                 gAAATCGGACTTTGTAACTCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTgccggc      <  2:851446/151‑3 (MQ=255)
                                                                                                                                 gAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTgccgcc      <  2:863222/151‑1 (MQ=255)
                                                                                                                                  aaaTCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGccc     <  1:291702/151‑1 (MQ=255)
                                                                                                                                   aaTCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCt    <  2:450989/151‑1 (MQ=255)
                                                                                                                                     tCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTAt  >  2:644317/1‑151 (MQ=255)
                                                                                                                                      cGGACTTGGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTAt  <  2:838468/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTAT  >  NC_000913/3823846‑3824129

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