Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,174 A→G 6.7% E159G (GAG→GGG)  PP_3452 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,912,1740AG 6.7% 231.4 / 4.4 85E159G (GAG→GGG) PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base G (3/3);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

CCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912026‑3912322
                                                                                                                                                    |                                                                                                                                                    
ccGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAg                                                                                                                                                     <  1:1911842/150‑1 (MQ=255)
ccGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACACGCATCGAg                                                                                                                                                     >  1:928204/1‑150 (MQ=255)
ccGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGGGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAg                                                                                                                                                     >  1:166919/1‑150 (MQ=255)
       cccTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATg                                                                                                                                              <  2:961600/150‑1 (MQ=255)
            cgcCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCCTGATCAGCACCTGGGCGGTTCCTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTc                                                                                                                                         <  2:1901936/150‑1 (MQ=255)
               ccGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGACCGATGCGCTCGaa                                                                                                                                      >  1:1711912/1‑150 (MQ=255)
                       gttggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCAt                                                                                                                              >  2:1147909/1‑150 (MQ=255)
                       gttggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCAt                                                                                                                              >  1:797612/1‑150 (MQ=255)
                        ttggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAg                                                                                                                                     >  1:987344/1‑142 (MQ=255)
                        ttggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCACCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAg                                                                                                                                     <  2:987344/142‑1 (MQ=255)
                         tggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGc                                                                                                                            <  2:2380932/150‑1 (MQ=255)
                          ggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCt                                                                                                                           <  1:2208696/150‑1 (MQ=255)
                           gtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTg                                                                                                                          <  2:166919/150‑1 (MQ=255)
                            tGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGa                                                                                                                         >  1:2031851/1‑150 (MQ=255)
                            tGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGa                                                                                                                         >  1:1876424/1‑150 (MQ=255)
                               tCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCa                                                                                                                      >  1:1219007/1‑150 (MQ=255)
                                     tgtCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATg                                                                                                                <  1:519930/150‑1 (MQ=255)
                                         cACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGGGCTCGATgcgggtg                                                                                                                                        <  2:2229464/122‑5 (MQ=255)
                                         cACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGGGCTCGATgcgggtg                                                                                                                                        >  1:2229464/1‑118 (MQ=255)
                                         cACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATGCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCGGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCgt                                                                                                            >  2:1231821/1‑150 (MQ=255)
                                           ccgcgCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCgtgt                                                                                                          <  1:763448/146‑1 (MQ=255)
                                           ccGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCgtgt                                                                                                          >  1:1554908/1‑150 (MQ=255)
                                            cGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGAGTCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGtt                                                                                                         >  1:1707125/1‑150 (MQ=255)
                                              tGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTc                                                                                                       >  1:193124/1‑150 (MQ=255)
                                                   tATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGc                                                                                                  <  2:1043381/150‑1 (MQ=255)
                                                      gCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGc                                                                                               <  2:2031851/150‑1 (MQ=255)
                                                              ggTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcag                                                                                       >  2:872969/1‑150 (MQ=255)
                                                              aattcTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcag                                                                                       >  1:683908/5‑150 (MQ=255)
                                                                     ccGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTGAGCg                                                                                                 <  2:1020542/133‑1 (MQ=255)
                                                                     ccGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCg                                                                                                 >  1:1020542/1‑133 (MQ=255)
                                                                     ccGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctacagc                                                                                >  2:1085989/1‑150 (MQ=255)
                                                                       ggTGCTGCCCGGGTCATTCTAATGCAACATGGTCATCACCTGGGGGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctacagcta                                                                              >  2:998803/1‑150 (MQ=255)
                                                                        gTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  1:1147909/150‑1 (MQ=255)
                                                                           cggcgCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTg                                                                          >  1:1104424/1‑150 (MQ=255)
                                                                             gcgcTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCa                                                                        <  2:2197551/150‑1 (MQ=255)
                                                                                   gTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCa                                                                        <  1:347098/144‑1 (MQ=255)
                                                                                   gTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCa                                                                        >  2:347098/1‑144 (MQ=255)
                                                                                   gTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCg                                                                  >  1:210494/1‑150 (MQ=255)
                                                                                   gTCATTCTTCTGCACCATGGTCATCACCTGGGCGATGCTTACTCCATTCCGCCAGACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACCGCTACCTGCAGATCCg                                                                  >  1:1006242/1‑150 (MQ=255)
                                                                                    tCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            <  2:849616/123‑1 (MQ=255)
                                                                                    tCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            >  1:849616/1‑123 (MQ=255)
                                                                                      aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:1846157/136‑1 (MQ=255)
                                                                                      aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  2:942750/1‑136 (MQ=255)
                                                                                      aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  1:942750/136‑1 (MQ=255)
                                                                                      aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:1846157/1‑136 (MQ=255)
                                                                                      aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCaa                                                               <  1:1020307/150‑1 (MQ=255)
                                                                                       ttcttcTGCGCCATGGTCATCTCCTGGGTGACGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAg                                                              <  2:27736/150‑1 (MQ=255)
                                                                                       ttcttcTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAg                                                              >  1:727415/1‑150 (MQ=255)
                                                                                           tctGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCa                                                          >  2:1537075/1‑150 (MQ=255)
                                                                                             tGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCTCACACGCATCGGGCCCGCTGCGCTCGAAGGCCCCATGCTGATCAGCTCTGCCTTCCTTCTATACGGCCTGGTCAGCTACTGCTACCTGCAGGTCCGCAAGGCCAAg                                                        <  2:928204/150‑1 (MQ=255)
                                                                                              gcgcCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCACCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAga                                                       >  2:341227/1‑149 (MQ=255)
                                                                                                  cATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCACATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:1647541/1‑124 (MQ=255)
                                                                                                  cATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCACATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:1647541/124‑1 (MQ=255)
                                                                                                  cATGGTCATCACCTGGGTGATGCTTAATCCGTGCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGca                                                   >  2:1348363/1‑150 (MQ=255)
                                                                                                  attaagcatcaCCCAGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGca                                                   >  1:869762/7‑150 (MQ=255)
                                                                                                      gTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAAt                                               >  1:2258302/1‑150 (MQ=255)
                                                                                                       tcatcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            <  1:2248764/104‑1 (MQ=255)
                                                                                                       tcatcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            >  2:2248764/1‑104 (MQ=255)
                                                                                                          tcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCt                                           <  2:797612/150‑1 (MQ=255)
                                                                                                          tcaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCt                                           >  1:2422225/1‑150 (MQ=255)
                                                                                                           gaCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGCGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTa                                          <  2:1707125/149‑1 (MQ=255)
                                                                                                            aCCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGACCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTAc                                         <  2:1711912/150‑1 (MQ=255)
                                                                                                                    gATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCaa                                 <  1:987571/150‑1 (MQ=255)
                                                                                                                     aTGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGGGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCacg                                >  1:2269593/1‑148 (MQ=255)
                                                                                                                       gCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGgt                              >  1:2295512/1‑150 (MQ=255)
                                                                                                                        cTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGcc                                                                                                              >  1:375489/1‑69 (MQ=255)
                                                                                                                        cTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGcc                                                                                                              <  2:375489/69‑1 (MQ=255)
                                                                                                                           aaTCCGTTCTACCACACCCACATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTg                          <  1:1717643/150‑1 (MQ=255)
                                                                                                                                gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctac                                                                                   <  1:417914/88‑1 (MQ=255)
                                                                                                                                gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctac                                                                                   >  2:417914/1‑88 (MQ=255)
                                                                                                                                  tCTACCCCCCCCCCAGCGAGCCCGATGCGCTCGAAGGCGCGAGGCGGATGAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGc                   <  2:1006242/150‑1 (MQ=255)
                                                                                                                                    tACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCaa                 <  2:869762/150‑1 (MQ=255)
                                                                                                                                       cacaCCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGGGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGcc              >  1:170118/1‑150 (MQ=255)
                                                                                                                                        acacCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCt             >  1:2025875/1‑150 (MQ=255)
                                                                                                                                           agcGCATCGGGCTCGATGCGCTAGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:188148/3‑83 (MQ=255)
                                                                                                                                           agcGCATCGGGCTCGATGCGCTAGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:188148/81‑1 (MQ=255)
                                                                                                                                            gcGCATCGGGCTCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  2:1162978/2‑150 (MQ=255)
                                                                                                                                            ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGc                                                                                               <  2:1420234/64‑1 (MQ=255)
                                                                                                                                            ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGc                                                                                               >  1:1420234/1‑64 (MQ=255)
                                                                                                                                            ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         <  1:1085989/150‑1 (MQ=255)
                                                                                                                                               cATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGa      <  1:2339429/150‑1 (MQ=255)
                                                                                                                                                aTCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTATGTCGAt     <  2:727415/150‑1 (MQ=255)
                                                                                                                                                 tCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATg    >  2:868300/1‑150 (MQ=255)
                                                                                                                                                  cGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGc   <  2:2422225/150‑1 (MQ=255)
                                                                                                                                                   gAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGcc  <  2:1876424/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGATAGCCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912026‑3912322

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: