Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,002,454 C→T 5.0% intergenic (‑98/+48) amaA ← / ← amaB L‑pipecolate oxidase/L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,002,4540CT 5.0% 283.7 / 2.5 91intergenic (‑98/+48)amaA/amaBL‑pipecolate oxidase/L‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base T (3/2);  total (46/45)
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.71e-01

ACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGAC  >  NC_002947/6002307‑6002601
                                                                                                                                                   |                                                                                                                                                   
aCGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTa                                                                                                                                                   >  1:1472910/1‑150 (MQ=255)
   tgtTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTAccc                                                                                                                                                >  1:813293/1‑150 (MQ=255)
        ccACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGaa                                                                                                                                           >  2:436617/1‑150 (MQ=255)
         cacaAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAAc                                                                                                                                          >  1:569880/1‑150 (MQ=255)
          acaAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAAcc                                                                                                                                         >  2:1529592/1‑150 (MQ=255)
          acaAGCATTGTTGACGCAGTTCGGACATGGCCGGCTACAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAAcc                                                                                                                                         <  1:996547/150‑1 (MQ=255)
             aGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCggg                                                                                                                                      <  1:1529592/150‑1 (MQ=255)
              gCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCgggg                                                                                                                                     <  2:813293/150‑1 (MQ=255)
                       aCGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGAccc                                                                                                                            >  2:915503/1‑150 (MQ=255)
                             ttCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGcc                                                                                                                      >  1:993247/1‑150 (MQ=255)
                             gtCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTGGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGcc                                                                                                                      >  2:1123486/2‑150 (MQ=255)
                                ggACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATg                                                                                                                   >  1:833036/1‑150 (MQ=255)
                                 gACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCTGGCCCCCggttcgcg                                                                                                                  <  2:569880/150‑7 (MQ=255)
                                 gACATGGCCGGCTACATCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCggttcgcg                                                                                                                  >  1:470659/1‑144 (MQ=255)
                                  aCATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgca                                                                                                                 <  1:436617/150‑1 (MQ=255)
                                  aCATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgca                                                                                                                 <  1:1009078/150‑1 (MQ=255)
                                  aCAGGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCTAATTTCAGCGTGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCggttcgcgg                                                                                                                 <  2:470659/150‑8 (MQ=255)
                                     tGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgcagca                                                                                                              >  2:1301487/1‑150 (MQ=255)
                                         cGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCaa                                                                                                          <  1:1084651/150‑1 (MQ=255)
                                          ggCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAAt                                                                                                         <  1:1301487/150‑1 (MQ=255)
                                           gCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATa                                                                                                        >  2:1288528/1‑150 (MQ=255)
                                            cTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATAc                                                                                                       >  1:466314/1‑150 (MQ=255)
                                            cTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATAc                                                                                                       <  2:466314/150‑1 (MQ=255)
                                               cAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGaa                                                                                                                                           <  1:278691/111‑1 (MQ=255)
                                               cAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGaa                                                                                                                                           <  1:278698/111‑1 (MQ=255)
                                               cAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGaa                                                                                                                                           >  2:278698/1‑111 (MQ=255)
                                               cAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGaa                                                                                                                                           >  2:278691/1‑111 (MQ=255)
                                                aGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGGCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCTGCATCAATACATCa                                                                                                   >  1:337310/1‑150 (MQ=255)
                                                     tttGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGa                                                                                              >  2:849199/1‑150 (MQ=255)
                                                         ttCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGc                                                                                                                              <  2:894673/114‑1 (MQ=255)
                                                         ttCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGc                                                                                                                              >  1:894673/1‑114 (MQ=255)
                                                         ttCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgcagc                                                                                                               >  2:535422/1‑129 (MQ=255)
                                                         ttCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATgcagc                                                                                                               <  1:535422/129‑1 (MQ=255)
                                                            ttgttgGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGAt                                                                                       >  1:1199120/1‑150 (MQ=255)
                                                            ttgttgGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGAt                                                                                       <  1:1127358/150‑1 (MQ=255)
                                                            ttgttgGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGAt                                                                                       >  2:74282/1‑150 (MQ=255)
                                                               ttgGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATAcc                                                                                    >  2:1505721/1‑150 (MQ=255)
                                                                 ggggCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCt                                                                                  <  1:1459934/150‑1 (MQ=255)
                                                                   ggCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTgc                                                                                <  2:1472910/150‑1 (MQ=255)
                                                                   ggCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTgc                                                                                <  2:1525409/150‑1 (MQ=255)
                                                                    gCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAAt                                                                                                         >  1:1122054/1‑124 (MQ=255)
                                                                    gCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTgcg                                                                               >  1:475361/1‑150 (MQ=255)
                                                                    gCCGCAGGGCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAAt                                                                                                         <  2:1122054/124‑1 (MQ=255)
                                                                       gCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCa                                                                            >  1:753020/1‑150 (MQ=255)
                                                                        cAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAg                                                                           <  2:993247/150‑1 (MQ=255)
                                                                        cAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAg                                                                           >  1:281317/1‑150 (MQ=255)
                                                                               gcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCggggg                                                                                                                                    >  1:1496946/1‑86 (MQ=255)
                                                                               gcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCggggg                                                                                                                                    <  2:1496946/86‑1 (MQ=255)
                                                                               gcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCa                                                                                                                     >  2:317693/1‑101 (MQ=255)
                                                                               gcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCa                                                                                                                     <  1:317693/101‑1 (MQ=255)
                                                                               gcCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGc                                                                    <  2:753020/150‑1 (MQ=255)
                                                                                cccTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCt                                                                   >  1:381882/1‑150 (MQ=255)
                                                                                 ccTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCGGAGCTGCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTc                                                                  >  1:15929/1‑150 (MQ=255)
                                                                                   tttCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCaa                                                                                                          >  2:659781/1‑108 (MQ=255)
                                                                                   tttCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCaa                                                                                                          <  1:659781/108‑1 (MQ=255)
                                                                                        cgGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATTCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGag                                                           <  1:1123486/150‑1 (MQ=255)
                                                                                        cgGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGag                                                           <  1:1505721/150‑1 (MQ=255)
                                                                                        cgGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGag                                                           >  1:959338/1‑150 (MQ=255)
                                                                                         gggTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGagt                                                          <  1:460197/150‑1 (MQ=255)
                                                                                               ACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCa                                                    >  1:74282‑M1/1‑150 (MQ=255)
                                                                                               ACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCa                                                    <  2:1240548‑M1/150‑1 (MQ=255)
                                                                                                 CCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTcgcgc                                                              <  1:1614586‑M1/138‑1 (MQ=255)
                                                                                                 CCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTcgcgc                                                              >  2:1614586‑M1/1‑138 (MQ=255)
                                                                                                 CCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACg                                                  <  1:1553391‑M1/150‑1 (MQ=255)
                                                                                                 CCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACg                                                  <  1:1177275‑M1/150‑1 (MQ=255)
                                                                                                       ctacagggttcgtccaagcgacaccgacccctgtaggagcgggtTTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTg                                            >  2:520399‑M2/45‑150 (MQ=255)
                                                                                                       ctacagggttcatccaagcgacaccgacccctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATAcc                                                                                    <  1:1270003‑M2/66‑1 (MQ=255)
                                                                                                       ctacagggttcatccaagcgacaccgacccctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATAcc                                                                                    >  2:1270003‑M2/45‑110 (MQ=255)
                                                                                                        tacagggttcgtccaagcgacaccgacccctgtaggagcgggtTTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTgg                                           >  1:1274686‑M2/44‑150 (MQ=255)
                                                                                                        tacagggttcgtccaagcgacaccgacccctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCa                                                                      >  2:753335‑M2/44‑123 (MQ=255)
                                                                                                        tacagggttcgtccaagcgacaccgacccctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCa                                                                      <  1:753335‑M2/80‑1 (MQ=255)
                                                                                                         acagggttcgtccaagcgacaccgacccctgtaggagcgggtTTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCa                                                                                                   >  2:1589142‑M2/43‑93 (MQ=255)
                                                                                                         acagggttcgtccaagcgacaccgacccctgtaggagcgggtTTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCa                                                                                                   <  1:1589142‑M2/51‑1 (MQ=255)
                                                                                                             GGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCcgc                                                                                                                                              <  1:523738‑M1/46‑1 (MQ=255)
                                                                                                             GGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCcgc                                                                                                                                              >  2:523738‑M1/1‑46 (MQ=255)
                                                                                                                 cGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTTTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGAc                                  <  2:362212/150‑1 (MQ=255)
                                                                                                                      tCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATg                             >  2:448695/1‑150 (MQ=255)
                                                                                                                                 aaCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGACCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCa                  <  1:915503/150‑1 (MQ=255)
                                                                                                                                  aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCATCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAg                 <  2:475361/150‑1 (MQ=255)
                                                                                                                                  aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAg                 <  2:381882/150‑1 (MQ=255)
                                                                                                                                  aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAg                 <  1:1288528/150‑1 (MQ=255)
                                                                                                                                  aCCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAg                 <  2:1199120/150‑1 (MQ=255)
                                                                                                                                   ccctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGa                                                            >  2:716141‑M2/17‑106 (MQ=255)
                                                                                                                                   ccctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGa                                                            <  1:716141‑M2/90‑1 (MQ=255)
                                                                                                                                     ctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCa                                                    >  2:509806‑M2/15‑112 (MQ=255)
                                                                                                                                     ctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCa                                                    <  1:509806‑M2/98‑1 (MQ=255)
                                                                                                                                     ctgtaggagcgggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCa              >  2:1191394‑M2/15‑150 (MQ=255)
                                                                                                                                                ggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGaa                                                                                             >  2:627587‑M2/4‑60 (MQ=255)
                                                                                                                                                ggtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGaa                                                                                             <  1:627587‑M2/57‑1 (MQ=255)
                                                                                                                                                 gtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGAc  <  2:1553391‑M2/148‑1 (MQ=255)
                                                                                                                                                 gtCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGAc  >  2:1474913‑M2/3‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACGTGTTCCCACAAGCATTGTTGACGCAGTTCGGACATGGCCGGCTCCAGCAATTTGTTCTTGTTGGGGCCGCAGGGCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGTCGGTGTCGCTTGGACGAACCCTGTAGGAGCGGGTCTACCCGCGAACCGGGGGCCAGGCCCCCGGCCATGCAGCATCAATACATCAGTCGAACACGATACCCTGCGCCAGCGGCAGCTCGCGCGAGTAGTTCACGGTGTTGGTCTGGCGACGCATGTAGCCTTTCCAGGCATCCGAGCCCGAC  >  NC_002947/6002307‑6002601

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