Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,758 G→C 5.5% 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,7580GC 5.5% 293.1 / 6.7 91intergenic (+57/+37)guaA/PP_1033glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (44/42);  new base C (2/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

TCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTT  >  NC_002947/1178610‑1178904
                                                                                                                                                      |                                                                                                                                                  
tCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGa                                                                                                                                                   >  1:703945/1‑150 (MQ=255)
   gCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTg                                                                                                                                                <  2:1931461/150‑1 (MQ=255)
    cGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGc                                                                                                                                               >  2:778571/1‑150 (MQ=255)
     ggTCGTATCATCCACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCg                                                                                                                                              <  1:677981/150‑1 (MQ=255)
       tCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCggg                                                                                                                                            >  2:1938246/1‑150 (MQ=255)
         gTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCt                                                                                                                                          >  1:1159788/1‑150 (MQ=255)
           atcatcAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCA‑TGCGGGCttt                                                                                                                                        <  1:1303313/150‑1 (MQ=255)
           atcatcAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCgg                                                                                                                                             <  1:1044217/145‑1 (MQ=255)
           atcatcAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCgg                                                                                                                                             >  2:1044217/1‑145 (MQ=255)
           atcatcAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCttt                                                                                                                                        >  2:704124/1‑150 (MQ=255)
            tcatcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGTTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTc                                                                                                                                       <  1:2163177/150‑1 (MQ=255)
             catcaACGAAAGCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGAGCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCt                                                                                                                                      <  1:103974/150‑1 (MQ=255)
              atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGCATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTg                                                                                                                                     >  1:391206/1‑150 (MQ=255)
                caACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCA‑TGCGGGCTTTCTGGt                                                                                                                                   <  1:505926/150‑1 (MQ=255)
                 aaCGAAACCGACGGTATCCGGCGCGTTACCGACGAGGTGGCGAGCAAGGCGCGTGGGAGTGGCGAG‑TGGGAGTGTGCGGGTAGCCAGCCCAGGCTGGCACCCGAAGGGCGGATACACCAGGAAGTCCGC‑ATGATTGCGGGCTTTCTGGtt                                                                                                                                  <  2:749054/150‑1 (MQ=255)
                   cGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCt                                                                                                                                >  1:2048338/1‑150 (MQ=255)
                      aaTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                             >  1:1942942/1‑150 (MQ=255)
                      aaTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                             >  2:1059261/1‑150 (MQ=255)
                          gACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgttcgt                                                                                                                         >  2:2020534/1‑150 (MQ=255)
                          gACGGTATCTCGCGCGTTACCTACGACGTGCCGAGCAAGCCGCCTGCCACTA‑CAAGCAGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgttcgt                                                                                                                         <  1:1194321/150‑1 (MQ=255)
                           aCGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgttcgtt                                                                                                                        <  1:2094292/150‑1 (MQ=255)
                                 tctcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGAtt                                                                                                                  >  2:1740484/1‑150 (MQ=255)
                                   tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                             >  2:104795/1‑137 (MQ=255)
                                   tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgt                                                                                                                             <  1:104795/137‑1 (MQ=255)
                                   ccgcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATTGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                                >  2:587044/2‑150 (MQ=255)
                                     gcgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTATAGTTCGGTAAGATTTCtt                                                                                                              >  2:1886071/1‑150 (MQ=255)
                                      cgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                             <  2:1159788/150‑1 (MQ=255)
                                         gTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgttcgttc                                                                                                                       <  2:465514/137‑1 (MQ=255)
                                         gTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTgttcgttc                                                                                                                       >  1:465514/1‑137 (MQ=255)
                                            aCCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAttcttc                                                                                                       <  1:2117281/150‑1 (MQ=255)
                                             ccTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                      <  1:778571/150‑1 (MQ=255)
                                             ccTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                      <  2:2048338/150‑1 (MQ=255)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCA‑TGCGGGCTTTCTGGt                                                                                                                                   >  1:58409/1‑120 (MQ=38)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCA‑TGCGGGCTTTCTGGt                                                                                                                                   <  2:58409/120‑1 (MQ=38)
                                               tACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTgg                                                                                                    >  2:1097899/1‑150 (MQ=255)
                                               tACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTgg                                                                                                    <  1:1097899/150‑1 (MQ=255)
                                                   acgTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTc                                                                                                                                 >  1:902566/1‑117 (MQ=255)
                                                   acgTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCGGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTc                                                                                                                                 <  2:902566/117‑1 (MQ=255)
                                                          gAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTg                                                                                                                                     <  2:813322/106‑1 (MQ=255)
                                                          gAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTg                                                                                                                                     >  1:813322/1‑106 (MQ=255)
                                                          gAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaa                                                                                         >  1:1271161/1‑150 (MQ=255)
                                                               aGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATa                                                                                      <  1:917779/148‑1 (MQ=255)
                                                               aGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATa                                                                                      >  2:917779/1‑148 (MQ=255)
                                                                  cgccTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCt                                                                                 >  1:2007505/1‑150 (MQ=255)
                                                                    ccTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCttt                                                                               <  1:1059261/150‑1 (MQ=255)
                                                                     cTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTc                                                                              <  1:2020534/150‑1 (MQ=255)
                                                                      tGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATGGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCt                                                                             >  1:187174/1‑150 (MQ=255)
                                                                       gCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTg                                                                            <  2:1942942/150‑1 (MQ=255)
                                                                        ccACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTg                                                                                                                              >  1:1060661/1‑99 (MQ=255)
                                                                        ccACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTg                                                                                                                              <  2:1060661/99‑1 (MQ=255)
                                                                        ccACTATCGAG‑TGGGAGTGATCGGGTAACGATCCGAAGGGCGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTgg                                                                           <  2:2007505/150‑1 (MQ=255)
                                                                          aCTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTa                                                                         >  1:1915500/1‑150 (MQ=255)
                                                                                gAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTgg                                                                                                    >  2:888736/1‑117 (MQ=255)
                                                                                gAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTgg                                                                                                    <  1:888736/117‑1 (MQ=255)
                                                                                  g‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtat                                                                 >  1:492588/1‑150 (MQ=255)
                                                                                     gggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATAAACCAGCAAGCCCGC‑ATGATTGAGGGCTATCTGGTTCTTGTTAGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTGTAg                                                               >  2:1007719/1‑150 (MQ=255)
                                                                                       gAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGa                                                             >  2:484068/1‑150 (MQ=255)
                                                                                       gAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACAACAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAgaa                                                             >  2:863533/1‑148 (MQ=255)
                                                                                         gTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                           <  1:1462547/150‑1 (MQ=255)
                                                                                           gATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                    >  2:708151/1‑139 (MQ=255)
                                                                                           gATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                    <  1:708151/139‑1 (MQ=255)
                                                                                                    gccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                    <  1:773390/130‑1 (MQ=255)
                                                                                                    gccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                    >  2:773390/1‑130 (MQ=255)
                                                                                                     ccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagca                                                        >  1:345489/1‑141 (MQ=255)
                                                                                                     ccagccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagca                                                        <  2:345489/141‑1 (MQ=255)
                                                                                                        gccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCa                                            >  1:2002789/1‑150 (MQ=255)
                                                                                                            aGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGc                                        <  2:703945/150‑1 (MQ=255)
                                                                                                              gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                      <  1:704124/150‑1 (MQ=255)
                                                                                                                  gCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGaa                                  <  1:587044/150‑1 (MQ=255)
                                                                                                                    aCCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTgg                                                                           <  2:1605537/107‑1 (MQ=255)
                                                                                                                    aCCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTgg                                                                           >  1:1605537/1‑107 (MQ=255)
                                                                                                                      ccGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                                >  2:210397/1‑132 (MQ=255)
                                                                                                                      ccGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                                <  1:210397/132‑1 (MQ=255)
                                                                                                                         cAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCa                           <  2:2002789/150‑1 (MQ=255)
                                                                                                                          aGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAgg                                                                       <  1:687119/105‑1 (MQ=255)
                                                                                                                          aGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAgg                                                                       >  2:687119/1‑105 (MQ=255)
                                                                                                                              cGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAACGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTc                      <  1:1886071/150‑1 (MQ=255)
                                                                                                                                   acacCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                           >  1:494902/1‑108 (MQ=255)
                                                                                                                                   acacCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                           <  2:494902/108‑1 (MQ=255)
                                                                                                                                      ccAGAAAGCCCGCAATCA‑TGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCT‑GTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTg             >  1:1948786/1‑150 (MQ=255)
                                                                                                                                        aGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg            >  1:233276/1‑150 (MQ=255)
                                                                                                                                        aGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgg            <  2:1915500/150‑1 (MQ=255)
                                                                                                                                          aaaGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg          <  1:1740484/150‑1 (MQ=255)
                                                                                                                                          aaaGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTgggg          >  2:690479/1‑150 (MQ=255)
                                                                                                                                               ccGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCt     >  2:181546/1‑150 (MQ=255)
                                                                                                                                                cGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtt    >  2:762427/1‑150 (MQ=255)
                                                                                                                                                cGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtt    >  2:1711063/1‑150 (MQ=255)
                                                                                                                                                cGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGGAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGGGGAACTCGTCACCCTCCACGCGTTGGGGGTGCtt    <  1:1938246/150‑1 (MQ=255)
                                                                                                                                                  c‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtttt  <  2:233276/150‑1 (MQ=255)
                                                                                                                                                     tGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGt        >  2:1207574/1‑142 (MQ=255)
                                                                                                                                                     tGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGt        <  1:1207574/142‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAG‑TGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGC‑ATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTTTT  >  NC_002947/1178610‑1178904

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: ‑