Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,226,762 G→A 6.2% S228F (TCC→TTC)  PP_1964 ← dGMP/dTMP/dhydroxyCMP kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,226,7620GA 6.2% 183.2 / 2.6 64S228F (TCC→TTC) PP_1964dGMP/dTMP/dhydroxyCMP kinase
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TAGTACTTCTTTTGTGAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGAGTAGGGGTTGA  >  NC_002947/2226613‑2226910
                                                                                                                                                     |                                                                                                                                                    
tAGTACTTCTTTTGTGAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATggg                                                                                                                                                      <  1:944617/150‑1 (MQ=255)
   tACTTCTTTTGTGAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGAtt                                                                                                                                                   >  1:364460/1‑150 (MQ=255)
         ttttGTGAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGTGAGTGTTCtt                                                                                                                                             >  2:48729/1‑150 (MQ=255)
         ttttGTGAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACtt                                                                                                                                             >  1:229810/1‑150 (MQ=255)
             gtgAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTAt                                                                                                                                         <  1:1174757/150‑1 (MQ=255)
                 aaTTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGACATATCGCCGTTTGGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCa                                                                                                                                     >  1:415860/1‑150 (MQ=255)
                 aaTTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCa                                                                                                                                     <  1:1266983/150‑1 (MQ=255)
                 aaTTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCa                                                                                                                                     <  2:839970/150‑1 (MQ=255)
                 aaTTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCa                                                                                                                                     <  1:1801735/150‑1 (MQ=255)
                 aaTTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCa                                                                                                                                     <  2:287479/150‑1 (MQ=255)
                  aTTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCAt                                                                                                                                    <  2:837419/150‑1 (MQ=255)
                                tCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGata                                                                                                                      >  2:1916772/1‑150 (MQ=255)
                                        aaCGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAAt                                                                                                              >  1:691328/1‑150 (MQ=255)
                                         aCGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGa                                                                                                                                                     >  1:128432/1‑110 (MQ=255)
                                         aCGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGa                                                                                                                                                     <  2:128432/110‑1 (MQ=255)
                                           gCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTAc                                                                                                           >  1:1241030/1‑150 (MQ=255)
                                                tCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCAt                                                                                                      >  1:518749/1‑150 (MQ=255)
                                                  cAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCTTTGCGGAATTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGTTATGCCAAATTACTCCATTa                                                                                                    >  2:662507/1‑150 (MQ=255)
                                                   aaGCGGATTCAACCTTTTTTCTTAGTTCTTCGATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTAt                                                                                                   >  2:1233773/1‑150 (MQ=255)
                                                            cAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATggagatcgg                                                                                          >  1:1851652/1‑143 (MQ=255)
                                                             aaCCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCa                                                                                         <  2:229810/150‑1 (MQ=255)
                                                             aaCCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTACTGAACTCAATTCGATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCa                                                                                         >  1:707550/1‑150 (MQ=255)
                                                              aCCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGa                                                                                                                                                     <  1:1756002/89‑1 (MQ=255)
                                                              aCCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGa                                                                                                                                                     >  2:1756002/1‑89 (MQ=255)
                                                                cTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGaa                                                                                      <  2:364460/150‑1 (MQ=255)
                                                                       cTTAGTTCTTCAATTTTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCt                                                                               <  1:509956/150‑1 (MQ=255)
                                                                        ttAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCCTTTGTTCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCtt                                                                              <  1:48729/150‑1 (MQ=255)
                                                                         tAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGa                                                                                                                            <  2:1473417/103‑1 (MQ=255)
                                                                         tAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGa                                                                                                                            >  1:1473417/1‑103 (MQ=255)
                                                                               ttcAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGttt                                                                         >  1:480211/1‑148 (MQ=255)
                                                                               ttcAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGttt                                                                         <  2:480211/148‑1 (MQ=255)
                                                                                tcAATTGTATCATTGTTTTCAATTGAAATATCGCCCTTGCGGAACTCAATTCCATCTTCTTAAGAATGGGATTGCACTTTTATTGCATTGTTGCCAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCTCTTCGTTTTCa                                                                      <  1:426607/150‑1 (MQ=255)
                                                                                  aaTTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGtttt                                                                        <  1:244035/146‑1 (MQ=255)
                                                                                  aaTTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGtttt                                                                        >  2:244035/1‑146 (MQ=255)
                                                                                  aaTTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGATGCGACGGAGTGCCGCTTTATTTTcttc                                                                    >  2:1356867/1‑147 (MQ=255)
                                                                                    ttGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATc                                                                  <  1:1356867/150‑1 (MQ=255)
                                                                                       tATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGa                                                               >  2:802732/1‑150 (MQ=255)
                                                                                        aTCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGat                                                              <  2:691328/150‑1 (MQ=255)
                                                                                         tCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGata                                                             >  1:659512/1‑150 (MQ=255)
                                                                                           aTTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATc                                                           >  1:637650/1‑150 (MQ=255)
                                                                                               ttttCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTta                                                       >  2:1217580/1‑150 (MQ=255)
                                                                                                tttCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTtaa                                                      <  1:1422907/150‑1 (MQ=255)
                                                                                                     aTTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATg                                                 >  2:1829396/1‑150 (MQ=255)
                                                                                                       tGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCa                                                                                         <  2:862966/108‑1 (MQ=255)
                                                                                                       tGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCa                                                                                         >  1:862966/1‑108 (MQ=255)
                                                                                                       tGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGaa                                               <  1:112153/150‑1 (MQ=255)
                                                                                                                gCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAg                                                                                                                           <  1:1300913/65‑1 (MQ=255)
                                                                                                                gCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAg                                                                                                                           >  2:1300913/1‑65 (MQ=255)
                                                                                                                  cGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTAc                                    <  1:802732/150‑1 (MQ=255)
                                                                                                                    ttGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCCGTGGTAGCAGTAACATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCtagatcgg                                  >  2:1418441/1‑143 (MQ=255)
                                                                                                                      gCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGc                                <  1:1201969/150‑1 (MQ=255)
                                                                                                                            cTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCAtt                          <  2:121350/150‑1 (MQ=255)
                                                                                                                               aaTTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGt                       <  2:1440334/150‑1 (MQ=255)
                                                                                                                               aaTTCCATCTTCTGAAGAATGGAATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGaa                                               >  2:1844857/1‑126 (MQ=255)
                                                                                                                               aaTTCCATCTTCTGAAGAATGGAATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGaa                                               >  2:1844076/1‑126 (MQ=255)
                                                                                                                               aaTTCCATCTTCTGAAGAATGGAATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGaa                                               <  1:1844857/126‑1 (MQ=255)
                                                                                                                               aaTTCCATCTTCTGAAGAATGGAATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGaa                                               <  1:1844076/126‑1 (MQ=255)
                                                                                                                                   ccATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGt                   <  1:1233773/150‑1 (MQ=255)
                                                                                                                                    cATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTg                  >  1:1855949/1‑150 (MQ=255)
                                                                                                                                     aTCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTgg                 >  2:1257112/1‑150 (MQ=255)
                                                                                                                                      tcttctGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCAGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTggg                <  2:850104/150‑1 (MQ=255)
                                                                                                                                         tctGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGa             <  2:637650/150‑1 (MQ=255)
                                                                                                                                           tGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGAga           >  2:764536/1‑149 (MQ=255)
                                                                                                                                                aaTGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGAGTAgggg      <  1:1257112/150‑1 (MQ=255)
                                                                                                                                                  tGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGAGTAGGGGtt    >  2:1061843/1‑150 (MQ=255)
                                                                                                                                                    ggATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGAGTAGGGGTTGa  <  2:1241030/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGTACTTCTTTTGTGAAATTCATTGCATTAATCGACGCGAACGCAGCTCCAAGCGGATTCAACCTTTTTTCTTAGTTCTTCAATTGTATCATTGTTTTCAATTGAAATATCGCCGTTGCGGAACTCAATTCCATCTTCTGAAGAATGGGATTGTACTTTTATTGCATTGTTGCGAGTGATATGCCAAATTACTCCATTATGGCTGCGCCAGAATGCCGCTTCGTTTTCATATCTGATATCCTTAATAATGAAGATGCTTTTACCAGCGGCATTTGTTAGTGGGCGAGTAGGGGTTGA  >  NC_002947/2226613‑2226910

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