Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,751 T→A 8.2% intergenic (+56/+38) 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,7510TA 8.2% 208.0 / 7.6 73intergenic (+56/+38)guaA/PP_1033glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (38/29);  new base A (2/4);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

TCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTC  >  NC_002947/1178604‑1178899
                                                                                                                                                     |                                                                                                                                                     
tCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAt                                                                                                                                                       <  5:121595/148‑1 (MQ=255)
tCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAt                                                                                                                                                       <  8:85761/148‑1 (MQ=255)
tCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAt                                                                                                                                                       <  4:121034/148‑1 (MQ=255)
    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                                                                                                                                     >  5:33707/1‑123 (MQ=255)
                                           tACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTGGt                                                                                                                                     <  6:33707/123‑1 (MQ=255)
                                              cTACGACTTGTCGAGCAAGCCGCCTGCCACTATCGATTGGGAGTGATCTGGTAGCCAGCCGATGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        <  8:197635/149‑1 (MQ=255)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        <  8:67125/149‑1 (MQ=255)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        <  7:350039/149‑1 (MQ=255)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                        >  2:15860/1‑149 (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                                                            <  8:294846/149‑1 (MQ=255)
                                                                                         tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                            <  2:283662/149‑1 (MQ=255)
                                                                                         tGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                            <  7:165517/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                                                                                                  >  3:346054/1‑92 (MQ=255)
                                                                                                             gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCttt                                                                                                  <  4:346054/92‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
TCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTTC  >  NC_002947/1178604‑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: