Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,774 A→C N459H (AAC→CAC)  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,7740AC98.8% 903.8 / NA 254N459H (AAC→CAC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (134/117);  total (135/119)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

TGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCA  >  NC_000913/3823677‑3823870
                                                                                                 |                                                                                                
tGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCAc                                                                                                <  1:1748622/100‑1 (MQ=255)
  aCCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACtt                                                                                              <  1:3862361/100‑1 (MQ=255)
    cGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTg                                                                                            <  2:2160156/100‑1 (MQ=255)
    cGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTg                                                                                            >  2:4015962/1‑100 (MQ=255)
      ccAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTc                                                                                          <  2:4454471/100‑1 (MQ=255)
         gCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGtt                                                                                       <  1:3581941/100‑1 (MQ=255)
         gCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGtt                                                                                       <  2:5388067/100‑1 (MQ=255)
         gCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGtt                                                                                       >  2:2970004/1‑100 (MQ=255)
          ccTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTa                                                                                      <  1:4013423/100‑1 (MQ=255)
          ccTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTa                                                                                      >  2:3597743/1‑100 (MQ=255)
            ttACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGc                                                                                    <  2:2376254/100‑1 (MQ=255)
            ttACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGc                                                                                    <  2:3037668/100‑1 (MQ=255)
             tACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCt                                                                                   >  1:4515226/1‑100 (MQ=255)
               cccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCGTGGCTGAACTTTGTCGTTAGCTCg                                                                                 >  2:1422526/1‑100 (MQ=255)
                ccGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGa                                                                                >  1:2038068/1‑100 (MQ=255)
                       cGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAgcgcg                                                                         >  1:4148616/1‑100 (MQ=255)
                           gCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGcc                                                                     >  2:3536097/1‑100 (MQ=255)
                           gCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGcc                                                                     >  2:1013004/1‑100 (MQ=255)
                           gCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGcc                                                                     >  2:3107373/1‑100 (MQ=255)
                             cGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCaa                                                                   >  2:1849401/1‑100 (MQ=255)
                              gCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCaaa                                                                  >  2:267865/1‑100 (MQ=255)
                                ttACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAAt                                                                <  2:2174577/100‑1 (MQ=255)
                                ttACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAAt                                                                <  2:1118900/100‑1 (MQ=255)
                                ttACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAAt                                                                <  2:460461/100‑1 (MQ=255)
                                  aCCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTc                                                              >  1:329028/1‑100 (MQ=255)
                                  aCCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCCTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTc                                                              >  2:4514319/1‑100 (MQ=255)
                                   ccAGCGGTGAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAACGCGCGCGCCAAAATTCg                                                             >  2:5275669/1‑100 (MQ=255)
                                   ccAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCg                                                             <  1:3389531/100‑1 (MQ=255)
                                   ccAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCg                                                             <  2:181258/100‑1 (MQ=255)
                                   ccAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCg                                                             >  2:673685/1‑100 (MQ=255)
                                   ccAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCg                                                             >  1:4790356/1‑100 (MQ=255)
                                    cAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGt                                                            >  1:2345372/1‑100 (MQ=255)
                                     aGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTc                                                           <  1:2046322/100‑1 (MQ=255)
                                        ggTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGt                                                        >  1:2603834/1‑100 (MQ=255)
                                         gTCAAACCGTTGAAATCATTACCGCTCGGGGCGCTCGCCCGAGTGCCGCTTGGCTGCACTTTGTCTTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGtt                                                       >  2:2420577/1‑100 (MQ=255)
                                         gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGtt                                                       <  1:2508633/100‑1 (MQ=255)
                                         gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGtt                                                       >  1:4491797/1‑100 (MQ=255)
                                         gTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGtt                                                       >  1:1950617/1‑100 (MQ=255)
                                           cAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGc                                                     <  1:3735410/100‑1 (MQ=255)
                                           cAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGc                                                     >  1:1848323/1‑100 (MQ=255)
                                           cAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGc                                                     >  1:773140/1‑100 (MQ=255)
                                            aaaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCt                                                    >  1:3503040/1‑100 (MQ=255)
                                             aaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTg                                                   <  2:5151410/100‑1 (MQ=255)
                                             aaCCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTg                                                   >  1:5301944/1‑100 (MQ=255)
                                               ccGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaa                                                 >  2:1610656/1‑100 (MQ=255)
                                                cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaa                                                <  2:1220363/100‑1 (MQ=255)
                                                cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaa                                                >  1:4268717/1‑100 (MQ=255)
                                                cGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCACGCGCCAAAATTCGTCAGTTGCTGaaa                                                >  1:4827307/1‑100 (MQ=255)
                                                 gTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAACGCGCGCGCCAAAATTCGTCAGTTGCTGaaaa                                               >  2:1270748/1‑100 (MQ=255)
                                                 gTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaaa                                               <  2:4736173/100‑1 (MQ=255)
                                                 gTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGaaaa                                               <  1:3612792/100‑1 (MQ=255)
                                                    gAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAAcc                                            <  2:3861797/100‑1 (MQ=255)
                                                    gAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAAcc                                            >  1:38190/1‑100 (MQ=255)
                                                     aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCt                                           >  2:1938938/1‑100 (MQ=255)
                                                     aaaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCt                                           >  1:1359339/1‑100 (MQ=255)
                                                      aaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTc                                          >  2:1998084/1‑100 (MQ=255)
                                                      aaTCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTc                                          <  1:2415613/100‑1 (MQ=255)
                                                         cATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTTAGTTGCTGAAAAACCTCAAg                                       <  1:1073334/100‑1 (MQ=255)
                                                         cATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAg                                       <  2:3286513/100‑1 (MQ=255)
                                                           ttACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCCGTTGCTGAAAAACCTCAAGCg                                     >  2:2707664/1‑100 (MQ=255)
                                                           ttACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCg                                     <  2:2808521/100‑1 (MQ=255)
                                                              ccGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtga                                  <  1:611686/100‑1 (MQ=255)
                                                              ccGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATGCGTCAGTTGCTGAAAAACCTCAAGCGtga                                  >  1:1537983/1‑100 (MQ=255)
                                                               gcattcGGGCGAGCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgat                                 >  1:3110619/6‑100 (MQ=255)
                                                               ctcTCCGGGCGCTCGCCCGAATGCCGCGTGGCTGCACTTTGTCGTTAGCTGGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgat                                 <  2:2274102/98‑1 (MQ=255)
                                                               cGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgat                                 <  1:4034700/100‑1 (MQ=255)
                                                               cGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgat                                 <  1:2969778/100‑1 (MQ=255)
                                                                cattcGGGCGAGCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgatg                                >  2:1509293/5‑100 (MQ=255)
                                                                 attcGGGCGAGCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGtgatga                               >  2:3647587/4‑100 (MQ=255)
                                                                  tCCGGGCGCTCGCCCGAATGCCGCTTGGGTGCACGTTGTCGTTGGGGCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCATGCGTgatgat                              <  2:5075518/100‑1 (MQ=255)
                                                                   ccGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGAtt                             >  2:4214109/1‑100 (MQ=255)
                                                                   ccGGGCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGAtt                             >  2:728242/1‑100 (MQ=255)
                                                                      ggCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                          <  2:1094444/100‑1 (MQ=255)
                                                                      ggCGCTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTg                          >  1:3024429/1‑100 (MQ=255)
                                                                       gcgcTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                         <  2:2085105/100‑1 (MQ=255)
                                                                       gcgcTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                         <  2:3235850/100‑1 (MQ=255)
                                                                       gcgcTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGt                         >  1:946849/1‑100 (MQ=255)
                                                                        cgcgcgCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTa                        <  2:2498997/100‑1 (MQ=255)
                                                                        cgcTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTa                        <  1:1016969/100‑1 (MQ=255)
                                                                        cgcTCGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTa                        >  2:601728/1‑100 (MQ=255)
                                                                            cGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGcc                    <  1:5521050/100‑1 (MQ=255)
                                                                            cGCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGcc                    >  2:4427832/1‑100 (MQ=255)
                                                                             gCCCGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCt                   >  1:2856274/1‑100 (MQ=255)
                                                                              cccGAATGGCGCTTGGCTGCACTCTGTCTTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGGTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTg                  <  2:4169076/100‑1 (MQ=255)
                                                                              cccGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTg                  >  2:3751470/1‑100 (MQ=255)
                                                                              cccGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTg                  >  1:4676043/1‑100 (MQ=255)
                                                                              cccGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTg                  >  2:2746870/1‑100 (MQ=255)
                                                                                cGAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTggg                <  1:4940759/100‑1 (MQ=255)
                                                                                 gAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGc               <  1:5005264/100‑1 (MQ=255)
                                                                                 gAATGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGc               <  1:215517/100‑1 (MQ=255)
                                                                                    tGCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgt            >  1:1456953/1‑100 (MQ=255)
                                                                                    ngCCGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgt            >  1:4427024/2‑100 (MQ=255)
                                                                                      ccGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgtcg          <  1:1416475/100‑1 (MQ=255)
                                                                                       cGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgtcgt         >  1:4104070/1‑100 (MQ=255)
                                                                                       cGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgtcgt         >  2:1381704/1‑100 (MQ=255)
                                                                                       cGCTTGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCcgccgt         >  2:3525673/1‑100 (MQ=255)
                                                                                          ttGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTg      >  2:4245935/1‑100 (MQ=255)
                                                                                          ttGGCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTg      <  1:4493016/100‑1 (MQ=255)
                                                                                            ggCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCt    <  1:3205958/100‑1 (MQ=255)
                                                                                             gCTGCACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTc   >  1:1599977/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTCGTCAGTTGCTGAAAAACCTCAAGCGTGATGATTCTGTAAGCCTGGGCCGTCGTCTGCTCA  >  NC_000913/3823677‑3823870

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