Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 5.8% intergenic (+19/+21) metE → / ← ysgA 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3330AT 5.8% 417.7 / 8.4 143intergenic (+19/+21)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (67/67);  new base T (4/5);  total (71/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 146 total aligned reads displayed.

CGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACAC  >  NC_000913/4015184‑4015482
                                                                                                                                                     |                                                                                                                                                     
cGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTaa                                                                                                                                                       >  2:397979/1‑150 (MQ=255)
 ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAt                                                                                                                                                      >  2:759697/1‑150 (MQ=255)
  gCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                     >  1:1758367/1‑150 (MQ=255)
  gCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                     >  1:313687/1‑150 (MQ=255)
  gCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                     >  1:398212/1‑150 (MQ=255)
  gCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                     <  2:1590759/150‑1 (MQ=255)
   cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                    >  2:338870/1‑150 (MQ=255)
    aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                                                                   >  2:348584/1‑150 (MQ=255)
    aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                                                                   >  1:1283518/1‑150 (MQ=255)
     aaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacca                                                                                                                                                  >  1:99007/1‑150 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                                                                             <  1:465894/150‑1 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                                                                             <  1:827204/150‑1 (MQ=255)
              ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCtt                                                                                                                                         >  2:610382/1‑150 (MQ=255)
                ggCCGAGCGCCTGGGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttgga                                                                                                                                       <  1:158452/150‑6 (MQ=255)
                ggCAGAGCGCCTGTGGGGCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCtttt                                                                                                                                       <  1:800199/150‑1 (MQ=255)
                  cAGAGCGCCTGGGGGTCAACCCGGACTGGGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGACCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtct                                                                                                                                     <  2:1741775/150‑1 (MQ=255)
                      gcgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCAtt                                                                                                                                 <  1:967186/150‑1 (MQ=255)
                      gcgcCTGTGGGGCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAATCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCAtt                                                                                                                                 <  2:418811/150‑1 (MQ=255)
                      gcgcCTGTGGGGCAACCCGGACTGTGGCCTGAAAACGCGCGGCGGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCAtt                                                                                                                                 <  2:1129555/150‑1 (MQ=255)
                      gcgcCGGTGGGTCAACCCGGACTGTGGCCTGAAATCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATTGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCAtt                                                                                                                                 <  2:711360/150‑1 (MQ=255)
                       cgcCGGTGGGGCACCCCGGACTGTGGCCTGAAAACGCGCGGCTGCCCAGAACCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGCAATACCACCCGGTCTTTTCTCATTa                                                                                                                                <  1:1097419/150‑1 (MQ=255)
                        gcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACAACCCGGTCTTGTCTCATTAc                                                                                                                               >  2:1366333/1‑150 (MQ=255)
                         ccTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGCACTTGCGTCGCGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACa                                                                                                                              >  2:1450094/1‑150 (MQ=255)
                           tgtgGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGc                                                                                                                            >  2:1595404/1‑150 (MQ=255)
                            gtgGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCg                                                                                                                           >  2:959397/1‑150 (MQ=255)
                                gTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAACCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGActt                                                                                                                       >  1:61016/1‑150 (MQ=255)
                                  cAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtct                                                                                                                     >  2:121617/1‑150 (MQ=255)
                                  cAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtct                                                                                                                     >  1:68391/1‑150 (MQ=255)
                                    aCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGCTGGTAAGACCCCCCGGGCTTTTCGCATTACAGCGActtctaa                                                                                                                   >  1:510994/1‑148 (MQ=255)
                                     cccGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACAACACGGTCTTTTCTCATTACAGCGActtatacc                                                                                                                  >  2:828974/1‑145 (MQ=255)
                                           cTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCAt                                                                                                            <  2:862653/150‑1 (MQ=255)
                                             gtgGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAc                                                                                                          <  2:300708/150‑1 (MQ=255)
                                             gtgGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATAc                                                                                                          >  2:566453/1‑150 (MQ=255)
                                              tgGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCAAACt                                                                                                         >  2:1340688/1‑150 (MQ=255)
                                               ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                        >  1:936263/1‑150 (MQ=255)
                                                gCCTGAACACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAACTCCAAACCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCACACCGTACTGc                                                                                                       >  1:98444/1‑150 (MQ=255)
                                                 ccTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACTGCGACTTCTTCCCACCATACTGCt                                                                                                      >  2:1550634/1‑150 (MQ=255)
                                                 ccTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                                                                      <  1:1081955/150‑1 (MQ=255)
                                                 ccTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                                                                      >  2:243551/1‑150 (MQ=255)
                                                 ccTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                                                                      <  2:261038/150‑1 (MQ=255)
                                                  cTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCtt                                                                                                     <  1:1202227/150‑1 (MQ=255)
                                                    gAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACGTCTTCCCACCATACTGCTTaa                                                                                                   >  1:671957/1‑150 (MQ=255)
                                                    gAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCCTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaa                                                                                                   >  1:1118111/1‑150 (MQ=255)
                                                     aaagCGCGCGGCTGTCCAGAAACCCGCGCGGATCTGGTGAACATGGTTCAGGCGGCGGAGAACTTGCGTCGGGGGTAAAATCCAATCCGGGTGGTAATACGTCCGGGTCTTTTCTGATTACAGCGACTTCTTCCCTCCATACTGCTTaaa                                                                                                  <  1:772474/150‑1 (MQ=255)
                                                       accgcgcgGCGGGCCATAAACCCGCGCGGAACTGGCTAACATGGTGCAGGCGGCGCATAAATTGCGTCGGGGTAAAAATCCAAACCGGGTGGAAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                                                                                <  1:1758702/148‑1 (MQ=255)
                                                        aCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCa                                                                                               >  2:292579/1‑150 (MQ=255)
                                                         cgcgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCAt                                                                                              >  1:1682165/1‑150 (MQ=255)
                                                                   gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGAAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACACGGTATTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATAc                                                                                    >  1:1588386/1‑150 (MQ=255)
                                                                    ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGCCTTCTTCCCACCATACTGCTTAACCCATTCCAGCATACg                                                                                   >  1:647853/1‑150 (MQ=255)
                                                                    ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACg                                                                                   >  2:1426008/1‑150 (MQ=255)
                                                                       gAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTg                                                                                >  2:207412/1‑150 (MQ=255)
                                                                         aaCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCCGTCTTTTCTCATAACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgcc                                                                              >  2:710121/1‑150 (MQ=255)
                                                                           cccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                            >  2:760366/1‑150 (MQ=255)
                                                                                   gCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttc                                                                    >  1:1615109/1‑150 (MQ=255)
                                                                                    cACTGGCGAACATGGGGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCAtcttct                                                                   <  2:406233/150‑1 (MQ=255)
                                                                                      cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCGCATTACAGCGACTTCTTCCCCCCATACTGCTGAACCCATTCCGGCATCCGCTGCCAGCCATCTTCTgg                                                                 >  1:863745/1‑149 (MQ=255)
                                                                                        ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGTCAGCCATCTTCTGCAg                                                               >  1:197694/1‑150 (MQ=255)
                                                                                        ggCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGCGGTAAAATCCAAACCGGGTGGTAATACAACCCGGTCTGTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAACCCATTCCAGCATAAGCTGACAGCCAGCGTCTGCAg                                                               >  1:1525107/1‑150 (MQ=255)
                                                                                          cGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGAt                                                             >  2:1407594/1‑150 (MQ=255)
                                                                                          cGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACAATAATGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGaa                                                             >  2:1414842/1‑149 (MQ=255)
                                                                                             aCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCCGCATACGCTGCCAGCCATCTTCTGCGGATGCg                                                          >  1:1066560/1‑150 (MQ=255)
                                                                                             aCATGGTCCAGGCGGCGCCGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg                                                          <  1:1390910/150‑1 (MQ=255)
                                                                                              cATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg                                                         >  1:1569834/1‑150 (MQ=255)
                                                                                               aTGGTGCAGGCGGCGCAGAACTTGCGTCGGGGTTAAAATCCAAACCGTGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGc                                                        <  1:56568/150‑1 (MQ=255)
                                                                                                tGGTGCAGGCGGCGCAGAACTTGCGGCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa                                                       <  2:1072427/150‑1 (MQ=255)
                                                                                                    gCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATATCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGAt                                                   <  2:952919/150‑1 (MQ=255)
                                                                                                    gCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGAt                                                   <  2:885675/150‑1 (MQ=255)
                                                                                                     cAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATa                                                  >  2:683446/1‑150 (MQ=255)
                                                                                                      aGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAg                                                 <  1:348584/150‑1 (MQ=255)
                                                                                                       ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGc                                                >  1:1267090/1‑150 (MQ=255)
                                                                                                           gcgcAGAACTTGCGTCGGGGGTAAAATACAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCgg                                            <  1:1278858/150‑1 (MQ=255)
                                                                                                            cgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCggg                                           <  1:397979/150‑1 (MQ=255)
                                                                                                            cgcAGAACTTCCGTCCGGCGTAAAATCCAAACCGCGTGGTAATACCACCCGGTCTTTTCTCATAACAGCGACTTCTTCCCACCATACTGCTTACACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCggg                                           >  2:847155/1‑150 (MQ=255)
                                                                                                             gcATAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCTGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGc                                          <  1:1326241/150‑1 (MQ=255)
                                                                                                             gcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGc                                          >  2:1289987/1‑150 (MQ=255)
                                                                                                                 aaCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATa                                      <  1:759697/150‑1 (MQ=255)
                                                                                                                         tCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGt                              <  2:1758367/150‑1 (MQ=255)
                                                                                                                         tCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGt                              <  2:398212/150‑1 (MQ=255)
                                                                                                                           gggggTAAAATCCAAACCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCTCCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTg                            <  2:1713689/150‑1 (MQ=255)
                                                                                                                            ggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATGCTGCTTACACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGAGTTGa                           >  1:75103/1‑150 (MQ=255)
                                                                                                                             gggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGaa                          <  1:1109127/150‑1 (MQ=255)
                                                                                                                              ggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAat                         <  2:171820/150‑1 (MQ=255)
                                                                                                                               ggaaaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatg                        <  2:196665/148‑1 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTGGAatg                        >  1:1522286/1‑150 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatg                        <  1:1566925/150‑1 (MQ=255)
                                                                                                                                   aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgcatg                    <  2:265104/150‑1 (MQ=255)
                                                                                                                                   aaTCCAAACCGGGTGGTAATACCACCAGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGACAGCTTCTGCAGATGCGGCATGATAGCGCGGGCGAGAATCAGCGTTGActgcatg                    >  1:463370/1‑150 (MQ=255)
                                                                                                                                      ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGccc                 <  1:134693/150‑1 (MQ=255)
                                                                                                                                      aaagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGccc                 >  1:60620/5‑150 (MQ=255)
                                                                                                                                       cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCg                >  1:1490687/1‑150 (MQ=255)
                                                                                                                                          aCCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCg             <  1:646251/150‑1 (MQ=255)
                                                                                                                                           ccGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGt            >  1:1076227/1‑150 (MQ=255)
                                                                                                                                           ccGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGt            >  1:782734/1‑150 (MQ=255)
                                                                                                                                           ccGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGt            <  2:1283518/150‑1 (MQ=255)
                                                                                                                                            cGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTc           <  2:1391342/150‑1 (MQ=255)
                                                                                                                                                tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTTCTTAAACCATTCCAGCATACGCTTCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCggg       <  1:1450094/150‑1 (MQ=255)
                                                                                                                                                tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCggg       <  2:1468675/150‑1 (MQ=255)
                                                                                                                                                gggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCggg       <  2:1610887/149‑1 (MQ=255)
                                                                                                                                                 ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGt      >  1:618757/1‑150 (MQ=255)
                                                                                                                                                   tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTac    <  1:1268963/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACAC  >  NC_000913/4015184‑4015482

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