Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,373 C→G 7.6% intergenic (+30/+13) crr → / ← pdxK glucose‑specific enzyme IIA component of PTS/pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,536,3730CG 7.6% 366.9 / 17.9 118intergenic (+30/+13)crr/pdxKglucose‑specific enzyme IIA component of PTS/pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (55/53);  new base G (3/6);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 118 total aligned reads displayed.

AGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACTGA  >  NC_000913/2536226‑2536519
                                                                                                                                                   |                                                                                                                                                  
aGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCgg                                                                                                                                                  >  1:1386610/1‑150 (MQ=255)
 gCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc                                                                                                                                                 <  1:1207532/150‑1 (MQ=255)
 gCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc                                                                                                                                                 <  1:294015/150‑1 (MQ=255)
 gCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc                                                                                                                                                 >  2:1552131/1‑150 (MQ=255)
 gCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc                                                                                                                                                 >  2:491703/1‑150 (MQ=255)
    aaGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGcc                                                                                                                                              >  2:1626135/1‑150 (MQ=255)
     aGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATTGCGCCCATCGGCGCCa                                                                                                                                             >  1:593587/1‑150 (MQ=255)
     aGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGGGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCa                                                                                                                                             <  1:1458466/150‑1 (MQ=255)
       tCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtt                                                                                                                                           <  1:974662/150‑1 (MQ=255)
        cTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttt                                                                                                                                          <  2:175095/150‑1 (MQ=255)
         tACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtttt                                                                                                                                         <  1:650153/150‑1 (MQ=255)
          aCCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTACCCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttttt                                                                                                                                        <  1:192398/150‑1 (MQ=255)
          aCCCTGACTCCAGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAttttt                                                                                                                                        <  2:1002018/150‑1 (MQ=255)
           cccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt                                                                                                                                       >  1:862404/1‑150 (MQ=255)
            ccTGACTCCGCTTGTTATCTCCAACATGGACGAAATCAAAGAACAGATCAAACAGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAttttttt                                                                                                                                      <  1:279787/150‑1 (MQ=255)
             cTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtttttcac                                                                                                                                     <  1:749396/150‑4 (MQ=255)
             cTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtttttcac                                                                                                                                     <  1:572576/150‑4 (MQ=255)
               gACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATg                                                                                                                                   <  1:407742/150‑1 (MQ=255)
                aCTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATgt                                                                                                                                  <  1:136561/150‑2 (MQ=255)
                   ccGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCTCCCATCGGCGCCATTTTTTTATGCTTc                                                                                                                               >  1:874458/1‑150 (MQ=255)
                     ggTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCg                                                                                                                             <  1:563194/150‑1 (MQ=255)
                      gttgttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGc                                                                                                                            >  1:1239828/1‑150 (MQ=255)
                            aTCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcgg                                                                                                                      >  1:204344/1‑150 (MQ=255)
                              ctcCAACATGGACGAAATCAAAGAACTGATCAAACTGTCAGTTAGCGTAAACGTGGGTGAAACCCAGGTTATCCGCATCAAGAAGTAATTATTGCCGCAGTGAAAAATGGCGCCCATAGGCGCCATTTTTTTATGCTTCCGCCAgcggcg                                                                                                                    >  2:761538/1‑150 (MQ=255)
                               tcaaaCATGGACGAAATCAAAGAACTGATCAAACTGTCAGGTAGCGTAACCGTGGGTGAAACCCCGGTTATACACATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTATATGCATCCTACAgcgtcga                                                                                                                   >  1:869272/4‑149 (MQ=255)
                                ccAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCggcggc                                                                                                                  >  2:1728538/1‑150 (MQ=255)
                                ccAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCggcggc                                                                                                                  >  2:101094/1‑150 (MQ=255)
                                  aaCATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaa                                                                                                                >  2:77299/1‑150 (MQ=255)
                                   aCATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCAGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaaa                                                                                                               >  1:265192/1‑150 (MQ=255)
                                      tGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATc                                                                                                            <  1:122044/150‑1 (MQ=255)
                                       ggACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATACTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCa                                                                                                           <  2:809373/150‑1 (MQ=255)
                                         aCGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAACACCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCCGCGCCAAATTTTTATGCTTCCGCCAGCGGCGGCAAAATCAAt                                                                                                         >  1:1277161/1‑150 (MQ=255)
                                          cGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAAtt                                                                                                        <  1:126167/150‑1 (MQ=255)
                                            aaaTCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCa                                                                                                      >  2:923919/1‑150 (MQ=255)
                                             aaTCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCAt                                                                                                     >  2:1480367/1‑150 (MQ=255)
                                              aTCAAAGAACTGATCAAACTGGCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGACGTAAATCATGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATc                                                                                                    <  2:321144/150‑1 (MQ=255)
                                                cAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGc                                                                                                  <  2:1016095/150‑1 (MQ=255)
                                                cAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGc                                                                                                  >  2:1074243/1‑150 (MQ=255)
                                                cAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGc                                                                                                  >  2:911230/1‑150 (MQ=255)
                                                cAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGc                                                                                                  <  1:1543277/150‑1 (MQ=255)
                                                 aaaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGct                                                                                                 <  1:1062708/150‑1 (MQ=255)
                                                 aaaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAATTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGct                                                                                                 >  1:1009974/1‑150 (MQ=255)
                                                      aCTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAt                                                                                            >  1:1630049/1‑150 (MQ=255)
                                                       cTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtg                                                                                           >  2:555415/1‑150 (MQ=255)
                                                        tGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATTGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgc                                                                                          <  1:826657/150‑1 (MQ=255)
                                                        tGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgc                                                                                          >  1:841120/1‑150 (MQ=255)
                                                         gATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTAATCGCTAACAtgct                                                                                         >  1:854275/1‑150 (MQ=255)
                                                             aaaCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGTCGCCCATCGGCGCCATTATTTTTTTCTTCCGCCAGCGGCGGCAAAATCAATTGATCGCTCTCATGctgctg                                                                                     >  1:182889/1‑150 (MQ=255)
                                                             aaaCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGctgctg                                                                                     >  2:1360093/1‑150 (MQ=255)
                                                               aCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTggg                                                                                   <  2:815791/150‑1 (MQ=255)
                                                                cTGTCCGGTAGCGTAACCGTGGGTGACACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt                                                                                  >  1:258057/1‑150 (MQ=255)
                                                                     cGGTAGCGTGACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTAAAAGCTCTCATGCTGCTGGGTGTAgc                                                                             >  2:1674876/1‑150 (MQ=255)
                                                                     cGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCTGATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgc                                                                             <  1:772337/150‑1 (MQ=255)
                                                                       gTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgcgc                                                                           <  1:491703/150‑1 (MQ=255)
                                                                        tAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCa                                                                          <  2:649859/150‑1 (MQ=255)
                                                                           cGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCa                                                                       <  1:524636/150‑1 (MQ=255)
                                                                              aaCCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACtt                                                                    >  1:545144/1‑150 (MQ=255)
                                                                              aaCCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACtt                                                                    >  1:176915/1‑150 (MQ=255)
                                                                                     ggTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTAc                                                             >  2:489804/1‑150 (MQ=255)
                                                                                             cccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAAcc                                                     <  1:1085342/150‑1 (MQ=255)
                                                                                                   tATCCGCATCAAGAAGTAATGCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGCGTGTAGCGCATCACTTCCAGGACGCGCAACCCCGCTc                                               <  1:1430808/150‑1 (MQ=255)
                                                                                                     tCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCATAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCgg                                             >  1:301731/1‑150 (MQ=255)
                                                                                                       ctcATAAAGAAGTAATTCTTGACGCAGTGTAAAATGGTGACCATCGGAGACATTTTTTTATGCTTCCGCCAGCTGAGGAAAAATCAATTCATCGCTCTCATGCTGCTGGTTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTg                                           <  2:102499/148‑1 (MQ=255)
                                                                                                       cGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTg                                           >  2:558063/1‑150 (MQ=255)
                                                                                                        gCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGc                                          >  1:591094/1‑150 (MQ=255)
                                                                                                          aTCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATAATTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATAACTTCCAGTACGCGCAACCCCGCTCGGTGCAc                                        >  2:614603/1‑150 (MQ=255)
                                                                                                           tCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACt                                       <  2:800263/150‑1 (MQ=255)
                                                                                                           tCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACt                                       <  1:642657/150‑1 (MQ=255)
                                                                                                               gaagTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCAt                                   >  1:719795/1‑150 (MQ=255)
                                                                                                                 agTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCg                                 >  1:1628008/1‑150 (MQ=255)
                                                                                                                 agTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCg                                 >  1:1618172/1‑150 (MQ=255)
                                                                                                                  gTACTTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTGCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGGAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                <  2:841153/150‑1 (MQ=255)
                                                                                                                  gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                >  2:1257852/1‑150 (MQ=255)
                                                                                                                  gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                <  1:1480367/150‑1 (MQ=255)
                                                                                                                  gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                <  2:353245/150‑1 (MQ=255)
                                                                                                                  gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTACAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                >  2:1542333/1‑150 (MQ=255)
                                                                                                                  gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTCGGTGTAGCGCATCACTTCCAGTACGCGCACCCCCGCTCGGTGCACTGCATAgg                                >  2:906479/1‑150 (MQ=255)
                                                                                                                   tAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGt                               >  1:1603225/1‑150 (MQ=255)
                                                                                                                        cTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACg                          >  2:970107/1‑150 (MQ=255)
                                                                                                                               cataaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTccc                   >  2:118679/6‑150 (MQ=255)
                                                                                                                               cAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTccc                   <  1:21514/150‑1 (MQ=255)
                                                                                                                                ataaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTATCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCt                  >  2:21054/5‑150 (MQ=255)
                                                                                                                                 gTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGTACTGCATCGGTTAACGCCTTCCCtt                 <  1:886037/150‑1 (MQ=255)
                                                                                                                                   taaaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               <  2:1230803/149‑1 (MQ=255)
                                                                                                                                   gAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               <  1:1284823/150‑1 (MQ=255)
                                                                                                                                   gAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               >  1:461581/1‑150 (MQ=255)
                                                                                                                                   gAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               <  2:545144/150‑1 (MQ=255)
                                                                                                                                   gAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               <  2:1691176/150‑1 (MQ=255)
                                                                                                                                     aaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAg             <  1:489804/150‑1 (MQ=255)
                                                                                                                                      aaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGc            <  1:1626135/150‑1 (MQ=255)
                                                                                                                                       aaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCa           >  1:915607/1‑150 (MQ=255)
                                                                                                                                          ggCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCAGTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGGTAACGCCTTCCCTTTCAGCAAGc        >  2:1206901/1‑150 (MQ=255)
                                                                                                                                           gcgcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGGTGGGTGTAGGGGATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGcc       <  2:258057/150‑1 (MQ=255)
                                                                                                                                           gcgcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCGTTCTGCAAGcc       >  2:1685846/1‑150 (MQ=255)
                                                                                                                                            cgcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCa      >  1:1375066/1‑150 (MQ=255)
                                                                                                                                             gcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCAc     >  2:1441008/1‑150 (MQ=255)
                                                                                                                                              cccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACt    <  1:1659871/150‑1 (MQ=255)
                                                                                                                                              cccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCGGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAACAAGCCACt    <  1:288640/150‑1 (MQ=255)
                                                                                                                                               ccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACtg   >  2:952713/1‑150 (MQ=255)
                                                                                                                                                cATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACtga  <  1:1441008/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACTGA  >  NC_000913/2536226‑2536519

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: