Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,913,035 A→G L206P (CTC→CCC)  relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,0350AG98.8% 314.0 / NA 82L206P (CTC→CCC) relA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (40/41);  minor base T (0/1);  total (40/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.93e-01

TGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCA  >  NC_000913/2912888‑2913181
                                                                                                                                                   |                                                                                                                                                  
tgtTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGgggg                                                                                                                                                  >  1:462807/1‑150 (MQ=255)
    tCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTaa                                                                                                                                               <  2:727581/149‑1 (MQ=255)
     cGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAAct                                                                                                                                             <  1:523860/150‑2 (MQ=255)
     cGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACg                                                                                                                                             <  2:462095/150‑1 (MQ=255)
     cGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACg                                                                                                                                             <  2:621952/150‑1 (MQ=255)
      ggACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACgg                                                                                                                                            >  2:986087/1‑150 (MQ=255)
         cgacCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGc                                                                                                                                        <  1:1046048/151‑1 (MQ=255)
             cATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGt                                                                                                                                     >  2:588616/1‑150 (MQ=255)
               tACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAAt                                                                                                                                  >  1:550255/1‑151 (MQ=255)
                   cTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTc                                                                                                                              >  2:836551/1‑151 (MQ=255)
                     tCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTcc                                                                                                                             <  2:193787/150‑1 (MQ=255)
                       cGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGt                                                                                                                          >  1:307092/1‑151 (MQ=255)
                         cTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTc                                                                                                                        >  2:179093/1‑151 (MQ=255)
                         cTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTc                                                                                                                        >  1:11985/1‑151 (MQ=255)
                           ttAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTc                                                                                                                        <  2:656380/149‑1 (MQ=255)
                            tatcGCCTTCAGATCTCATCCCTGTGCGCAGATGACCAACGAACTCTTCGCTGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGCTTGTATTCGGTTGGATGCGGGTAACGGAAGCAGTAATCTTCCAGTTcc                                                                                                                       >  2:165917/4‑149 (MQ=255)
                             aaCGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGACCTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGTGGTAACGGCAGCAGTAATCTTCCAGtcacca                                                                                                                     <  2:1198139/150‑6 (MQ=255)
                                 ccTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCa                                                                                                                <  2:462792/151‑1 (MQ=255)
                                     cAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGtt                                                                                                             <  1:986088/150‑1 (MQ=255)
                                     cAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTg                                                                                                            <  2:71858/151‑1 (MQ=255)
                                     cAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTg                                                                                                            <  1:822692/151‑1 (MQ=255)
                                         tttCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCg                                                                                                        >  1:879545/1‑151 (MQ=255)
                                         tttCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCg                                                                                                        >  1:1041230/1‑151 (MQ=255)
                                            cATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGAtt                                                                                                     >  1:1205289/1‑151 (MQ=255)
                                             aTCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTc                                                                                                    >  2:1059734/1‑151 (MQ=255)
                                              tctcAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTcc                                                                                                   >  1:827247/1‑151 (MQ=255)
                                              tctcAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTcc                                                                                                   >  1:251844/1‑151 (MQ=255)
                                                   gcgcgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGAc                                                                                              <  1:790498/151‑1 (MQ=255)
                                                    cgcgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACg                                                                                             <  1:631080/151‑1 (MQ=255)
                                                     gcgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACgg                                                                                            >  1:1112209/1‑151 (MQ=255)
                                                     gcgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACgg                                                                                            >  2:1187918/1‑151 (MQ=255)
                                                      cgcAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGt                                                                                           <  1:969427/151‑1 (MQ=255)
                                                        cAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTa                                                                                         <  1:465514/151‑1 (MQ=255)
                                                        cAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTa                                                                                         <  2:248807/151‑1 (MQ=255)
                                                        cAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTa                                                                                         <  1:440084/151‑1 (MQ=255)
                                                        cAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTa                                                                                         <  1:352431/151‑1 (MQ=255)
                                                         aGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGtt                                                                                          >  2:427396/1‑149 (MQ=255)
                                                          gATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGc                                                                                       <  1:377815/151‑1 (MQ=255)
                                                            tGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGcc                                                                                      >  2:695335/1‑150 (MQ=255)
                                                              aCCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCGGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGc                                                                                   <  2:1038912/151‑1 (MQ=255)
                                                              aCCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGc                                                                                   >  2:996487/1‑151 (MQ=255)
                                                                cAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCgg                                                                                 >  2:517154/1‑151 (MQ=255)
                                                                 aaCGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCgg                                                                                 >  1:745121/1‑150 (MQ=255)
                                                                  aCGCACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTg                                                                               >  2:302598/1‑151 (MQ=255)
                                                                    gAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCg                                                                             >  1:431862/1‑151 (MQ=255)
                                                                    gAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCg                                                                             >  1:572023/1‑151 (MQ=255)
                                                                    gAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCg                                                                             >  2:424780/1‑151 (MQ=255)
                                                                    aaaCTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCg                                                                             >  1:847734/2‑151 (MQ=255)
                                                                           tCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATg                                                                       <  1:554999/150‑1 (MQ=255)
                                                                            cGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGtt                                                                     <  1:846678/151‑1 (MQ=255)
                                                                             gATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGAt                                                                        <  2:701840/147‑1 (MQ=255)
                                                                             gATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTg                                                                    <  1:663957/151‑1 (MQ=255)
                                                                              aTGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTgg                                                                   <  2:753363/151‑1 (MQ=255)
                                                                                  aGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTaca                                                               <  2:193346/151‑1 (MQ=255)
                                                                                        cgcgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCtt                                                          <  2:1008114/150‑1 (MQ=255)
                                                                                         gcgGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCtttt                                                        >  2:282904/1‑151 (MQ=255)
                                                                                           ggTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTg                                                       <  2:1205303/150‑1 (MQ=255)
                                                                                              cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccgc                                                   <  2:935745/151‑1 (MQ=255)
                                                                                              cGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccg                                                    <  1:314482/150‑1 (MQ=255)
                                                                                               gagaCGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTgccg                                                    <  1:1225433/149‑1 (MQ=255)
                                                                                                         tCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGtt                                        >  2:715955/1‑151 (MQ=255)
                                                                                                               ggcagTGTGGGAATTCGTTTGTGGTAGGTTGGATGTGGCTAACGGCAGCAGTAGTGTTCCCGGTCCCATTTAAGTTATACGATTCCGAGACGGTTAGCCAGCGGAGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTAATCtt                                  >  2:84086/2‑151 (MQ=255)
                                                                                                               agcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCtt                                  <  1:821825/151‑1 (MQ=255)
                                                                                                                gcagTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTc                                 <  1:179080/151‑1 (MQ=255)
                                                                                                                  agTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCg                               <  2:1085276/151‑1 (MQ=255)
                                                                                                                   gTTTGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCgg                              <  2:787784/151‑1 (MQ=255)
                                                                                                                    tttGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCgg                              >  1:52834/1‑150 (MQ=255)
                                                                                                                      tGGCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGgcgc                           >  1:806959/1‑151 (MQ=255)
                                                                                                                        gCAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCAt                         >  1:435173/1‑151 (MQ=255)
                                                                                                                         cAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATc                        >  2:256641/1‑151 (MQ=255)
                                                                                                                         cAATTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATc                        <  1:880359/151‑1 (MQ=255)
                                                                                                                          aaTTCGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCt                       >  1:142427/1‑151 (MQ=255)
                                                                                                                              cGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTa                   <  2:863871/151‑1 (MQ=255)
                                                                                                                              cGTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTa                   <  2:718305/151‑1 (MQ=255)
                                                                                                                               gTTTGTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTAc                  >  1:520456/1‑151 (MQ=255)
                                                                                                                                   gTATTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcg              <  2:827241/151‑1 (MQ=255)
                                                                                                                                     aTTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcgcg            >  1:179944/1‑151 (MQ=255)
                                                                                                                                     aTTCGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTcgcg            >  1:433405/1‑151 (MQ=255)
                                                                                                                                        cGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc           >  2:1163539/1‑149 (MQ=255)
                                                                                                                                        cGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc           >  1:514333/1‑149 (MQ=255)
                                                                                                                                        cGGTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCa          <  1:33895/150‑1 (MQ=255)
                                                                                                                                         ggTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCgcgc           >  1:960720/1‑148 (MQ=255)
                                                                                                                                          gTTGGATGGGGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGAt       <  2:383512/151‑1 (MQ=255)
                                                                                                                                            tggatgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATg      <  2:709394/150‑1 (MQ=255)
                                                                                                                                            tggatgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATg      <  1:583074/150‑1 (MQ=255)
                                                                                                                                            tggatgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATg      <  2:864314/150‑1 (MQ=255)
                                                                                                                                            tggatgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATg      >  1:566727/1‑150 (MQ=255)
                                                                                                                                                tgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCa  >  1:219341/1‑150 (MQ=255)
                                                                                                                                                tgggGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCa  >  1:379437/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTTTCGGACGACCATACACTTCCGCTTTAACGCCTTCAGCTTTCATCTCAGCGCGCAGATGACCAACGAACTCTTCGATGTAGTGTTCGCGGTCGAGACGCCGTTCATGCAGCAGTTTGGCAATTCGTTTGTATTCGGTTGGATGGAGGTAACGGAAGCAGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCA  >  NC_000913/2912888‑2913181

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