Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,912,177 C→T 6.0% P160L (CCC→CTC)  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,1770CT 6.0% 235.3 / 3.4 82P160L (CCC→CTC) PP_3452GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base T (3/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912033‑3912322
                                                                                                                                                |                                                                                                                                                 
cccTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATg                                                                                                                                              <  2:961600/150‑1 (MQ=255)
     cgcCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCCTGATCAGCACCTGGGCGGTTCCTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTc                                                                                                                                         <  2:1901936/150‑1 (MQ=255)
        ccGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGACCGATGCGCTCGaa                                                                                                                                      >  1:1711912/1‑150 (MQ=255)
                gttggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCAt                                                                                                                              >  1:797612/1‑150 (MQ=255)
                gttggtGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCAt                                                                                                                              >  2:1147909/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                                                                                            >  1:849616/1‑123 (MQ=255)
                                                                             tCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTg                                                                                            <  2:849616/123‑1 (MQ=255)
                                                                               aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  2:942750/1‑136 (MQ=255)
                                                                               aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:1846157/1‑136 (MQ=255)
                                                                               aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  1:942750/136‑1 (MQ=255)
                                                                               aTTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             <  2:1846157/136‑1 (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                                                                                                              <  2:375489/69‑1 (MQ=255)
                                                                                                                 cTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGcc                                                                                                              >  1:375489/1‑69 (MQ=255)
                                                                                                                    aaTCCGTTCTACCACACCCACATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTg                          <  1:1717643/150‑1 (MQ=255)
                                                                                                                         gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctac                                                                                   >  2:417914/1‑88 (MQ=255)
                                                                                                                         gTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTcagctac                                                                                   <  1:417914/88‑1 (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                                                                             <  2:188148/81‑1 (MQ=255)
                                                                                                                                    agcGCATCGGGCTCGATGCGCTAGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCagctacagctac                                                                             >  1:188148/3‑83 (MQ=255)
                                                                                                                                     gcGCATCGGGCTCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGt         >  2:1162978/2‑150 (MQ=255)
                                                                                                                                     ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGc                                                                                               >  1:1420234/1‑64 (MQ=255)
                                                                                                                                     ccGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGc                                                                                               <  2:1420234/64‑1 (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)
                                                                                                                                                |                                                                                                                                                 
CCCTGCGCCCGGTGTGGTTGGTGCTCGGGGTGTCCACCGTGCTGTATGCCATGTCGGTGCTACCGGTGCGGCGCTGGTCATTCTTCTGCGCCATGGTCATCACCTGGGTGATGCTTAATCCGTTCTACCACACCCGCATCGAGCCCGATGCGCTCGAAGGCGCCATGCTGATCAGCTATGCCGTGTTTCTAAGCGGCCTGGTCAGCTACAGCTACCTGCAGATCCGCAAGGCCAAGCTGCACAATTTCTACCTGTCCAAGGTGTTGCTGGAGCAAGCCTACGTCGATGCC  >  NC_002947/3912033‑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: