Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,703,552 A→T 5.4% L553H (CTT→CAT)  kdpA ← potassium‑transporting ATPase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,703,5520AT 5.4% 227.7 / 3.5 74L553H (CTT→CAT) kdpApotassium‑transporting ATPase subunit A
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base T (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GACGAACGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGTATC  >  NC_002947/4703407‑4703701
                                                                                                                                                  |                                                                                                                                                     
gACGAACGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGc                                                                                                                                                   >  2:588583/1‑150 (MQ=255)
  cGAACGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGCTCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCa                                                                                                                                                 <  1:252466/150‑1 (MQ=255)
   gAACGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAg                                                                                                                                                <  1:991113/150‑1 (MQ=255)
     aCGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCTTGGGCCaatc                                                                                                                                              <  2:28451/150‑4 (MQ=255)
     aCGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCTTGGGCCaat                                                                                                                                               >  2:1330457/1‑147 (MQ=255)
     aCGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCTTGGGCCaat                                                                                                                                               <  1:1330457/149‑3 (MQ=255)
            cACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAg                                                                                                                                       <  2:493703/150‑1 (MQ=255)
               cAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAg                                                                                                                                                <  1:854728/138‑1 (MQ=255)
               cAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAg                                                                                                                                                >  2:854728/1‑138 (MQ=255)
               cAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGaa                                                                                                                                    <  2:1095692/150‑1 (MQ=255)
               cAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGGCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCTTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGaa                                                                                                                                    <  1:1384487/150‑1 (MQ=255)
                aGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAg                                                                                                                                   >  2:230241/1‑150 (MQ=255)
                  cgcCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGt                                                                                                                                 >  2:815763/1‑150 (MQ=255)
                     cGGGGGCCACCGCGCGGCGCAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAg                                                                                                                              <  2:307829/150‑1 (MQ=255)
                     cGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAg                                                                                                                              <  2:639019/150‑1 (MQ=255)
                          gcCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCgg                                                                                                                                          >  2:383298/1‑133 (MQ=255)
                          gcCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCgg                                                                                                                                          <  1:383298/133‑1 (MQ=255)
                          gcCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCAc                                                                                                                         <  1:815763/150‑1 (MQ=255)
                            cacaGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCg                                                                                                                       >  2:1285829/1‑150 (MQ=255)
                              caGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTTACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCCCAATCTTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGAc                                                                                                                     <  1:1170431/150‑1 (MQ=255)
                                 cgcgGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACcag                                                                                                                  >  2:1146945/1‑150 (MQ=255)
                                   cggcgAAGCGGGTCTGGTCCTTGGCAACGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACcagca                                                                                                                >  1:616464/1‑150 (MQ=255)
                                    ggcgAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACcagcag                                                                                                               <  2:1424227/150‑1 (MQ=255)
                                     gcgAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCAAAGACcaggcgg                                                                                                              >  2:865238/1‑146 (MQ=255)
                                              gTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCAcc                                                                                                     >  2:466194/1‑150 (MQ=255)
                                               tCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACca                                                                                                    >  1:1258149/1‑150 (MQ=255)
                                                    tCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGt                                                                                                         <  1:960850/140‑1 (MQ=255)
                                                    tCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGt                                                                                                         >  2:960850/1‑140 (MQ=255)
                                                       ttGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAg                                                                                            >  1:416306/1‑150 (MQ=255)
                                                       ttGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAg                                                                                            >  2:607619/1‑150 (MQ=255)
                                                           cACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACca                                                                                                    >  1:583805/1‑138 (MQ=255)
                                                           cACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACca                                                                                                    <  2:583805/138‑1 (MQ=255)
                                                            aCTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCCCGTACCAGCAGGATGGTCACCAGCAGCAAGCCGGt                                                                                       <  1:275203/150‑1 (MQ=255)
                                                                  ccGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCg                                                                                         <  1:892000/142‑1 (MQ=255)
                                                                  ccGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCg                                                                                         >  2:892000/1‑142 (MQ=255)
                                                                    ggCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCg                                                                               <  2:1258149/150‑1 (MQ=255)
                                                                    ggCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCg                                                                               <  1:466194/150‑1 (MQ=255)
                                                                     gCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCgg                                                                              <  1:607619/150‑1 (MQ=255)
                                                                     gCTATTCACTTCAGGAATGGGCATGTCCATGATGGTTCATCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGCTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGCCACCAGCAGCAAGGCGGTGAACAGCgg                                                                              >  2:100810/1‑150 (MQ=255)
                                                                        aTTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGcc                                                                           >  1:1158875/1‑150 (MQ=255)
                                                                        aTTAACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCATGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGcc                                                                           <  1:76900/150‑1 (MQ=255)
                                                                          tCACTTCAGGAATGCGCATGTTACCATCCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGt                                                                         <  1:714200/150‑1 (MQ=255)
                                                                            aCTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGAGGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGa                                                                       >  1:373394/1‑150 (MQ=255)
                                                                               tCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTc                                                                    <  1:1404717/150‑1 (MQ=255)
                                                                                   gaagtggCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGa                                                                >  2:356376/6‑150 (MQ=255)
                                                                                      tGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGc                                                             <  1:1209952/150‑1 (MQ=255)
                                                                                         gCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACcag                                                                                                                  >  2:1138918/1‑94 (MQ=255)
                                                                                         gCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACcag                                                                                                                  <  1:1138918/94‑1 (MQ=255)
                                                                                          cATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGtt                                                         <  2:1283881/150‑1 (MQ=255)
                                                                                           aTGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTc                                                        >  1:70222/1‑150 (MQ=255)
                                                                                              ttCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATg                                                                                                           <  2:972250/96‑1 (MQ=255)
                                                                                              ttCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATg                                                                                                           >  1:972250/1‑96 (MQ=255)
                                                                                                 aTGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGTGTCGGCAGGAAg                                                                                                                                   <  2:845328/69‑1 (MQ=255)
                                                                                                 aTGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGTGTCGGCAGGAAg                                                                                                                                   >  1:845328/1‑69 (MQ=255)
                                                                                                       cGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGAGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAgcgg                                             >  1:629296/1‑150 (MQ=255)
                                                                                                          tCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGGAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGgcgc                                          >  1:1184678/1‑150 (MQ=255)
                                                                                                          tCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGgcgc                                          >  1:903082/1‑150 (MQ=255)
                                                                                                          tCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGgcgc                                          >  2:1237314/1‑150 (MQ=255)
                                                                                                            cTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGgcgcgc                                        >  1:934977/1‑150 (MQ=255)
                                                                                                             tCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGCCACATGGGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCt                                       >  1:1384270/1‑150 (MQ=255)
                                                                                                                 acaCCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGACGGCCAGGCCTACCTTCATCAGTAACGTTACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCttttt                                   <  1:100810/148‑1 (MQ=255)
                                                                                                                  aaCCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCtttttt                                  >  1:1412998/1‑150 (MQ=255)
                                                                                                                       aaGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTcgcc                             <  1:1146945/150‑1 (MQ=255)
                                                                                                                          cTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCg                                                                                         >  2:330234/1‑87 (MQ=255)
                                                                                                                          cTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCg                                                                                         <  1:330234/87‑1 (MQ=255)
                                                                                                                           tAAGGTGTTCGGCGCTTGCCCACCCGGCCAGGCGCGGTAGGGATGTTATGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCa                         <  2:373394/150‑1 (MQ=255)
                                                                                                                           tAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTc                                <  1:1413294/143‑1 (MQ=255)
                                                                                                                           tAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTc                                >  2:1413294/1‑143 (MQ=255)
                                                                                                                                 gtTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCttc                   <  2:1417888/150‑3 (MQ=255)
                                                                                                                                         cctGGCCCGTGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTGTTTTTCGCCGCCAGGCTGCCGGCCAAc           >  1:889531/3‑150 (MQ=255)
                                                                                                                                             gCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCCACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCa       <  1:802536/150‑1 (MQ=255)
                                                                                                                                               ccAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGt     >  1:306535/1‑150 (MQ=255)
                                                                                                                                               ccAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGt     >  1:879838/1‑150 (MQ=255)
                                                                                                                                               ccAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGCTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGt     >  1:1291171/1‑150 (MQ=255)
                                                                                                                                                 aaGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGTAt   >  1:960937/1‑150 (MQ=255)
                                                                                                                                                 aaGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGTAt   <  2:903082/150‑1 (MQ=255)
                                                                                                                                                  aGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGTATc  <  1:356376/150‑1 (MQ=255)
                                                                                                                                                  aGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGTATc  <  1:1285829/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GACGAACGCCTGCACCAGCGCCGGGCGCCACAGCGCGGCGAAGCGGGTCTGGTCCTTGGCACTGTGCCGGCTATTCACTTCAGGAATGGGCATGTTCATGATCGGTCCTCAAAAACCCAAGCTAAGGTGTTCGGCGATTGGCCCAAGGGCCAGGGTCGGCAGGAAGGTCAGGCCACCGACCAGCAGGATGGTCACCAGCAGCAAGCCGGTGAACAGCGGGCCGTGAGTCGGGAAGCTGTTCTGCCCCAGCGGCGCGCTTTTTTTCGCCGCCAGGCTGCCGGCCAACGCCAGTATC  >  NC_002947/4703407‑4703701

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