Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,025 G→A R209H (CGT→CAT)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,0250GA96.0% 363.6 / 1.9 99R209H (CGT→CAT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (36/59);  total (36/63)
Fisher's exact test for biased strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

GTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAT  >  NC_000913/3822878‑3823172
                                                                                                                                                   |                                                                                                                                                   
gTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATaa                                                                                                                                                  <  2:1631984/151‑1 (MQ=255)
 tGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATaa                                                                                                                                                  >  3:909452/1‑150 (MQ=255)
 tGAAACTCTCGAAATGTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATaa                                                                                                                                                  <  3:1089909/150‑1 (MQ=255)
  gAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATaaa                                                                                                                                                 <  3:862717/150‑1 (MQ=255)
    aaCTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAagag                                                                                                                                              <  3:1371188/151‑1 (MQ=255)
     aCTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAag                                                                                                                                                >  3:870761/1‑148 (MQ=255)
        ctcGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatg                                                                                                                                           >  3:1617173/1‑150 (MQ=255)
          cGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatgat                                                                                                                                         <  3:1613082/150‑1 (MQ=255)
          cGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATc                                                                                                                                        >  3:1044474/1‑151 (MQ=255)
          cGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATc                                                                                                                                        >  3:1151670/1‑151 (MQ=255)
           gAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATcc                                                                                                                                       >  2:62490/1‑151 (MQ=255)
            aaaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATcc                                                                                                                                       >  2:654865/1‑150 (MQ=255)
             aaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATcc                                                                                                                                       >  2:564746/1‑149 (MQ=255)
             aaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCag                                                                                                                                     <  2:1636551/151‑1 (MQ=255)
              aTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCaga                                                                                                                                    <  2:532456/151‑1 (MQ=255)
                ttATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCagaa                                                                                                                                   <  3:1036261/150‑1 (MQ=255)
                  ataGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCaga                                                                                                                                    >  2:1568787/1‑147 (MQ=255)
                  ataGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAt                                                                                                                                <  3:1135091/151‑1 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttc                                                                                                                              <  2:360468/151‑1 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttc                                                                                                                              >  2:741796/1‑151 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAtt                                                                                                                               <  3:989563/150‑1 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAtt                                                                                                                               <  3:511669/150‑1 (MQ=255)
                     gCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttct                                                                                                                             <  3:850140/151‑1 (MQ=255)
                      cccGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttc                                                                                                                              <  3:840833/149‑1 (MQ=255)
                        cGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttctttc                                                                                                                          <  3:1530775/151‑1 (MQ=255)
                               gcACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATc                                                                                                                   >  2:1609512/1‑151 (MQ=255)
                                cACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCg                                                                                                                  <  2:1475601/151‑1 (MQ=255)
                                  ccGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGaa                                                                                                                <  2:625647/151‑1 (MQ=255)
                                   cGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAg                                                                                                               <  3:587177/151‑1 (MQ=255)
                                     tttAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAggg                                                                                                             >  3:183108/1‑151 (MQ=255)
                                       tAGGTATCCACCGCATTAAAACCGAACTCGACGAGCTGGGTTTTGAGGCGCTGTGACCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAGGATGCGTTCTGAAACCGAAGGGc                                                                                                            <  3:1078221/150‑1 (MQ=255)
                                       tAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGc                                                                                                            <  2:1418207/150‑1 (MQ=255)
                                        aGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGt                                                                                                          >  3:197627/1‑151 (MQ=255)
                                            aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTg                                                                                                       <  2:1088388/150‑1 (MQ=255)
                                              ccaccaCATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGc                                                                                                      >  2:387060/1‑149 (MQ=255)
                                              ccaccaCATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAg                                                                                                    <  2:1077711/151‑1 (MQ=255)
                                                 ccacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGaa                                                                                                 <  3:42674/151‑1 (MQ=255)
                                                       tAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCggg                                                                                            >  3:643135/1‑150 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgc                                                                               >  2:692377/1‑151 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgc                                                                               <  3:588637/151‑1 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcg                                                                                <  3:1136516/150‑1 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcg                                                                                <  2:328175/150‑1 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcg                                                                                <  2:124212/150‑1 (MQ=255)
                                                                    gAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgcg                                                                              >  3:1464059/1‑151 (MQ=255)
                                                                       gagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCa                                                                           <  2:745048/151‑1 (MQ=255)
                                                                       gagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCa                                                                           <  2:1313368/151‑1 (MQ=255)
                                                                        agCTGGTTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAGTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCa                                                                           <  3:1048449/150‑1 (MQ=255)
                                                                        agCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTc                                                                            >  3:55321/1‑149 (MQ=255)
                                                                         gCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGt                                                                         <  2:974847/151‑1 (MQ=255)
                                                                            gggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTg                                                                        <  2:1082145/149‑1 (MQ=255)
                                                                            gggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTg                                                                        <  2:404648/149‑1 (MQ=255)
                                                                             ggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTgg                                                                       <  2:775046/149‑1 (MQ=255)
                                                                                  tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCga                                                                 <  2:1271577/150‑1 (MQ=255)
                                                                                   gAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTcgcg                                                                  >  2:909810/1‑148 (MQ=255)
                                                                                   gAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCgaga                                                               <  3:185461/151‑1 (MQ=255)
                                                                                    aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGaa                                                              <  2:414071/151‑1 (MQ=255)
                                                                                      gcgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGc                                                            <  2:1141453/151‑1 (MQ=255)
                                                                                      gcgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGc                                                            <  2:1171681/151‑1 (MQ=255)
                                                                                       cgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGc                                                            >  2:726989/1‑150 (MQ=255)
                                                                                        gcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCAt                                                          >  2:991383/1‑151 (MQ=255)
                                                                                         cTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATc                                                         >  3:926436/1‑151 (MQ=255)
                                                                                             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCttt                                                      <  3:783786/150‑1 (MQ=255)
                                                                                             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTa                                                     <  3:755947/151‑1 (MQ=255)
                                                                                               cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAt                                                    <  2:1122038/150‑1 (MQ=255)
                                                                                                 cAACCGTTATCACGTAATCAAAGAAGTGGTGAAAGCCGCGCACGGCAACCATAAAGAGATGAGCCAGAAGATTCTTTCTGAAATCGAAGGGCGTGTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCAGCTTGATTCg                                                 <  3:868043/151‑1 (MQ=255)
                                                                                                   aCCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAt                                               <  2:411275/151‑1 (MQ=255)
                                                                                                         aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTg                                         >  2:651260/1‑151 (MQ=255)
                                                                                                           cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCa                                       <  3:1413905/151‑1 (MQ=255)
                                                                                                            gcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCa                                       >  2:385907/1‑150 (MQ=255)
                                                                                                             cgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaa                                     <  3:796536/151‑1 (MQ=255)
                                                                                                             cgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaa                                      >  2:947622/1‑150 (MQ=255)
                                                                                                              gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaa                                     >  2:1162874/1‑150 (MQ=255)
                                                                                                               tAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                   >  3:1288974/1‑151 (MQ=255)
                                                                                                                aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                   <  3:949501/150‑1 (MQ=255)
                                                                                                                aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                   <  3:422312/150‑1 (MQ=255)
                                                                                                                aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATg                                  >  3:1328761/1‑151 (MQ=255)
                                                                                                                   cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGt                                <  2:1449207/150‑1 (MQ=255)
                                                                                                                   cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTg                               <  2:272989/151‑1 (MQ=255)
                                                                                                                    aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGc                              <  2:813357/151‑1 (MQ=255)
                                                                                                                        aagTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                          <  2:1454962/151‑1 (MQ=255)
                                                                                                                            ggtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAaga                       <  2:1582799/150‑1 (MQ=255)
                                                                                                                             gtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     <  2:1446418/151‑1 (MQ=255)
                                                                                                                             gtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     >  2:76681/1‑151 (MQ=255)
                                                                                                                             gtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     <  2:260801/151‑1 (MQ=255)
                                                                                                                               gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                    <  2:946618/150‑1 (MQ=255)
                                                                                                                               gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                    <  2:782683/150‑1 (MQ=255)
                                                                                                                               gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                    <  2:1393006/150‑1 (MQ=255)
                                                                                                                                 aaGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCg                 <  2:331556/151‑1 (MQ=255)
                                                                                                                                     cgcgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGtttt             >  4:9739/1‑151 (MQ=255)
                                                                                                                                       cgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc            >  3:656755/1‑150 (MQ=255)
                                                                                                                                       cgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCa           <  2:934122/151‑1 (MQ=255)
                                                                                                                                        gcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCAc          >  3:1474299/1‑151 (MQ=255)
                                                                                                                                          gcgGCAACCATAAAGAGATGATCCAGAAGATTCCTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCAcg         <  3:1277109/150‑2 (MQ=255)
                                                                                                                                            ggCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCg       >  3:605904/1‑150 (MQ=255)
                                                                                                                                             gCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGAt     <  3:881482/151‑1 (MQ=255)
                                                                                                                                             gCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGAt     >  2:249112/1‑151 (MQ=255)
                                                                                                                                               aaCCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATc    >  2:1603200/1‑150 (MQ=255)
                                                                                                                                                aCCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt  <  3:867669/151‑1 (MQ=255)
                                                                                                                                                aCCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt  >  4:5360/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAT  >  NC_000913/3822878‑3823172

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