Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,015,997 A→T 5.8% S437C (AGC→TGC)  PP_3543 → iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,015,9970AT 5.8% 206.6 / 3.4 69S437C (AGC→TGC) PP_3543iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base T (2/2);  total (31/38)
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.06e-01

ACAAGACCGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGAAGGCCAGAAT  >  NC_002947/4015848‑4016144
                                                                                                                                                     |                                                                                                                                                   
acaaGACCGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGagcagca                                                                                                                                                     <  2:1006788/150‑1 (MQ=255)
 caaGACCGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGagcagcag                                                                                                                                                    >  1:1316723/1‑150 (MQ=255)
  aaGACCGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGagcagcagc                                                                                                                                                   >  1:823664/1‑150 (MQ=255)
  aaGACCGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGagcagcagc                                                                                                                                                   >  1:1523143/1‑150 (MQ=255)
      ccGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCTGCTGCTCgg                                                                                                                                               <  1:2431280/150‑2 (MQ=255)
       cGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCATGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCa                                                                                                                                              >  2:1178189/1‑150 (MQ=255)
       cGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCa                                                                                                                                              >  2:2201570/1‑150 (MQ=255)
       cGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCa                                                                                                                                              <  1:1608086/150‑1 (MQ=255)
         aaCGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTCACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGa                                                                                                                                            <  1:30374/150‑1 (MQ=255)
             tACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGCCGCt                                                                                                                                        <  2:1523143/150‑1 (MQ=255)
               cccAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGa                                                                                                                                      <  2:974957/150‑1 (MQ=255)
                ccAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAg                                                                                                                                     >  2:1951412/1‑150 (MQ=255)
                   gCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCtt                                                                                                                                  >  2:2271273/1‑150 (MQ=255)
                      cTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGa                                                                                                                               <  1:1148430/150‑1 (MQ=255)
                       tctCAATTGCGCCTGCGGGGACATGCAGGCTTTACCGTGCAGGGCAGTAGCGCACTTGTCATCCCAGCTGGCGAAATGAGCGAACTGGCGGTGTCGGTGGCAATGCTCGCCGAACGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGaa                                                                                                                              <  2:300588/148‑1 (MQ=255)
                       caTCAATTGCGCCTGTCGGACGCTGACGGTTTCACGCTGCAGAGCAGTAGCGCACTGGTCACCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGaa                                                                                                                              <  2:186429/149‑1 (MQ=255)
                         tCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAAt                                                                                                                            >  1:2153938/1‑150 (MQ=255)
                           aaTTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCTCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCa                                                                                                                          <  2:413125/150‑1 (MQ=255)
                                  ccTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCagcgagcg                                                                                                                   >  1:1554277/1‑150 (MQ=255)
                                  ccTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCagcgagcg                                                                                                                   >  1:2250889/1‑150 (MQ=255)
                                   ctgctgGACGCTGACGGTTTCACCCTGCAGGGCAGTCGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGt                                                                                                                  >  1:660119/1‑150 (MQ=255)
                                   ctgctgGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGt                                                                                                                  <  2:1735280/150‑1 (MQ=255)
                                    tgctgGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCTGCTGCTCgggcgc                                                                                                                                           <  1:2085069/124‑6 (MQ=255)
                                    tgctgGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCTGCTGCTCgggcgc                                                                                                                                           >  2:2085069/1‑119 (MQ=255)
                                     gctgGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTga                                                                                                                <  1:1178189/150‑1 (MQ=255)
                                     gctgGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTCCGCGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTga                                                                                                                <  1:756572/150‑1 (MQ=255)
                                      gagtaCGCCGGCGGTTTCACATGGCAGGTCAGAAGCGAACGGGTCAGCCCATGTGGCGCAATGAGCGAATTGCCGGTGTCGGGGGAAATGCTCGCCGGCCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTGGAAATCAGCGAGCGTgac                                                                                                               <  2:925614/146‑1 (MQ=255)
                                        ggACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTgacga                                                                                                             <  2:1085385/150‑1 (MQ=255)
                                         gACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTgacgac                                                                                                            >  1:1367092/1‑150 (MQ=255)
                                            gCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTc                                                                                                                                                 >  1:592014/1‑110 (MQ=255)
                                            gCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTc                                                                                                                                                 <  2:592014/110‑1 (MQ=255)
                                             cTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCg                                                                                                        <  2:69728/150‑1 (MQ=255)
                                             cTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCg                                                                                                        >  1:2357338/1‑150 (MQ=255)
                                               gACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCag                                                                                                                         <  1:1587530/131‑1 (MQ=255)
                                               gACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCag                                                                                                                         >  2:1587530/1‑131 (MQ=255)
                                                 cGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAGCTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTAt                                                                                                    <  1:1208434/150‑1 (MQ=255)
                                                         cccTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCg                                                                                            >  1:1769194/1‑150 (MQ=255)
                                                         cccTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTGCGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCg                                                                                            >  1:1308682/1‑150 (MQ=255)
                                                          ccTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGt                                                                                           >  1:1780984/1‑150 (MQ=255)
                                                          ccTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGt                                                                                           >  2:1770195/1‑150 (MQ=255)
                                                           cTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTg                                                                                          >  1:29638/1‑150 (MQ=255)
                                                             gCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGgc                                                                                        <  2:29638/150‑1 (MQ=255)
                                                                gggCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGca                                                                                     <  2:1367092/150‑1 (MQ=255)
                                                                   cAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGt                                                                                                                  >  1:1524961/1‑118 (MQ=255)
                                                                   cAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGt                                                                                                                  <  2:1524961/118‑1 (MQ=255)
                                                                        gcgcACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCt                                                                             >  2:800485/1‑150 (MQ=255)
                                                                         aagccTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCATTGGCACTGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCtt                                                                            >  2:515629/5‑150 (MQ=255)
                                                                          gcACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTtc                                                                           <  1:2201570/150‑1 (MQ=255)
                                                                                     cccAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATcc                                                                                                  <  2:1766991/116‑1 (MQ=255)
                                                                                     cccAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATcc                                                                                                  >  1:1766991/1‑116 (MQ=255)
                                                                                        aGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGaa                                                             >  2:1742285/1‑150 (MQ=255)
                                                                                         gCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAAc                                                            <  2:2153938/150‑1 (MQ=255)
                                                                                             gCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTTGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCt                                                        <  1:1770195/150‑1 (MQ=255)
                                                                                             gCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCt                                                        <  2:2250889/150‑1 (MQ=255)
                                                                                                           cTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCCTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTc                                          >  2:480633/1‑150 (MQ=255)
                                                                                                                 gtgtCGGTGGCAATGCTCGCCGGGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCa                                    <  1:1918786/150‑1 (MQ=255)
                                                                                                                  tgtCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGCCGCTCAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCACCATGAACCGCTGATCCCTTACACTCACGCCAc                                   <  2:1246168/150‑1 (MQ=255)
                                                                                                                      ggTGGCAATGCTCGCCGAGAGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCACCCCCATGAACCGCTGATGCCTTACACTCACGCCACCTtg                               >  2:1488370/1‑150 (MQ=255)
                                                                                                                       gTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCCTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTtgc                              <  1:480633/150‑1 (MQ=255)
                                                                                                                         ggCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTtgcct                            <  1:2271273/150‑1 (MQ=255)
                                                                                                                           cAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGTGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCAc                                   <  2:1228867/141‑1 (MQ=255)
                                                                                                                           cAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGTGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCAc                                   >  1:1228867/1‑141 (MQ=255)
                                                                                                                            aaTGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTtgcctgcc                         >  1:537800/1‑150 (MQ=255)
                                                                                                                                cTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGc                     <  2:537800/150‑1 (MQ=255)
                                                                                                                                cTCGCCGAGCACCCGAGCATCAGCGCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGc                     >  2:635462/1‑150 (MQ=255)
                                                                                                                                  cGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCtt                   <  1:515629/150‑1 (MQ=255)
                                                                                                                                     cGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGa                <  2:159793/150‑1 (MQ=255)
                                                                                                                                        gtctgcGAGCTGCTGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGAACTCGGTCACTCACGCCACCTTGCCTGCCAGGCTTGGACAg             >  1:723990/6‑150 (MQ=255)
                                                                                                                                          gcCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGaa           <  1:1951412/150‑1 (MQ=255)
                                                                                                                                          gcCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGaa           <  2:823664/150‑1 (MQ=255)
                                                                                                                                            ccGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGAAgg         >  2:129507/1‑150 (MQ=255)
                                                                                                                                             cGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGaa           <  2:1787995/147‑1 (MQ=255)
                                                                                                                                             cGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGaa           >  1:1787995/1‑147 (MQ=255)
                                                                                                                                                   gcagcTCGCAGACGCTGAGCTTTGAAATCAGCGCCCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGAAGGCCAGAAt  >  1:2018672/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACAAGACCGAACGTACCCAGCACTATCAATTGCGCCTGCTGGACGCTGACGGTTTCACCCTGCAGGGCAGTAGCGCACTGGTCATCCCAGCTGGCGAAATGAGCGAACTGCCGGTGTCGGTGGCAATGCTCGCCGAGCGCCCGAGCAGCAGCTCGCAGACGCTGAGCTTTGAAATCAGCGAGCGTGACGACCCCGCTATCCGCAGCGTGGCGCACAGTCGCTTCGTCGCCCCCATGAACCGCTGATCCCTTACACTCACGCCACCTTGCCTGCCAGGCTTGGACAGAAGGCCAGAAT  >  NC_002947/4015848‑4016144

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