Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,069,320 A→G S463G (AGT→GGT)  purH → bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/inosine monophosphate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,069,3200AG99.8% 2888.5 / NA 1094S463G (AGT→GGT) purHbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/inosine monophosphate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (1064/27);  minor base C (2/0);  total (1067/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1099 total aligned reads displayed.

CAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGAGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAATTTTAAA  >  NC_007793/1069223‑1069419
                                                                                                 |                                                                                                   
cAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGt                                                                                                   >  1:2677777/1‑100 (MQ=255)
       cATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGa                                                                                            >  1:6517136/1‑100 (MQ=255)
       cATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTAGATCgg                                                                                            >  1:1538540/1‑99 (MQ=255)
       cATGTAGCGTTAGTATCTGATGGATTTCTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGa                                                                                            >  1:5069853/1‑100 (MQ=255)
        aTGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGAt                                                                                           >  1:4541356/1‑100 (MQ=255)
          gTAGCGTTAGTATCTGATGGGTTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTa                                                                                         >  1:6287054/1‑100 (MQ=255)
          atAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTa                                                                                         >  1:3796807/2‑100 (MQ=255)
           gaGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTaa                                                                                        <  1:3975526/99‑1 (MQ=255)
                         gATGGGTTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:835486/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATcggaagagca                                                                          >  1:2861422/1‑91 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:1035149/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:3496589/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:5265093/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:5682620/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:5736820/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:617428/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:6467832/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:1357433/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:1288529/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:1212286/1‑100 (MQ=255)
                                  ttCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATgg                                                                 >  1:1965834/1‑100 (MQ=255)
                                  ttCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATgg                                                                 >  1:427977/1‑100 (MQ=255)
                                    ctcTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCt                                                               >  1:2367635/3‑100 (MQ=255)
                                     aagaTGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCta                                                              >  1:2996775/4‑100 (MQ=255)
                                       tctGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAGTTGATATGGCtaat                                                            >  1:4805699/3‑100 (MQ=255)
                                       tATGGGAGATACAGTTGGACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaat                                                            >  1:1855700/1‑100 (MQ=255)
                                        aTGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaata                                                           >  1:4373685/1‑100 (MQ=255)
                                        aTGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaata                                                           >  1:4033679/1‑100 (MQ=255)
                                        aTGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaata                                                           >  1:1154929/1‑100 (MQ=255)
                                         tGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaataa                                                          >  1:5314566/1‑100 (MQ=255)
                                             agaTACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACAt                                                      >  1:5140410/1‑100 (MQ=255)
                                                           ttGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAGTGGt                                        <  1:4615903/100‑1 (MQ=255)
                                                           atGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGt                                        >  1:5918807/2‑100 (MQ=255)
                                                           aggcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATgat                                        >  1:4095353/3‑98 (MQ=255)
                                                           aggcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGt                                        >  1:518294/3‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATggga                                       >  1:1680053/1‑98 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:2433563/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:2764871/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:5838176/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:4128776/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:596341/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3664988/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3605527/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:2815490/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:146793/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:1078545/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3048265/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3088574/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGCATAAAGGCAATTATCCAACCGGGTGGTCCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:655918/1‑100 (MQ=255)
                                                            cccagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3903555/3‑100 (MQ=255)
                                                             gcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAg                                      >  1:5642866/1‑100 (MQ=255)
                                                             ccagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAg                                      >  1:73470/2‑100 (MQ=255)
                                                              cagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCATTGGTAga                                     >  1:1426251/1‑99 (MQ=255)
                                                              cagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGt                                     >  1:3845329/1‑100 (MQ=255)
                                                              cagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATACAATTGATATGGCTAATAAACATGGTATTGCAATGGCAGt                                     >  1:781063/1‑100 (MQ=255)
                                                               agtgcaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCGATGGTAgat                                    >  1:6071978/5‑98 (MQ=255)
                                                                gcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTca                                   >  1:3556625/1‑100 (MQ=255)
                                                                 cacaACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcac                                  >  1:4222300/1‑100 (MQ=255)
                                                                 cacaACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcac                                  >  1:1004986/1‑100 (MQ=255)
                                                                  gcaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcaca                                 >  1:4476205/2‑100 (MQ=255)
                                                                  acaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcaca                                 >  1:6410014/1‑100 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAg                                >  1:6616079/1‑100 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAg                                >  1:4885572/1‑100 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAg                                >  1:497211/1‑100 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAg                                >  1:4738430/1‑100 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAg                                >  1:4416402/1‑100 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAACGGTAgatcgga                                >  1:2941140/1‑94 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGAATTGCAATGGTAGTCACAg                                >  1:3710812/1‑100 (MQ=255)
                                                                     acaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAgatcggaag                              >  1:3288707/1‑92 (MQ=255)
                                                                     acaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGc                              >  1:2502103/1‑100 (MQ=255)
                                                                      actGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCa                             >  1:6340229/3‑100 (MQ=255)
                                                                       ctGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCAc                            >  1:33392/2‑100 (MQ=255)
                                                                       aTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCAc                            >  1:5491220/1‑100 (MQ=255)
                                                                         ggTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGGATTGCAATGGTAGTCACAGGCACTc                          >  1:3418510/1‑100 (MQ=255)
                                                                           acactgggCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGa                        >  1:5771642/7‑100 (MQ=255)
                                                                            ataAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGAc                       >  1:2166167/1‑100 (MQ=255)
                                                                              caaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAt                     >  1:4172711/2‑100 (MQ=255)
                                                                               aaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACGtt                    >  1:6213569/1‑100 (MQ=255)
                                                                               aaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtt                    >  1:2288358/1‑100 (MQ=255)
                                                                                aGGCAATTATCCAATCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAttt                   <  1:204513/100‑1 (MQ=255)
                                                                                aGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAttt                   >  1:5564740/1‑100 (MQ=255)
                                                                                 ggCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtttt                  >  1:1600284/1‑100 (MQ=255)
                                                                                 ggCAATTATCCAACCGGGCGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtttt                  >  1:3136889/1‑100 (MQ=255)
                                                                                  ccAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTa                 >  1:2107781/2‑100 (MQ=255)
                                                                                        acccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACt           >  1:1925043/3‑100 (MQ=255)
                                                                                        aTCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACt           >  1:6128570/1‑100 (MQ=255)
                                                                                        aTCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACt           >  1:3370929/1‑100 (MQ=255)
                                                                                         cccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTa          >  1:5424874/2‑100 (MQ=255)
                                                                                          ccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTaa         >  1:2213522/1‑100 (MQ=255)
                                                                                          ccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTaa         >  1:2548491/1‑100 (MQ=255)
                                                                                           cAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAt        >  1:677488/1‑100 (MQ=255)
                                                                                            gaCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:6021277/2‑100 (MQ=255)
                                                                                            aaCCGGGTGTCTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGGATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:4930432/1‑100 (MQ=255)
                                                                                             gccGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:5011597/2‑100 (MQ=255)
                                                                                             aCCGGGTGGTTCGTTTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:910110/1‑100 (MQ=255)
                                                                                             aCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:1723534/1‑100 (MQ=255)
                                                                                             aCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:4666869/1‑100 (MQ=255)
                                                                                             aCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:2599359/1‑100 (MQ=255)
                                                                                             aCCGGGTGGCTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:3228915/1‑100 (MQ=255)
                                                                                              ccGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtttt     >  1:1814336/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGAGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAATTTTAAA  >  NC_007793/1069223‑1069419

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