Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,377,768 C→T 6.9% intergenic (‑348/+47) sspA ← / ← rpsI stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,377,7680CT 6.9% 311.7 / 6.9 101intergenic (‑348/+47)sspA/rpsIstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9
Reads supporting (aligned to +/- strand):  ref base C (54/40);  new base T (5/2);  total (59/42)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 107 total aligned reads displayed.

GGGGAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGC  >  NC_000913/3377619‑3377914
                                                                                                                                                      |                                                                                                                                                  
ggggAAAAAGAGTTGCGCTTTAGTGTGAAATAGTTGACTTAGTCCCTAATTGGCGATGTGGTTTTTGTGTTACCTGTCTGTCAGGTGGCAGCAAAAAGCTACTTTCCTGTTTTTACGGTGATTCTGATTTTAGCTAT‑AAAAAAACacgcc                                                                                                                                                    >  2:329917/1‑150 (MQ=255)
 gggAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTTTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACccgccg                                                                                                                                                   <  1:95894/150‑1 (MQ=255)
 gggAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACccgccg                                                                                                                                                   <  1:732254/150‑1 (MQ=255)
  ggAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGa                                                                                                                                                  >  1:1059627/1‑150 (MQ=255)
  ggAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCATTATTTGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGa                                                                                                                                                  >  2:1085081/1‑150 (MQ=255)
   gAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGaa                                                                                                                                                 >  2:1270209/1‑150 (MQ=255)
     aaaaGAGATGCGCTTTAGTCTGAAATATTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGc                                                                                                                                               <  1:861053/150‑1 (MQ=255)
        agagATGCGCTTTAGTCTGAAATAGTTGTCTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCggg                                                                                                                                            <  1:12656/150‑1 (MQ=255)
        agagATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCggg                                                                                                                                            >  2:1556505/1‑150 (MQ=255)
        agagATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCggg                                                                                                                                            <  1:1325746/150‑1 (MQ=255)
         gagaTGCGCTTTAGTCTTAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGt                                                                                                                                           <  2:237891/150‑1 (MQ=255)
         gagaTGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGt                                                                                                                                           <  2:855843/150‑1 (MQ=255)
         gagaTGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGt                                                                                                                                           >  2:1204991/1‑150 (MQ=255)
             tGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGttttt                                                                                                                                       >  1:854433/1‑150 (MQ=255)
             tGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGttttt                                                                                                                                       <  2:1228902/150‑1 (MQ=255)
               cgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCTTCGGCGGGttttttt                                                                                                                                     <  2:56436/150‑2 (MQ=255)
               cgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTc                                                                                                                                     <  1:1621217/150‑1 (MQ=255)
               cgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATAGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTc                                                                                                                                     <  1:1081929/150‑1 (MQ=255)
                gcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCAGACGAAGCGGGTTTTTTCg                                                                                                                                    >  1:1099144/1‑150 (MQ=255)
                 cTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGa                                                                                                                                   <  2:470974/150‑1 (MQ=255)
                    tAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCTTCGGCGGGtttttttatagc                                                                                                                                <  1:1354542/150‑7 (MQ=255)
                      gTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAAt                                                                                                                               >  1:585437/1‑150 (MQ=255)
                          gAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                          >  1:939417/1‑150 (MQ=255)
                             aTAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTAATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCt                                                                                                                       <  1:827045/150‑1 (MQ=255)
                                  tGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCgg                                                                                                                  <  2:729114/150‑1 (MQ=255)
                                   gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGa                                                                                                                 >  2:416801/1‑150 (MQ=255)
                                       tattCCCTTCTTGGCGATGTGGTGTTGGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGATTTTACGCTGATTCAGTTTTTAGATAT‑AAAAAAACCCGCCGATGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAg                                                                                                             <  1:236484/147‑1 (MQ=255)
                                       tAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAg                                                                                                             >  1:125287/1‑150 (MQ=255)
                                           cccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑GAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                         <  1:1486911/150‑1 (MQ=255)
                                            ccTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGcc                                                                                                        <  1:1327645/150‑1 (MQ=255)
                                             cTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGcc                                                                                                        <  2:585437/150‑1 (MQ=255)
                                               tATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAAt                                                                                                     <  2:139014/150‑1 (MQ=255)
                                               tATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCAGGAGCAGAAGCCAAt                                                                                                     >  1:538329/1‑150 (MQ=255)
                                                aTTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAAtt                                                                                                    >  1:1140685/1‑150 (MQ=255)
                                                 ttGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTa                                                                                                   <  1:416801/150‑1 (MQ=255)
                                                     cGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGt                                                                                               >  2:1029215/1‑150 (MQ=255)
                                                      gATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGtt                                                                                              <  2:854433/150‑1 (MQ=255)
                                                          tgGTTTTTGTTTTACCTGTCTGTCAGGTTGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTATCGTTTGga                                                                                          <  1:1085081/150‑1 (MQ=255)
                                                          tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                          <  2:1375707/150‑1 (MQ=255)
                                                          tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                          <  1:248622/150‑1 (MQ=255)
                                                           ggTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgag                                                                                         >  2:472094/1‑150 (MQ=255)
                                                           ggTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgag                                                                                         >  1:1419607/1‑150 (MQ=255)
                                                            gTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAATCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgaga                                                                                        <  1:468015/150‑1 (MQ=255)
                                                            gTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGgaga                                                                                        <  2:4839/150‑1 (MQ=255)
                                                            gTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGGTTTTACGCTGATTCAGATTTTTGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCTGAGCAGAAGCCAATTAACGTTTgggta                                                                                        >  1:674736/1‑147 (MQ=255)
                                                             tttttGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTTACGTTTGGAGaa                                                                                       >  2:588844/1‑150 (MQ=255)
                                                              ttttGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTATTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAAc                                                                                      <  1:1052424/150‑1 (MQ=255)
                                                               tttGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACt                                                                                     <  2:1206344/150‑1 (MQ=255)
                                                                  gTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCg                                                                                  >  1:1639066/1‑150 (MQ=255)
                                                                  gTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGACGAAGCAGAAGCAAATTAACGTTTGGAGAACTgag                                                                                  >  1:331960/1‑148 (MQ=255)
                                                                   ttttACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCgg                                                                                 >  2:1300157/1‑150 (MQ=255)
                                                                      tACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacg                                                                              <  2:6740/150‑1 (MQ=255)
                                                                       aCCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacga                                                                             >  1:826171/1‑150 (MQ=255)
                                                                        ccTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgac                                                                            >  2:108943/1‑150 (MQ=255)
                                                                           gtctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAAGTGCGgacgacgac                                                                         >  2:1523176/1‑150 (MQ=255)
                                                                              tgtcAGGTGGCAGCAAAAAGCAACTTTCCAGGTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                      <  1:1067727/150‑1 (MQ=255)
                                                                                  aGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCttt                                                                  >  2:1266260/1‑150 (MQ=255)
                                                                                   ggTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTa                                                                 >  1:736478/1‑150 (MQ=255)
                                                                                   ggTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTa                                                                 >  1:433111/1‑150 (MQ=255)
                                                                                      ggCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGc                                                              >  1:1408517/1‑150 (MQ=255)
                                                                                       gcagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACgaa                                                             >  1:149842/1‑148 (MQ=255)
                                                                                       gcagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCa                                                             <  2:1059627/150‑1 (MQ=255)
                                                                                           caAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCTGAGCAGAAGCCAATTAACGTTTGGAGAACTGCTGACGACGACGTGCTTTACGCAGAcc                                                         >  1:1106008/1‑150 (MQ=255)
                                                                                              aaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGAc                                                      >  2:1643939/1‑150 (MQ=255)
                                                                                                aGCAACTTTCCAGTATTTACGCTTATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCAGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACtt                                                    >  2:372285/1‑150 (MQ=255)
                                                                                                  cAACATTCCATTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTc                                                  <  1:620346/150‑1 (MQ=255)
                                                                                                   aaCTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCt                                                 >  1:1010016/1‑150 (MQ=255)
                                                                                                     cTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTATTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTa                                               >  2:802238/1‑150 (MQ=255)
                                                                                                             tttttACGCTGATTCAGATTTTAGCTCT‑AAAAAAACCCGCCGACGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTACCGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAc                                       <  1:270583/150‑1 (MQ=255)
                                                                                                              ttttACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAcc                                      <  1:583310/150‑1 (MQ=255)
                                                                                                               tttACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCt                                     >  2:1454131/1‑150 (MQ=255)
                                                                                                                   cGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACg                                 >  1:644797/1‑150 (MQ=255)
                                                                                                                   cGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACg                                 >  2:1441047/1‑150 (MQ=255)
                                                                                                                   cGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACg                                 >  1:739862/1‑150 (MQ=255)
                                                                                                                       gATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCg                             >  2:1548489/1‑150 (MQ=255)
                                                                                                                        aTTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGt                            >  1:696073/1‑150 (MQ=255)
                                                                                                                         ttCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                           >  1:779910/1‑150 (MQ=255)
                                                                                                                          tCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGCCGACGTGCTTTACGCAGCCCGACTTTCTTTCGTTCATCCTGACGAGCGTCa                          >  2:1520666/1‑150 (MQ=255)
                                                                                                                           cAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGTAGCAGAGGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACTTTCAACCTGACGAGCGTCAc                         >  2:1504657/1‑150 (MQ=255)
                                                                                                                           cAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCAc                         <  1:1029215/150‑1 (MQ=255)
                                                                                                                            aGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACg                        >  1:146209/1‑150 (MQ=255)
                                                                                                                            aGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACg                        <  1:1655479/150‑1 (MQ=255)
                                                                                                                             gATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGa                       >  1:651239/1‑150 (MQ=255)
                                                                                                                               tttaatcTAT‑AAAAAAACCCGCCGAAGAGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGAAGAAGCCAATTAACGTTTGGAGAAATGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGt                     <  2:645046/144‑1 (MQ=255)
                                                                                                                                tttAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTa                    >  1:1469713/1‑150 (MQ=255)
                                                                                                                                 ttAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTaa                   <  1:1548489/150‑1 (MQ=255)
                                                                                                                                   aGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACg                 >  2:144370/1‑150 (MQ=255)
                                                                                                                                          gaaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCa           >  2:1285536/2‑150 (MQ=255)
                                                                                                                                          gaaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCa           >  1:1633560/2‑150 (MQ=255)
                                                                                                                                          aaaaacaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCa           <  1:449858/150‑1 (MQ=255)
                                                                                                                                           aaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGTGGACGCCGACGTGCTTTACGCAGACCCACTTTCTTACGTTCAACCTGACGCGCGTCACGAGTAACgccgcccg          >  2:106443/1‑143 (MQ=255)
                                                                                                                                           aaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAg          >  1:241759/1‑150 (MQ=255)
                                                                                                                                             aaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCt        >  1:142279/1‑150 (MQ=255)
                                                                                                                                             aaaaCCCGCCGAAGCGGGTATTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGGACTGCGCACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACCAGCGTCACGAGTAACGAAGCCAGCt        >  1:1551782/1‑150 (MQ=255)
                                                                                                                                              aaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGCCCGACTTTCTTACGTTCAACCTGCCGAGCGTCACGAGTAACGAAGCCAGCtt       >  1:329271/1‑150 (MQ=255)
                                                                                                                                               aaCCCGCCGAAGCGGGGTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGCCCTGCGGACGACGACGTGCGTTAAGCAGACCGACTTTCTTACGTTCATCCTGACGCGCGTCACGAGTAACGAAGCCAGCgtt      >  1:1446475/1‑150 (MQ=255)
                                                                                                                                                 cccGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTAc    >  1:666526/1‑150 (MQ=255)
                                                                                                                                                   cgccgAAGCGGGTTTTTTCGAAAATTTTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGc  <  1:802238/150‑1 (MQ=255)
                                                                                                                                                   cgccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGc  <  2:736478/150‑1 (MQ=255)
                                                                                                                                                   cgccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGAAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGCCGTGCTTTACGCAGACCGACTTTCTTACGTTCACCCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGc  >  1:1276416/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GGGGAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTAT‑AAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGC  >  NC_000913/3377619‑3377914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑