Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,752 G→T 8.1% intergenic (+57/+37) guaA → / ← PP_1033 glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,7520GT 8.1% 226.3 / 8.9 75intergenic (+57/+37)guaA/PP_1033glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base T (2/4);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

CGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTC  >  NC_002947/1178608‑1178899
                                                                                                                                                  |                                                                                                                                                    
cGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTg                                                                                                                                                  <  7:26384/149‑1 (MQ=255)
 ggTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGc                                                                                                                                                 <  8:295519/149‑1 (MQ=255)
    cGTATCATAAACTAAATCGACTGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGTTCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGAGGGATACACCAGAAAGCCCGCATGTTTGCgg                                                                                                                                               <  2:369129/148‑1 (MQ=255)
      tATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCt                                                                                                                                            <  6:353974/149‑1 (MQ=255)
          atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                        >  5:209754/1‑149 (MQ=255)
          atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTc                                                                                                                                         <  6:224439/148‑1 (MQ=255)
           tcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                       >  7:215220/1‑149 (MQ=255)
           tcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                       <  6:386822/149‑1 (MQ=255)
             aaCGAAATCGACGGTATCTCGCTCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGt                                                                                                                                     >  7:67125/1‑149 (MQ=255)
             aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGt                                                                                                                                     >  7:270871/1‑149 (MQ=255)
               cGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTc                                                                                                                                   <  4:58313/149‑1 (MQ=255)
                   aTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                               <  8:215220/149‑1 (MQ=255)
                    tCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGAATGCGGCCTTACTGGATCTTgtt                                                                                                                              >  7:23932/1‑149 (MQ=255)
                      gACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcg                                                                                                                            >  2:100624/1‑149 (MQ=255)
                       aCGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcgt                                                                                                                           >  7:294846/1‑149 (MQ=255)
                         ggTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcgttc                                                                                                                         <  1:100624/149‑1 (MQ=255)
                                cgcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                  >  7:148558/1‑149 (MQ=255)
                                  cgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCt                                                                                                                 >  1:298510/1‑148 (MQ=255)
                                       tACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTGGt                                                                                                                                     <  6:33707/123‑1 (MQ=255)
                                       tACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTGGt                                                                                                                                     >  5:33707/1‑123 (MQ=255)
                                          cTACGACTTGTCGAGCAAGCCGCCTGCCACTATCGATTGGGAGTGATCTGGTAGCCAGCCGATGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        <  8:197635/149‑1 (MQ=255)
                                          cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        <  7:350039/149‑1 (MQ=255)
                                          cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        >  2:15860/1‑149 (MQ=255)
                                          cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        <  8:67125/149‑1 (MQ=255)
                                             cgtcgTCTCGAGCAAGCCGCCTGCCACTCTCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTTGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTgg                                                                                                      <  8:109340/145‑1 (MQ=255)
                                                  tgtCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCt                                                                                                <  6:209754/149‑1 (MQ=255)
                                                    tCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTgg                                                                                              >  3:98540/1‑149 (MQ=255)
                                                      gAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAATCTTCTGGCTTTCTGGc                                                                                             >  1:96855/1‑148 (MQ=255)
                                                      gAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCAGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGc                                                                                             >  3:68462/1‑148 (MQ=255)
                                                       tgCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                               <  7:325951/112‑1 (MQ=255)
                                                        gCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                          <  2:298510/149‑1 (MQ=255)
                                                        gCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                          >  6:182896/1‑149 (MQ=255)
                                                         cAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAAt                                                                                         <  8:228917/149‑1 (MQ=255)
                                                          aaGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATAATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATa                                                                                        <  4:81987/149‑1 (MQ=255)
                                                           aGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAg                                                                                       >  8:165517/1‑149 (MQ=255)
                                                                ccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCtt                                                                                  >  7:160979/1‑149 (MQ=255)
                                                                 cTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCttt                                                                                 >  1:216235/1‑149 (MQ=255)
                                                                  tGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTc                                                                                >  4:236707/1‑149 (MQ=255)
                                                                    ccACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTg                                                                              <  5:22337/149‑1 (MQ=255)
                                                                      aCTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTgg                                                                             >  5:176962/1‑148 (MQ=255)
                                                                       cTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTgg                                                                             >  2:342587/1‑147 (MQ=255)
                                                                        tATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAg                                                                          >  3:365457/1‑149 (MQ=255)
                                                                          tCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGc                                                                        <  8:348729/149‑1 (MQ=255)
                                                                            gAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCg                                                                       <  4:365457/148‑1 (MQ=255)
                                                                            gAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                      <  4:98540/149‑1 (MQ=255)
                                                                              gTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGtt                                                                     <  8:148558/148‑1 (MQ=255)
                                                                                 ggAGTGATCGGGTAGCCAGCCCAGGCGTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTta                                                                    >  1:225670/1‑147 (MQ=255)
                                                                                 ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAg                                                                 <  6:168294/149‑1 (MQ=255)
                                                                                 ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAg                                                                 >  3:272059/1‑149 (MQ=255)
                                                                                   aGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAgg                                                                >  7:319437/1‑148 (MQ=255)
                                                                                     tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                            <  2:283662/149‑1 (MQ=255)
                                                                                     tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                            <  7:165517/149‑1 (MQ=255)
                                                                                     tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                            <  8:294846/149‑1 (MQ=255)
                                                                                      gATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAag                                                           >  4:397497/1‑149 (MQ=255)
                                                                                       aTCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAag                                                           >  8:160045/1‑148 (MQ=255)
                                                                                           ggTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGcc                                                      >  6:73494/1‑149 (MQ=255)
                                                                                            ggtgccagccTGGGCTGGCTACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTAAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGGTGGCGTTTTAGGATAGAAGcc                                                      >  6:300403/4‑149 (MQ=255)
                                                                                                    cccAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCa                                             >  8:5526/1‑149 (MQ=255)
                                                                                                         gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCttt                                                                                                  <  4:346054/92‑1 (MQ=255)
                                                                                                         gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCttt                                                                                                  >  3:346054/1‑92 (MQ=255)
                                                                                                         gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGc                                         >  3:271865/1‑148 (MQ=255)
                                                                                                          cTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                       <  1:15860/149‑1 (MQ=255)
                                                                                                              gCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGAt                                      <  8:270871/147‑1 (MQ=255)
                                                                                                                   ccGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTc                                <  4:272059/149‑1 (MQ=255)
                                                                                                                     gCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGt                              >  3:238114/1‑149 (MQ=255)
                                                                                                                     gCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCCTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGt                              >  6:160560/1‑149 (MQ=255)
                                                                                                                         gACGGATACACCAGAAAGCCCGCATGATTGCGTGCTTTCTGGTTCTTGTTCGTTCATATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAc                           <  4:234551/148‑1 (MQ=255)
                                                                                                                                  acCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgcg                 >  8:113358/1‑149 (MQ=255)
                                                                                                                                   ccAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGt                >  5:306581/1‑149 (MQ=255)
                                                                                                                                    cAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGtt               >  5:403894/1‑149 (MQ=255)
                                                                                                                                           cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTg        >  5:333131/1‑149 (MQ=255)
                                                                                                                                            ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGc       >  8:106654/1‑149 (MQ=255)
                                                                                                                                               cATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCttt    >  5:79133/1‑149 (MQ=255)
                                                                                                                                                 tGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGACAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtttt   <  1:342587/149‑1 (MQ=255)
                                                                                                                                                 tGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTc  <  3:397497/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTC  >  NC_002947/1178608‑1178899

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