Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,330,246 A→G 5.9% Q200R (CAG→CGG)  PP_2048 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,330,2460AG 5.9% 204.7 / 2.3 68Q200R (CAG→CGG) PP_2048acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base G (2/2);  total (34/34)
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.54e-01

GTCAGTGGCGAGTGGGCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGAGCACAA  >  NC_002947/2330101‑2330393
                                                                                                                                                 |                                                                                                                                                   
gTCAGTGGCGAGTGGGCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAgg                                                                                                                                                   <  2:52013/148‑1 (MQ=255)
gTCAGTGGCGAGTGGGCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGa                                                                                                                                                  <  4:211665/149‑1 (MQ=255)
 tCAGTGGCGAGTGGGCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGAc                                                                                                                                                 <  8:144388/149‑1 (MQ=255)
       gCGAGTGGGCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGAcc                                                                                                                                           >  8:130414/1‑149 (MQ=255)
               gCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACAt                                                                                                                                   <  7:130414/149‑1 (MQ=255)
               gCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCAGTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCCGGACCTGACCGACACCAt                                                                                                                                   >  8:131/1‑149 (MQ=255)
                    cACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTc                                                                                                                              <  2:280595/149‑1 (MQ=255)
                     aCCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCa                                                                                                                             >  8:274048/1‑149 (MQ=255)
                        aTGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCt                                                                                                                          <  6:369728/149‑1 (MQ=255)
                                   cGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCGGGACCTGACCGACAACATTGTTCATCTGGTACTcgcgc                                                                                                               >  6:35775/1‑149 (MQ=255)
                                     aGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCgcgcgc                                                                                                             >  8:103324/1‑149 (MQ=255)
                                      gCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTTTTCATCTGGTCATcgcgcgag                                                                                                            >  6:265057/1‑147 (MQ=255)
                                      gCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACAGCATTACCGTCAGCAATATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTcgcgc                                                                                                               <  5:265057/146‑1 (MQ=255)
                                       ccGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCgcgcgc                                                                                                             <  2:128799/147‑1 (MQ=255)
                                       ccGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCgcgcgc                                                                                                             >  1:128799/1‑147 (MQ=255)
                                             gCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGa                                                                                                     >  2:26102/1‑149 (MQ=255)
                                             gCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGa                                                                                                     >  3:61419/1‑149 (MQ=255)
                                             gCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGa                                                                                                     >  6:346575/1‑149 (MQ=255)
                                             gCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGa                                                                                                     >  8:243524/1‑149 (MQ=255)
                                              ccGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCgcctgcc                                                                                                        <  3:275858/145‑1 (MQ=255)
                                              ccGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCgcctgcc                                                                                                        >  4:275858/1‑145 (MQ=255)
                                                 gCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATgcg                                                                                                 >  5:262959/1‑149 (MQ=255)
                                                   aCCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCgcc                                                                                               >  1:328770/1‑149 (MQ=255)
                                                     cGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCgccgg                                                                                             >  5:36654/1‑149 (MQ=255)
                                                       aCCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTcgc                                                                                                                 >  6:181903/1‑127 (MQ=255)
                                                       aCCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTcgc                                                                                                                 <  5:181903/127‑1 (MQ=255)
                                                               ttGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAg                                                                                   <  6:36654/149‑1 (MQ=255)
                                                                 gctgcGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAgg                                                                                  <  3:332725/148‑1 (MQ=255)
                                                                  ctgcGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCGGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAggg                                                                                 <  5:35775/148‑1 (MQ=255)
                                                                         cccGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGc                                                                         <  6:262959/149‑1 (MQ=255)
                                                                          ccGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGcc                                                                        >  4:123319/1‑149 (MQ=255)
                                                                                  aGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTg                                                                <  3:123319/149‑1 (MQ=255)
                                                                                          gCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGccc                                                          >  1:322197/1‑147 (MQ=255)
                                                                                              aCGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAgctg                                                    >  5:217612/1‑149 (MQ=255)
                                                                                                  cAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATATGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCtgctgc                                                >  4:230964/1‑149 (MQ=255)
                                                                                                      tACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCgg                                            >  6:252033/1‑149 (MQ=255)
                                                                                                        ccGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTg                                          <  5:252033/149‑1 (MQ=255)
                                                                                                        ccGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTg                                          <  3:230964/149‑1 (MQ=255)
                                                                                                         cGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAgg                                                                                  >  1:139806/1‑108 (MQ=255)
                                                                                                         cGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAgg                                                                                  <  2:139806/108‑1 (MQ=255)
                                                                                                         cGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTgg                                         <  2:328770/149‑1 (MQ=255)
                                                                                                             ttACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAg                                     >  5:154378/1‑149 (MQ=255)
                                                                                                             ttACCGGAAGCAATTTATTCTTGAGTGTTGGCGAGCAGGACCTTACCGACAAAATTGTTCATCTGGTAGTCGCGCGCCTGCCTGATGCGCCGGACGGTGTCAGGGGCATCAGCCTGTTCATGGTTCCCAAGCTGCTGCCCGGTGGCGAg                                     <  6:92966/149‑1 (MQ=255)
                                                                                                                       caAGATATTCATCAGTGGTGGCGAGCGGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGTCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATgcgc                           >  7:321135/1‑149 (MQ=255)
                                                                                                                        aaGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCt                          >  1:82562/1‑149 (MQ=255)
                                                                                                                         aGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTg                         <  1:26102/149‑1 (MQ=255)
                                                                                                                         aGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTg                         <  7:274048/149‑1 (MQ=255)
                                                                                                                           atatTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCAt                                                                             >  5:77815/1‑95 (MQ=255)
                                                                                                                           atatTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCAt                                                                             <  6:77815/95‑1 (MQ=255)
                                                                                                                           atatTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTggg                       >  5:371845/1‑149 (MQ=255)
                                                                                                                           atatTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTggg                       >  7:29520/1‑149 (MQ=255)
                                                                                                                             atTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCt                     >  6:387052/1‑149 (MQ=255)
                                                                                                                                 atcaGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGc                                                                         >  4:204472/1‑93 (MQ=255)
                                                                                                                                 atcaGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGc                                                                         <  3:204472/93‑1 (MQ=255)
                                                                                                                                 atcaGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCg                  <  4:295523/148‑1 (MQ=255)
                                                                                                                                    aGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCgggg               <  7:103324/148‑1 (MQ=255)
                                                                                                                                      tggtggCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGgcgc            >  5:104279/1‑149 (MQ=255)
                                                                                                                                        gtggCGAGCGGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGTCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTg          <  8:321135/149‑1 (MQ=255)
                                                                                                                                        gtggCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTg          <  6:217612/149‑1 (MQ=255)
                                                                                                                                        gtggCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTg          >  6:130461/1‑149 (MQ=255)
                                                                                                                                          ggCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGa        <  5:346575/149‑1 (MQ=255)
                                                                                                                                           gCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCTTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGa        <  4:61419/148‑1 (MQ=255)
                                                                                                                                           gCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGa        <  5:130461/148‑1 (MQ=255)
                                                                                                                                            cGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGAGc      <  8:29520/149‑1 (MQ=255)
                                                                                                                                             gAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGAGca     >  4:299296/1‑149 (MQ=255)
                                                                                                                                                cAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCgg                 >  7:359112/1‑134 (MQ=255)
                                                                                                                                                cAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCgg                 <  8:359112/134‑1 (MQ=255)
                                                                                                                                                cAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGAGCACaa  <  5:387052/149‑1 (MQ=255)
                                                                                                                                                 aGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGAGcac    >  1:194729/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTCAGTGGCGAGTGGGCCGGCACCATGTGCCTGACCGAGCCGCAGGCCGGCACCGACCTGGCGTTGCTGCGTACCCGCGCCCAGCCGGAAGCGGACGGCAGCTACCGCATTACCGGCAGCAAGATATTCATCAGTGGTGGCGAGCAGGACCTGACCGACAACATTGTTCATCTGGTACTCGCGCGCCTGCCTGATGCGCCGGCCGGTGTCAGGGGCATCAGCCTGTTCCTGGTGCCCAAGCTGCTGCCCGGTGGCGAGCGCAATGCGCTGGGCTGCGGGGCGCTGGAGCACAA  >  NC_002947/2330101‑2330393

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