Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,250,427 G→T 6.3% F274L (TTC→TTA PP_3725 ← acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,250,4270GT 6.3% 196.7 / 2.3 63F274L (TTC→TTAPP_3725acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base T (2/2);  total (31/32)
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.76e-01

CTTGGCCATGCCGGCCTCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCACCTTGC  >  NC_002947/4250280‑4250575
                                                                                                                                                   |                                                                                                                                                    
cTTGGCCATGCCGGCCTCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCCGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGaa                                                                                                                                                    <  1:1073940/150‑1 (MQ=255)
    gCCATGCCGGCCTCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGc                                                                                                                                                <  2:1123751/150‑1 (MQ=255)
       aTGTCGGCCTCACTGTCGGTCATTTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTg                                                                                                                                             <  2:913096/150‑1 (MQ=255)
       aTGCCGGCCTCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTg                                                                                                                                             <  2:91839/150‑1 (MQ=255)
             gCCTCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCGGCAGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTcg                                                                                                                                       >  1:655320/1‑150 (MQ=255)
                tCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTcggcg                                                                                                                                    >  1:999869/1‑150 (MQ=255)
                  aCCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTTGCGGTGGTCTGTCCGAACTGCTTGCGGTCggcggc                                                                                                                                  >  2:67413/1‑150 (MQ=255)
                      tCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATAGCGGCCAGCTTGAACGACACCCCNTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCAt                                                                                                                              <  1:819248/150‑1 (MQ=255)
                        ggTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTg                                                                                                                            >  1:322013/1‑150 (MQ=255)
                          tCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCa                                                                                                                          >  1:282048/1‑150 (MQ=255)
                           cATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAg                                                                                                                         <  2:167476/150‑1 (MQ=255)
                              gTTGCCCTGGTCCATCTTCCAGGCGGAGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCGTGATAGCTGCCGATGGGCTGGCCGAAATGCTTGCGGTCGGCGGCCCACTGAATCGa                                                                                                                      >  2:354513/1‑150 (MQ=255)
                                       gTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCTCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCTAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAg                                                                                                             >  1:590690/1‑150 (MQ=255)
                                         ccATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCa                                                                                                                          >  2:784183/1‑135 (MQ=255)
                                         ccATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCa                                                                                                                          <  1:784183/135‑1 (MQ=255)
                                              ttCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGgcgcgc                                                                                                      >  1:6732/1‑150 (MQ=255)
                                                 cAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTgg                                                                                                   <  1:354513/150‑1 (MQ=255)
                                                  aGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACCCCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTggg                                                                                                  >  2:243061/1‑150 (MQ=255)
                                                     cGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCt                                                                                               <  2:6732/150‑1 (MQ=255)
                                                       gtgtGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTgg                                                                                             <  1:959384/150‑1 (MQ=255)
                                                         gtgCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGcc                                                                                           >  1:1215160/1‑150 (MQ=255)
                                                               ggTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCa                                                                                     <  1:201213/150‑1 (MQ=255)
                                                                 tGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCTAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGt                                                                                   >  2:1098351/1‑150 (MQ=255)
                                                                     cATCTCGGCCGCTCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCTAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTggc                                                                               <  2:590690/150‑1 (MQ=255)
                                                                     cATCTCGGCCGCGCGGATCTGCGTGGCCATAACGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTggc                                                                               >  1:47324/1‑150 (MQ=255)
                                                                     cATCGCGGCCGCGCGGATATGCGTGGCACTCTCGGCCAGTCGGTACGCATTCGCCGTATAGCGGCCGAGGGCCTGTCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTggc                                                                               <  2:655320/150‑1 (MQ=255)
                                                                          cGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGa                                                                          <  2:673302/150‑1 (MQ=255)
                                                                            gCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGAcc                                                                        <  1:5132/150‑1 (MQ=255)
                                                                            gCCACGCAGATCCGCGCGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGAcc                                                                        >  2:856012/1‑150 (MQ=255)
                                                                               gcgcgGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCAACGCAGTTGGCGGCGACCATc                                                                     >  1:488034/1‑150 (MQ=255)
                                                                                 gcgGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCAc                                                                   <  2:320708/150‑1 (MQ=255)
                                                                                  cgGATCTGCGTGGCCATATCGGCGAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCAcc                                                                  <  2:414776/150‑1 (MQ=255)
                                                                                           gTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCAACGCAGTTGGCGGCGACCATCACCCGGCcggcg                                                         <  2:488034/150‑1 (MQ=255)
                                                                                              gCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTa                                                      >  1:828405/1‑150 (MQ=255)
                                                                                                 atatCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTACGc                                                   <  1:1283637/150‑1 (MQ=255)
                                                                                                   atCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCa                                                 >  1:375137/1‑150 (MQ=255)
                                                                                                    ccgatctgCTTGAACGACACCCCCTGATAGCTGCCGATACAGGGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAg                                                <  1:654989/143‑1 (MQ=255)
                                                                                                         cAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCTAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGtt                                           <  1:1098351/150‑1 (MQ=255)
                                                                                                           gCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTgg                                         >  2:242136/1‑150 (MQ=255)
                                                                                                           gCTGGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTAGGCGGCCCATTGCAGCGACAGGGCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTgg                                         <  2:362430/150‑1 (MQ=255)
                                                                                                            cTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGc                                        >  1:689952/1‑150 (MQ=255)
                                                                                                               gAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCAc                                     >  1:482423/1‑150 (MQ=255)
                                                                                                                 aCGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTgg                                                                                             <  1:929590/92‑1 (MQ=255)
                                                                                                                 aCGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTgg                                                                                             >  2:929590/1‑92 (MQ=255)
                                                                                                                   gACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTcc                                 >  2:252748/1‑150 (MQ=255)
                                                                                                                     cacCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGc                                                   >  1:206105/1‑130 (MQ=255)
                                                                                                                     cacCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGc                                                   <  2:206105/130‑1 (MQ=255)
                                                                                                                      acCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATTGAACACCTGTCGGGAAGCCAGGCGTTGGCCACTTCCCAg                              >  2:744224/1‑150 (MQ=255)
                                                                                                                        ccccTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGcc                            >  1:57254/1‑150 (MQ=255)
                                                                                                                          ccTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCt                          >  2:1140909/1‑150 (MQ=255)
                                                                                                                           cTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCtt                         <  2:57254/150‑1 (MQ=255)
                                                                                                                           cTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCtt                         >  1:960329/1‑150 (MQ=255)
                                                                                                                             gATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGc                       <  1:1181440/150‑1 (MQ=255)
                                                                                                                              aTAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGcc                      <  2:47324/150‑1 (MQ=255)
                                                                                                                                          ggcctggccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAg          <  1:242136/150‑1 (MQ=255)
                                                                                                                                            cctggccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCa        >  2:104113/1‑150 (MQ=255)
                                                                                                                                            cctggccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCa        >  2:917404/1‑150 (MQ=255)
                                                                                                                                                gccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGc                       <  2:551287/131‑1 (MQ=255)
                                                                                                                                                gccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGc                       >  1:551287/1‑131 (MQ=255)
                                                                                                                                                 ccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCACCTTg   <  2:960329/150‑1 (MQ=255)
                                                                                                                                                 ccGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCACCTTg   <  2:375137/150‑1 (MQ=255)
                                                                                                                                                  cGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCACCTTGc  <  1:917404/150‑1 (MQ=255)
                                                                                                                                                  cGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCACCTTGc  >  1:1048516/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTGGCCATGCCGGCCTCACCGTCGGTCATGTTGCCCTGGTCCATCTTCCAGGCGGTGTGCAGGGTGAGCATCTCGGCCGCGCGGATCTGCGTGGCCATATCGGCCAGCTTGAACGACACCCCCTGATAGCTGCCGATGGCCTGGCCGAACTGCTTGCGGTCGGCGGCCCATTGCAGCGACAGGTCGAGGGCGCGCTGGGCCTGGCCCACGCAGTTGGCGGCGACCATCACCCGGCCGGCGGTAAGCCAGGCGTTGGCCACTTCCCAGCCCTTGCCCACTTCGCCCAGCACCTTGC  >  NC_002947/4250280‑4250575

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