Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,252,594 T→A 6.8% intergenic (+81/+162) ureG → / ← PP_2851 urease accessory protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,252,5940TA 6.8% 176.2 / 5.0 59intergenic (+81/+162)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/21);  new base A (2/2);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

GGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGC  >  NC_002947/3252448‑3252743
                                                                                                                                                  |                                                                                                                                                     
ggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTAc                                                                                                                                                    >  1:867814/1‑150 (MQ=255)
      gggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGAGCCGCGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGGGa                                                                                                                                              >  1:977953/1‑150 (MQ=255)
       ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                                                                                             >  1:1001822/1‑150 (MQ=255)
       ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                                                                             <  1:729574/150‑1 (MQ=255)
             aGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                                >  1:7132/1‑141 (MQ=255)
             aGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                                <  2:7132/141‑1 (MQ=255)
                ccGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                                                                                    <  2:867814/150‑1 (MQ=255)
                   tGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGc                                                                                                                                 >  1:1238764/1‑150 (MQ=255)
                    gATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGAACAGAGACCAACGCGGTGGCTGGAACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGcc                                                                                                                                >  1:1095193/1‑150 (MQ=255)
                        gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGTGGGCATGTCCGCACAGAGACAAACGCGGTGGCTGGCACCGGCTCTGCCGGTGATCGCGGGTTGACCCGCGAATACGGTAGCGGCGGcaa                                                                                                                            >  2:163165/1‑150 (MQ=255)
                            ttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGTTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                                                                                        <  1:92647/150‑1 (MQ=255)
                            ttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                                                                                        <  2:381524/150‑1 (MQ=255)
                             tATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACACTg                                                                                                                       >  2:306103/1‑150 (MQ=255)
                              aTCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGa                                                                                                                      >  1:271177/1‑150 (MQ=255)
                                      gcgcGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCggg                                                                                                              <  1:936730/150‑1 (MQ=255)
                                        gcgGGATGCTGGGGGTGCTAGGCTGACAGATCCAGGACTCTGTGCGAGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTCTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTg                                                                                                            <  1:163165/150‑1 (MQ=255)
                                         cgGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCAGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTgg                                                                                                           >  1:919897/1‑150 (MQ=255)
                                          gggATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGa                                                                                                          <  2:67952/150‑1 (MQ=255)
                                                cTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTgg                                                                                                    >  1:14813/1‑150 (MQ=255)
                                                 tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                                                                                   <  2:619398/150‑1 (MQ=255)
                                                   ggggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCca                                                                                                 >  1:238051/1‑150 (MQ=255)
                                                              cTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccggc                                                                                      >  2:760455/1‑150 (MQ=255)
                                                               tGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcc                                                                                                  <  1:638256/137‑1 (MQ=255)
                                                               tGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcc                                                                                                  >  2:638256/1‑137 (MQ=255)
                                                                 aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATg                                                                                                                       <  2:563965/114‑1 (MQ=255)
                                                                 aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATg                                                                                                                       >  1:563965/1‑114 (MQ=255)
                                                                      tCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAACAGGGCCGGCCTCTTcgc                                                                              >  2:256652/1‑150 (MQ=255)
                                                                             ctctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTaaa                                                                       >  1:466433/1‑150 (MQ=255)
                                                                              tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                                                                      <  2:238051/150‑1 (MQ=255)
                                                                              tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                                                                      <  1:115637/150‑1 (MQ=255)
                                                                              tctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                                                                      >  1:740039/1‑150 (MQ=255)
                                                                                tgtgCGGGCATGTCCGCACAGAGCCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaa                                                                                                                                             <  1:996962/77‑1 (MQ=18)
                                                                                tgtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaa                                                                                                                                             >  2:996962/1‑77 (MQ=25)
                                                                                            tCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACgggg                                                        >  2:926913/1‑150 (MQ=255)
                                                                                               gCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGa                                                     >  1:589401/1‑150 (MQ=255)
                                                                                                   agagACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACcgc                                                 <  2:589401/150‑1 (MQ=255)
                                                                                                    gagaCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGGTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACcgcg                                                >  1:1078140/1‑150 (MQ=255)
                                                                                                       aCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCAGCAACAATTAATTGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAg                                             <  2:368691/150‑1 (MQ=255)
                                                                                                         cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc                                                                  >  1:1138407/1‑127 (MQ=255)
                                                                                                         cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc                                                                  <  2:1138407/127‑1 (MQ=255)
                                                                                                         cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGc                                           >  2:227232/1‑150 (MQ=255)
                                                                                                         cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGc                                           >  1:29379/1‑150 (MQ=255)
                                                                                                          aaCGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGcc                                          >  2:937563/1‑150 (MQ=255)
                                                                                                          aaCGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGcc                                          >  1:880080/1‑150 (MQ=255)
                                                                                                             gcgGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTc                                       <  2:271177/150‑1 (MQ=255)
                                                                                                                 tggctggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCAGGTTATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAAt                                   <  1:12798/150‑1 (MQ=255)
                                                                                                                    ctggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGtt                                <  2:919897/150‑1 (MQ=255)
                                                                                                                    ctggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGtt                                <  1:760455/150‑1 (MQ=255)
                                                                                                                    ctggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGtt                                <  1:497165/150‑1 (MQ=255)
                                                                                                                     tggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTg                               >  1:99400/1‑150 (MQ=255)
                                                                                                                          ccGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAg                          >  1:31868/1‑150 (MQ=255)
                                                                                                                           cGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGATGACCTGGCAGGCCGTCGAATGTTGAGCAgg                         >  1:643215/1‑150 (MQ=255)
                                                                                                                              ctctGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACa                      >  1:624360/1‑150 (MQ=255)
                                                                                                                               tctGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAg                     >  1:778199/1‑150 (MQ=255)
                                                                                                                                     ggTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGccc               >  1:422422/1‑150 (MQ=255)
                                                                                                                                     ggTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGccc               <  2:466433/150‑1 (MQ=255)
                                                                                                                                      gtattCGCGGGTAAACCCGCGAACACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGcccc              >  2:907863/4‑150 (MQ=25)
                                                                                                                                            gcgGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAGATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAg        >  1:93003/1‑150 (MQ=255)
                                                                                                                                            gcgGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAg        >  2:177566/1‑150 (MQ=255)
                                                                                                                                                  ttACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGc  <  2:624360/150‑1 (MQ=255)
                                                                                                                                                  ttACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGc  <  2:1212180/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGC  >  NC_002947/3252448‑3252743

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