Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,688,230 A→C 5.8% K102Q (AAG→CAG)  PP_3253 → carboxylate‑amine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,688,2300AC 5.8% 198.3 / 4.2 69K102Q (AAG→CAG) PP_3253carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base A (39/26);  new base C (2/2);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

CGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTT  >  NC_002947/3688082‑3688378
                                                                                                                                                     |                                                                                                                                                     
cGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGaa                                                                                                                                                      <  1:1107967/150‑1 (MQ=255)
  aGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGTCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGc                                                                                                                                                    >  2:981157/1‑150 (MQ=255)
  aGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGc                                                                                                                                                    >  2:1284209/1‑150 (MQ=255)
      ggCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCTTGTGCCGCG‑GGCTTCTGCAGCAGccacgc                                                                                                                                               <  2:1374182/150‑5 (MQ=255)
       gCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑cgcg                                                                                                                                              <  2:607369/150‑1 (MQ=255)
              cGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAg                                                                                                                                       <  2:1242015/150‑1 (MQ=255)
                gtgtTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGGTTGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAggg                                                                                                                                     <  1:604547/150‑1 (MQ=255)
                     cACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACta                                                                                                                                >  1:1607459/1‑150 (MQ=255)
                        cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGTCGAAGAGGGGATGGGGCCATACTCTGCGGCAA‑GCCACCCCTGCGCAGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATaa                                                                                                                             >  1:1754135/1‑150 (MQ=255)
                        cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                               <  2:818224/148‑1 (MQ=255)
                        cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACtat                                                                                                                               >  1:818224/1‑148 (MQ=255)
                        cAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATaa                                                                                                                             >  1:1013538/1‑150 (MQ=255)
                          aCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAag                                                                                                                            <  2:1795360/149‑1 (MQ=255)
                          aCCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAag                                                                                                                            >  1:1795360/1‑149 (MQ=255)
                            cTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         <  2:1134571/150‑1 (MQ=255)
                             tGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCGTTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcagc                                                                                                                        >  1:48504/1‑150 (MQ=255)
                               tATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                             <  1:103743/127‑1 (MQ=255)
                               tATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCG‑GGCTTCTGCAGCAGCCACGCgg                                                                                                                                             >  2:103743/1‑127 (MQ=255)
                                aTGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         >  2:918223/1‑146 (MQ=255)
                                aTGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         <  1:918223/146‑1 (MQ=255)
                                   aaGCACGTGAGTTTCTCCAGCGCAGCCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcg                                                                                                                  >  1:599578/1‑150 (MQ=255)
                                    aGCACGTGAGTTTCTCCAGCGCAACAGTAAGAGGCTGAGAGTGTCATTGGTCGAAGAGGGGATGGGGCCATACTCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                                 <  2:1754135/150‑1 (MQ=255)
                                         gTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTAGACGACt                                                                                                            >  1:392814/1‑150 (MQ=255)
                                          tGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                                 >  1:1769645/1‑144 (MQ=255)
                                          tGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTcga                                                                                                                 <  2:1769645/144‑1 (MQ=255)
                                               ttCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGcc                                                                                                      >  2:1756039/1‑150 (MQ=255)
                                               ttCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGcc                                                                                                      >  2:224497/1‑150 (MQ=255)
                                                 ctcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCAt                                                                                                    <  1:1658318/150‑1 (MQ=255)
                                                  tcCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCAtg                                                                                                   >  1:225009/1‑150 (MQ=255)
                                                    cAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACGCCGTCTCTATAAGCAGCTGTTCGACGACTACCGCCAtgtg                                                                                                 >  1:1039649/1‑150 (MQ=255)
                                                      gcgcAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGc                                                                                               >  2:1371756/1‑150 (MQ=255)
                                                        gcAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAAGGTTACGGCTG‑GCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCaa                                                                                             <  1:797812/150‑1 (MQ=255)
                                                                  gCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         >  1:1629104/1‑112 (MQ=255)
                                                                  gCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         <  2:1629104/112‑1 (MQ=255)
                                                                     gCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCt                                                                                <  1:1011127/150‑1 (MQ=255)
                                                                       tGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCa                                                                              >  1:735168/1‑150 (MQ=255)
                                                                            gTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGc                                                                         >  2:378501/1‑150 (MQ=255)
                                                                               gCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTg                                                                      >  2:300294/1‑150 (MQ=255)
                                                                                  ttGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCAt                                                                   <  2:735168/150‑1 (MQ=255)
                                                                                        gAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                                            <  2:1544947/119‑1 (MQ=255)
                                                                                        gAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAg                                                                                            >  1:1544947/1‑119 (MQ=255)
                                                                                                   ggggCCATACGCTGCGGCAA‑GCCACCGCGGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgcc                                                  <  2:439797/150‑1 (MQ=255)
                                                                                                   ggggCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTgcc                                                  >  2:176104/1‑150 (MQ=255)
                                                                                                     ggCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                                     >  1:785029/1‑113 (MQ=255)
                                                                                                     ggCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGc                                                                                     <  2:785029/113‑1 (MQ=255)
                                                                                                          tACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGTACTATAAGCAGCCGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCt                                           >  2:373455/1‑150 (MQ=255)
                                                                                                                   gCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         >  2:605525/1‑63 (MQ=255)
                                                                                                                   gCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagcag                                                                                                                         <  1:605525/63‑1 (MQ=255)
                                                                                                                         ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACCAGGGCATTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCa                             >  1:1883918/1‑150 (MQ=255)
                                                                                                                         ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                                          <  1:363477/57‑1 (MQ=255)
                                                                                                                         ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAagca                                                                                                                          >  2:363477/1‑57 (MQ=255)
                                                                                                                         ccACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCa                             >  1:502776/1‑150 (MQ=255)
                                                                                                                          cACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                            >  2:1863199/1‑150 (MQ=255)
                                                                                                                          cACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCaa                            <  2:48504/150‑1 (MQ=255)
                                                                                                                           aCCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAAc                           >  1:259837/1‑150 (MQ=255)
                                                                                                                             ccttgtgccgcG‑GGCTTCTGCAGCAGCCACGCGGCGCAGGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTg                         >  2:1673461/7‑150 (MQ=37)
                                                                                                                              ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                                      <  1:121763/136‑1 (MQ=255)
                                                                                                                              ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGa                                      >  2:121763/1‑136 (MQ=255)
                                                                                                                              ccTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGa                        <  1:247529/150‑1 (MQ=255)
                                                                                                                                   gccgcGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAAc                   >  2:686074/1‑150 (MQ=255)
                                                                                                                                      gcgTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGc                <  1:300294/150‑1 (MQ=255)
                                                                                                                                            ctgctgCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTg          >  2:575005/1‑150 (MQ=255)
                                                                                                                                              gctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGcc        <  2:1619407/150‑1 (MQ=255)
                                                                                                                                               ctgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGccc       >  2:184093/1‑150 (MQ=255)
                                                                                                                                                tgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGc                                            <  2:311157/112‑1 (MQ=255)
                                                                                                                                                tgcAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGc                                            >  1:311157/1‑112 (MQ=255)
                                                                                                                                                  cAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTgg    >  1:1010944/1‑150 (MQ=255)
                                                                                                                                                   aGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGt                                                     >  2:131741/1‑100 (MQ=255)
                                                                                                                                                   aGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGt                                                     <  1:131741/100‑1 (MQ=255)
                                                                                                                                                    gAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGtt  >  1:1515938/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAGCTGGCGTCGCCGGTGTTCACCAACCTGTATGAAGCACGTGAGTTTCTCCAGCGCAACCGTCAGCGGCTGCGCGTGGCATTGGCCGAAGAGGGGATGGGGCCATACGCTGCGGCAA‑GCCACCCCTGCGCCGCGTGGCTGCTGCAGAAGCC‑CGCGGCACAAGGGCACTATAAGCAGCTGTTCGACGACTACCGCCATGTGGCAAGGCGTAGCCTGCTCAATGGCCTGCATGTGCACGTTGGCGTGCCGCCAGCTTGTGACCGCATGCAACTGATCAACCGCCTTCTGCCCTGGTT  >  NC_002947/3688082‑3688378

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