Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 446,384 C→G 8.4% A357G (GCC→GGC)  PP_0368 → acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947446,3840CG 8.4% 219.4 / 8.3 71A357G (GCC→GGC) PP_0368acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (28/37);  new base G (3/3);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

TGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGT  >  NC_002947/446240‑446529
                                                                                                                                                 |                                                                                                                                                 
tGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGaa                                                                                                                                               <  8:193709/149‑1 (MQ=255)
tGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGaa                                                                                                                                               <  5:315331/149‑1 (MQ=255)
  gCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAgg                                                                                                                                             >  5:184933/1‑149 (MQ=255)
  gCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAgg                                                                                                                                             <  6:128006/149‑1 (MQ=255)
       ccAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAAcc                                                                                                                                        >  6:173078/1‑149 (MQ=255)
           ggTGGCCTGAATCAAGCCCGTGAGGGCCTGCTGATTCGTTCCCTGTCCGGCCCGACAGCACCGGA‑CAACGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGACGGCCTTCGGCGAAGGccttcat                                                                                                                                      <  7:203354/147‑7 (MQ=25)
             tGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGTCTTCGCCGAAGGCAACCGCGCGa                                                                                                                                  <  8:176387/149‑1 (MQ=255)
                 cTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGt                                                                                                                              >  5:209760/1‑149 (MQ=255)
                    aaGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGttt                                                                                                                            <  3:205608/148‑1 (MQ=255)
                    aaGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTa                                                                                                                           <  3:290890/149‑1 (MQ=255)
                      gTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACt                                                                                                                         <  4:56763/149‑1 (MQ=255)
                       tACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGGCGAAGGcccttca                                                                                                                                      <  8:235573/135‑7 (MQ=255)
                       tACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGGCGAAGGcccttca                                                                                                                                      >  7:235573/1‑129 (MQ=255)
                              gTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGcc                                                                                                                 <  4:156257/149‑1 (MQ=255)
                                  gcgcCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAg                                                                                                              >  1:186541/1‑148 (MQ=255)
                                  gcgcCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGc                                                                                                             <  3:124363/149‑1 (MQ=255)
                                         cAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTgg                                                                                                       <  6:184933/148‑1 (MQ=255)
                                         cAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTg                                                                                                        <  6:268438/147‑1 (MQ=255)
                                                gTTCCCTGACCGGCCCGAAAGCACCGGAGCAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGt                                                                                                >  1:19983/1‑149 (MQ=255)
                                                gTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGt                                                                                                >  5:162283/1‑148 (MQ=255)
                                                     cTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTa                                                                                          >  2:179120/1‑149 (MQ=255)
                                                      tGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTa                                                                                          >  5:125897/1‑148 (MQ=255)
                                                            gCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAg                                                                                   <  3:72241/149‑1 (MQ=255)
                                                            gCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAg                                                                                   <  1:179120/149‑1 (MQ=255)
                                                             cccGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAgg                                                                                  <  5:180735/149‑1 (MQ=255)
                                                              ccGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGa                                                                                 >  3:245366/1‑149 (MQ=255)
                                                                 aaaGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACga                                                                              <  2:19983/149‑1 (MQ=255)
                                                                 aaaGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACg                                                                               >  7:41879/1‑148 (MQ=255)
                                                                      aCCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACa                                                                                                    <  1:164923/122‑1 (MQ=255)
                                                                      aCCGGA‑CAAGGACGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACa                                                                                                    >  2:164923/1‑122 (MQ=255)
                                                                       ccGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTAc                                                                                         >  6:132686/1‑132 (MQ=255)
                                                                       ccGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTAc                                                                                         <  5:132686/132‑1 (MQ=255)
                                                                       ccGGA‑CAAGGCCGCTGACCCGATCATCGCCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaagaa                                                                          >  6:295440/1‑147 (MQ=255)
                                                                              aaGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                             >  7:120527/1‑122 (MQ=255)
                                                                              aaGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                             <  8:120527/122‑1 (MQ=255)
                                                                              aaGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAgcgc                                                                    >  6:307955/1‑147 (MQ=255)
                                                                                 gCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCaagaa                                                               <  5:173078/149‑1 (MQ=255)
                                                                                 gCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCACGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCaag                                                                 >  2:123096/1‑147 (MQ=255)
                                                                                   cGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTgg                                                                                                       <  4:39495/107‑1 (MQ=255)
                                                                                   cGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAgg                                                             >  8:262573/1‑149 (MQ=255)
                                                                                   cGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGACAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAgg                                                             <  8:231464/149‑1 (MQ=255)
                                                                                    gCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGc                                                                                      <  5:207520/123‑1 (MQ=255)
                                                                                    gCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGc                                                                                      >  6:207520/1‑123 (MQ=255)
                                                                                     cTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCt                                                           <  1:282816/149‑1 (MQ=255)
                                                                                        aCCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGa                                                         <  6:125897/148‑1 (MQ=255)
                                                                                          ccGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCaa                                                                                             <  6:101062/110‑1 (MQ=255)
                                                                                                ctcGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCctgc                                                 <  2:186728/147‑1 (MQ=255)
                                                                                                atcGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCctgct                                                <  6:33678/149‑1 (MQ=255)
                                                                                                       aCCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTc                                         >  5:188622/1‑149 (MQ=255)
                                                                                                       aCCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTc                                         <  3:123003/149‑1 (MQ=255)
                                                                                                        cccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTcc                                        >  8:272998/1‑149 (MQ=255)
                                                                                                         ccGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCt                                       <  4:245366/149‑1 (MQ=255)
                                                                                                              cGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGAccc                                  >  2:160913/1‑149 (MQ=255)
                                                                                                                  cGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGAt                              >  8:92895/1‑149 (MQ=255)
                                                                                                                    tCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCg                            <  1:160913/149‑1 (MQ=255)
                                                                                                                       cATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCa                         <  5:295440/149‑1 (MQ=255)
                                                                                                                         tgctgcTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAg                       <  6:188622/149‑1 (MQ=255)
                                                                                                                         tgctgcTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAg                       <  7:92895/149‑1 (MQ=255)
                                                                                                                             gctgACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCAt                   >  6:81739/1‑149 (MQ=255)
                                                                                                                                 aCCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGAGGCCCTGCTGGCGTTCCTGACCCAGATCGCCAAGGCATTCAt               >  5:279729/1‑149 (MQ=255)
                                                                                                                                 aCCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTc                 <  6:162283/147‑1 (MQ=255)
                                                                                                                                  ccATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATg              <  7:262573/149‑1 (MQ=255)
                                                                                                                                   cATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaa                                                                             >  2:191789/1‑85 (MQ=255)
                                                                                                                                   cATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACgaa                                                                             >  6:228673/1‑85 (MQ=255)
                                                                                                                                     tGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGAcc           <  2:186541/149‑1 (MQ=255)
                                                                                                                                      ggttgCCTTCGGCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGAcc           >  2:231047/5‑148 (MQ=255)
                                                                                                                                      ggttgCCTTCGGCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGAcc           >  1:296354/5‑148 (MQ=255)
                                                                                                                                      gAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGAcc           >  4:140452/1‑148 (MQ=255)
                                                                                                                                       gttgCCTTCGGCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGa         <  1:231047/146‑1 (MQ=255)
                                                                                                                                          gCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGt      <  7:181587/149‑1 (MQ=255)
                                                                                                                                             ttCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCgg   >  6:167430/1‑149 (MQ=255)
                                                                                                                                                 ccGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGt  <  2:324806/146‑1 (MQ=255)
                                                                                                                                                 ccGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGt  >  1:324806/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGGCGTTCCAGGGTGGCCTGAAGTACGCCCGTGAGCGCCTGCAGATGCGTTCCCTGACCGGCCCGAAAGCACCGGA‑CAAGGCCGCTGACCCGATCATCGTCCACCCGGACGTGCGTCGCATGCTGCTGACCATGAAGGCCTTCGCCGAAGGCAACCGCGCGATGGTTTACTTCACCGCCAAGCAGGTGGACATCGTCAAGTACAGCCAGGACGAAGAAGAGCGCAAGAAGGCTGACGCCCTGCTGGCGTTCCTGACCCCGATCGCCAAGGCATTCATGACCGAAGTCGGT  >  NC_002947/446240‑446529

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: ‑