Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,860,748 C→T 7.0% M77I (ATG→ATA) ‡ purN ← phosphoribosylglycinamide formyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,860,7480CT 7.0% 474.7 / 16.4 156M77I (ATG→ATA) ‡purNphosphoribosylglycinamide formyltransferase
Reads supporting (aligned to +/- strand):  ref base C (79/66);  new base T (7/4);  total (86/70)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

GTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCC  >  NC_002947/1860605‑1860896
                                                                                                                                                |                                                                                                                                                    
gTATGCAGGCCCTTGTACTTGGGCTGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGTATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAgc                                                                                                                                                <  2:2692938/150‑1 (MQ=255)
    gCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCgg                                                                                                                                            >  1:351112/1‑150 (MQ=255)
      aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCCGCGCTAATGGcgcgca                                                                                                                                           <  2:4620359/150‑3 (MQ=255)
      aGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCa                                                                                                                                          >  1:2172735/1‑150 (MQ=255)
       ggCCCTTGGACTGGGGCAGCAGCGCCGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCTAATGGcgcgcatc                                                                                                                                         <  2:3459258/150‑5 (MQ=255)
        gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                        >  2:2124156/1‑150 (MQ=255)
        gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                        <  1:3114975/150‑1 (MQ=255)
        gCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATc                                                                                                                                        <  1:2108587/150‑1 (MQ=255)
         cccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCTAATGGCGCGCATCGa                                                                                                                                      <  1:2979634/150‑1 (MQ=255)
          ccTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGa                                                                                                                                      <  1:1483231/150‑1 (MQ=255)
             tGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGTCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAgg                                                                                                                                   <  2:3066039/150‑1 (MQ=255)
             tGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAgg                                                                                                                                   <  2:1121533/150‑1 (MQ=255)
                aCTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCtt                                                                                                                                >  1:432033/1‑150 (MQ=255)
                 cTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTc                                                                                                                               >  1:1875362/1‑150 (MQ=255)
                 cTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTc                                                                                                                               >  1:313858/1‑150 (MQ=255)
                  ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAATCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                              >  2:212113/1‑150 (MQ=255)
                  ttGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCa                                                                                                                              >  1:1892078/1‑150 (MQ=255)
                     ggCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACgg                                                                                                                           >  1:4638864/1‑150 (MQ=255)
                       cagcagCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGcc                                                                                                                         <  1:11655/150‑1 (MQ=255)
                             cGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGaa                                                                                                                   <  2:3242001/150‑1 (MQ=255)
                              gACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAAc                                                                                                                  >  1:3689316/1‑150 (MQ=255)
                              gACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAAc                                                                                                                  >  1:4276141/1‑150 (MQ=255)
                                cGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTg                                                                                                                <  2:4225157/150‑1 (MQ=255)
                                 gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                               >  1:4491600/1‑150 (MQ=255)
                                 gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                               >  1:828862/1‑150 (MQ=255)
                                 gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                               <  2:3861481/150‑1 (MQ=255)
                                 gggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTgg                                                                                                               <  2:1354868/150‑1 (MQ=255)
                                  ggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTggg                                                                                                              >  2:3672454/1‑150 (MQ=255)
                                  ggTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTggg                                                                                                              >  2:1091861/1‑150 (MQ=255)
                                      gATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                          <  2:3524538/150‑1 (MQ=255)
                                      gATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGtgtg                                                                                                          <  1:389120/150‑1 (MQ=255)
                                       aTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTgg                                                                                                         <  2:3373887/150‑1 (MQ=255)
                                       aTGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTgg                                                                                                         <  2:4708018/150‑1 (MQ=255)
                                         gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                       <  2:2507940/150‑1 (MQ=255)
                                         gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                       <  1:1237328/150‑1 (MQ=255)
                                         gTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTc                                                                                                       <  1:2758684/150‑1 (MQ=255)
                                                 ggcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacga                                                                                               <  2:1133871/150‑1 (MQ=255)
                                                 ggcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacga                                                                                               >  1:3808721/1‑150 (MQ=255)
                                                  gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                              >  1:2269583/1‑150 (MQ=255)
                                                  gcggcCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCacgac                                                                                              >  1:4054962/1‑150 (MQ=255)
                                                      cccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTg                                                                                          <  2:1529158/150‑1 (MQ=255)
                                                       ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                         >  2:3951341/1‑150 (MQ=255)
                                                       ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                         <  2:4694526/150‑1 (MQ=255)
                                                       ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                         >  2:60024/1‑150 (MQ=255)
                                                       ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                         >  1:3723844/1‑150 (MQ=255)
                                                       ccTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGt                                                                                         <  2:4826386/150‑1 (MQ=255)
                                                        cTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTc                                                                                        >  1:1415336/1‑150 (MQ=255)
                                                               tGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccg                                                                                 >  1:63242/1‑150 (MQ=255)
                                                                   gcACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATgccggccg                                                                             >  1:3279343/1‑150 (MQ=255)
                                                                       gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCCGCACCACCAGGTCCGGGGCCAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGCAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                         >  2:4103480/1‑150 (MQ=255)
                                                                       gAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgc                                                                         >  1:4289677/1‑150 (MQ=255)
                                                                          gccgccACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTgcggc                                                                      <  2:2707179/150‑1 (MQ=255)
                                                                               cACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGtt                                                                 >  1:719873/1‑150 (MQ=255)
                                                                                  tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                              <  2:2787580/150‑1 (MQ=255)
                                                                                  tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                              >  2:2383802/1‑150 (MQ=255)
                                                                                  tGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCa                                                              <  2:3274033/150‑1 (MQ=255)
                                                                                    aGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAgg                                                            <  2:4215246/150‑1 (MQ=255)
                                                                                         gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                       <  2:1278031/150‑1 (MQ=255)
                                                                                         gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                       <  2:4904340/150‑1 (MQ=255)
                                                                                         gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                       <  1:1999577/150‑1 (MQ=255)
                                                                                         gcgcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTa                                                       >  1:4201806/1‑150 (MQ=255)
                                                                                           gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                     >  2:4283546/1‑150 (MQ=255)
                                                                                           gcATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAgg                                                     >  1:655532/1‑150 (MQ=255)
                                                                                              tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                  >  1:1109252/1‑150 (MQ=255)
                                                                                              tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                  >  1:2450594/1‑150 (MQ=255)
                                                                                              tGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCAt                                                  >  1:4953372/1‑150 (MQ=255)
                                                                                               gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                 <  1:3532936/150‑1 (MQ=255)
                                                                                               gAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATc                                                 <  1:3532880/150‑1 (MQ=255)
                                                                                                   ccGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgcg                                             >  1:2516350/1‑150 (MQ=255)
                                                                                                    cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgcgc                                            >  2:531816/1‑150 (MQ=255)
                                                                                                    cGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGgcgc                                            >  1:919698/1‑150 (MQ=255)
                                                                                                       ccAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                         >  1:904358/1‑150 (MQ=255)
                                                                                                       ccAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGt                                         >  1:1810691/1‑150 (MQ=255)
                                                                                                        cAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGtt                                        <  2:2003793/150‑1 (MQ=255)
                                                                                                         aGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTg                                       >  1:1033044/1‑150 (MQ=255)
                                                                                                            accaccAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaa                                    <  1:5049851/150‑1 (MQ=255)
                                                                                                             ccaccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGaaa                                   <  2:749669/150‑1 (MQ=255)
                                                                                                              caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                  >  2:3178691/1‑150 (MQ=255)
                                                                                                              caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAg                                  >  2:3178233/1‑150 (MQ=255)
                                                                                                              caccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAAGCATCGGCGCGGTTGGAAAg                                  >  1:2733374/1‑150 (MQ=255)
                                                                                                                ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                >  2:3782237/1‑150 (MQ=255)
                                                                                                                ccaGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCa                                <  2:2633359/150‑1 (MQ=255)
                                                                                                                 caGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAc                               >  2:4224874/1‑150 (MQ=255)
                                                                                                                 caGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAc                               >  1:2224326/1‑150 (MQ=255)
                                                                                                                  aGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCAcc                              <  2:2865399/150‑1 (MQ=255)
                                                                                                                    gTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgc                            >  1:266602/1‑150 (MQ=255)
                                                                                                                       cGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACcgcgcg                         <  1:266029/150‑1 (MQ=255)
                                                                                                                         gggCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGa                       <  2:5112432/150‑1 (MQ=255)
                                                                                                                         gggCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGa                       <  2:943967/150‑1 (MQ=255)
                                                                                                                             aaaGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACg                   <  2:9179/150‑1 (MQ=255)
                                                                                                                                gCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc                >  2:1672244/1‑150 (MQ=38)
                                                                                                                                gCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc                >  1:943088/1‑150 (MQ=38)
                                                                                                                                gCCTTCGATGCCGCGCTAATGGCGCGCATCGACGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCAc                >  1:617104/1‑150 (MQ=38)
                                                                                                                                  cGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCg              >  1:3254959/1‑150 (MQ=255)
                                                                                                                                   ttCGATGCCGCGCTAATGGCGCGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg             >  2:1036536/2‑150 (MQ=255)
                                                                                                                                   gTCGATGCGCGCCCTTCGCGCGGCAGCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGCTGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg             <  1:3848037/150‑1 (MQ=255)
                                                                                                                                   gTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg             >  2:4041511/1‑150 (MQ=255)
                                                                                                                                   gTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCgg             >  2:2438363/1‑150 (MQ=255)
                                                                                                                                      gatgccgcgcTAATGGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGc           >  2:2981580/6‑150 (MQ=255)
                                                                                                                                               ccATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGcc  <  2:69581/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GTATGCAGGCCCTTGTACTTGGGCAGCAGCGACGGGTGGATGTTGAGCAGGCGGCCCTGATAGTGGCGCACGAAGCCGCCACTGAGGATGCGCATGAAACCGGCCAGCACCACCAGGTCCGGGGCAAAGCCGTCGATGCGCGCCATTAGCGCGGCATCGAAGGCTTCACGGCCGTCGAACTGGGTGTGGTCGAGCACGACGCTGTCGATGCCGGCCGCTGCGGCGCGTTGCAGGCCGTAGGCATCGGCGCGGTTGGAAAGCACCGCGCGGATACGCACCGGGCTGTCCTGCC  >  NC_002947/1860605‑1860896

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