Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,947 C→G 5.6% G375A (GGC→GCC) ‡ PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,9470CG 5.6% 399.9 / 11.0 124G375A (GGC→GCC) ‡PP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (56/61);  new base G (3/4);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 126 total aligned reads displayed.

TGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCGAGCTTC  >  NC_002947/4248798‑4249096
                                                                                                                                                     |                                                                                                                                                     
tGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCgcgc                                                                                                                                                       >  2:1803049/1‑150 (MQ=255)
    aggcaACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCcgc                                                                                                                                                   <  2:1694886/150‑1 (MQ=255)
           aGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGcgccg                                                                                                                                            >  2:1307516/1‑150 (MQ=255)
           aGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTCGCGTTCGCCGGACATCGCCGCCGCGTAGTCCAGCACCAGCATCGCCGCGCCCGCATcgcct                                                                                                                                            <  2:1529715/150‑2 (MQ=255)
            gTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCGGCGCCGCGGGCgcggcgatgc                                                                                                                                           <  1:742481/150‑5 (MQ=255)
             ttCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGcgcc                                                                                                                                             <  1:225349/147‑1 (MQ=255)
             ttCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGGAGTCCAGCACCAGCATCGCCGCGCCCGCGGcgcc                                                                                                                                             >  2:225349/1‑147 (MQ=255)
               ccGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCa                                                                                                                                        <  1:748102/150‑1 (MQ=255)
                  ggAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGCTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAggg                                                                                                                                     >  2:1434023/1‑150 (MQ=255)
                    aaaGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGgtg                                                                                                                                   <  1:586445/150‑1 (MQ=255)
                           aGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcgat                                                                                                                                             <  2:1323370/133‑3 (MQ=255)
                           aGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcgat                                                                                                                                             >  1:1323370/1‑131 (MQ=255)
                           aGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGa                                                                                                                            <  2:41196/150‑1 (MQ=255)
                               gACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCAt                                                                                                                        >  1:1210079/1‑150 (MQ=255)
                                  gAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCg                                                                                                                     >  2:513793/1‑150 (MQ=255)
                                    aGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGtt                                                                                                                   <  1:1029289/150‑1 (MQ=255)
                                     gCCGATCTTGCCGTGGGGTTCAACGGGATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTg                                                                                                                  >  1:263409/1‑150 (MQ=255)
                                     gCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCACCCGCTTTGCCGCTCAGGGTGTTGTCGAGCATGCGTTg                                                                                                                  >  2:645999/1‑150 (MQ=255)
                                              gCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCTGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCCGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGggcggg                                                                                                         <  2:605826/150‑2 (MQ=255)
                                              gCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAgg                                                                                                              >  1:1621743/1‑145 (MQ=255)
                                              gCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAgg                                                                                                              <  2:1621743/145‑1 (MQ=255)
                                              gCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGggcggc                                                                                                         <  1:1451912/150‑1 (MQ=255)
                                               ccGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCt                                                                                                        >  1:678820/1‑150 (MQ=255)
                                               ccGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCt                                                                                                        >  1:1618432/1‑150 (MQ=255)
                                                cGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGgcgc                                                                                                                                              <  1:1188035/111‑1 (MQ=255)
                                                cGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGgcgc                                                                                                                                              >  2:1188035/1‑111 (MQ=255)
                                                  tGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCa                                                                                                     >  1:366641/1‑150 (MQ=255)
                                                  tGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCa                                                                                                     >  1:1895974/1‑150 (MQ=255)
                                                      gTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGc                                                                                                 >  1:949484/1‑150 (MQ=255)
                                                      gTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGc                                                                                                 <  1:1292647/150‑1 (MQ=255)
                                                       ttCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGcc                                                                                                <  2:1210079/150‑1 (MQ=255)
                                                        tCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGccc                                                                                               <  2:1049749/150‑1 (MQ=255)
                                                            cGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGa                                                                                           <  2:231484/150‑1 (MQ=255)
                                                               cATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATgg                                                                                        <  1:1068176/150‑1 (MQ=255)
                                                                aTCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCa                                                                                                                         <  1:871505/116‑1 (MQ=255)
                                                                aTCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCa                                                                                                                         >  2:871505/1‑116 (MQ=255)
                                                                 tCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGgtg                                                                                                                                   >  2:1626799/1‑105 (MQ=255)
                                                                 tCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGgtg                                                                                                                                   <  1:1626799/105‑1 (MQ=255)
                                                                   gCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgt                                                                                    <  2:1491153/150‑1 (MQ=255)
                                                                    cAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtg                                                                                   <  2:1737958/150‑1 (MQ=255)
                                                                    cAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtg                                                                                   <  2:1618432/150‑1 (MQ=255)
                                                                       tACAGCCCCAGCAGCGGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCTCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGAGTTCCTGGGCGGCTTCATCGCCCCAGATGGCAgtggtg                                                                                <  1:1424965/150‑1 (MQ=255)
                                                                       tACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtggtg                                                                                >  2:941452/1‑150 (MQ=255)
                                                                         cAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGaa                                                                              <  1:1803049/150‑1 (MQ=255)
                                                                          aGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAg                                                                             >  2:1199039/1‑150 (MQ=255)
                                                                            ccccAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtggtg                                                                                <  2:615374/145‑1 (MQ=255)
                                                                            ccccAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAgtggtg                                                                                >  1:615374/1‑145 (MQ=255)
                                                                               cagcagcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATc                                                                                                   >  1:1671225/1‑123 (MQ=255)
                                                                               cagcagcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATc                                                                                                   <  2:1671225/123‑1 (MQ=255)
                                                                                agcagcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGt                                                                       <  2:678820/150‑1 (MQ=255)
                                                                                    gcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGt                                                                   <  1:1292273/150‑1 (MQ=255)
                                                                                    gcagGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGt                                                                   >  1:2633/1‑150 (MQ=255)
                                                                                        gTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTg                                                                      <  1:1390419/143‑1 (MQ=255)
                                                                                          cGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGt                                                                       >  2:1452050/1‑140 (MQ=255)
                                                                                          cGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGt                                                                       <  1:1452050/140‑1 (MQ=255)
                                                                                           gCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAg                                                            >  1:1119461/1‑150 (MQ=255)
                                                                                           gCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCACCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAg                                                            >  2:386301/1‑150 (MQ=255)
                                                                                            cATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGccc                                                                                               <  1:1443569/114‑1 (MQ=255)
                                                                                            cATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGccc                                                                                               >  2:1443569/1‑114 (MQ=255)
                                                                                             tcgccTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGt                                                          <  1:284416/147‑1 (MQ=255)
                                                                                             aTTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGt                                                          <  2:1287284/150‑1 (MQ=255)
                                                                                               tCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTg                                                        <  1:1144026/150‑1 (MQ=255)
                                                                                                 cTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGaa                                                                              >  1:1779935/1‑126 (MQ=255)
                                                                                                 cTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGc                                                      <  2:1895974/150‑1 (MQ=255)
                                                                                                  ttGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGcccc                                                                                              <  2:1289923/109‑1 (MQ=255)
                                                                                                  ttGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGcccc                                                                                              >  1:1289923/1‑109 (MQ=255)
                                                                                                  ttGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTAGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCa                                                     >  2:140836/1‑150 (MQ=255)
                                                                                                      gTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATgg                                                                                        >  2:1799651/1‑111 (MQ=255)
                                                                                                      gTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATgg                                                                                        <  1:1799651/111‑1 (MQ=255)
                                                                                                        tCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGccc                                               <  2:1119461/150‑1 (MQ=255)
                                                                                                          gCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCgg                                             >  1:1045117/1‑150 (MQ=255)
                                                                                                           ccGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcga                                                                                                                                              >  1:732171/1‑51 (MQ=255)
                                                                                                           ccGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGGGCgcggcga                                                                                                                                              <  2:732171/52‑2 (MQ=255)
                                                                                                            cGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGGGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCa                                           >  1:1860035/1‑150 (MQ=255)
                                                                                                             ggTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCa                                                                         >  1:1205574/1‑119 (MQ=255)
                                                                                                             ggTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCa                                                                         <  2:1205574/119‑1 (MQ=255)
                                                                                                                   cgccgcCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTc                                    <  2:1110468/150‑1 (MQ=255)
                                                                                                                   cgccgcCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGGTGTCGATCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGGCTAGGTTGGCATAGCCCCGCCACACCTc                                    >  2:330322/1‑150 (MQ=255)
                                                                                                                      cgcggGGTAGACCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCGCAGGGTTGTGTTGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTCGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGc                                 <  2:53198/150‑1 (MQ=255)
                                                                                                                           ggTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCATGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTg                            >  2:1594060/1‑150 (MQ=255)
                                                                                                                               gTCCAGCACCAGCATCGCCGCGCCCTCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTgccgc                        <  1:1420402/150‑1 (MQ=255)
                                                                                                                                tCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCg                                   >  2:1487511/1‑138 (MQ=255)
                                                                                                                                tCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCg                                   <  1:1487511/138‑1 (MQ=255)
                                                                                                                                tCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCACCAGATGGCAGTTTTGAATTCCAGTGGGGTGAACAGGTTGGCATAGCCCGGCAAAACCTCGACCAGTgccgcc                       >  1:1406008/1‑150 (MQ=255)
                                                                                                                                 ccagcaccagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCa                      >  1:879431/1‑150 (MQ=255)
                                                                                                                                  cagcaccagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCAc                     >  1:1897342/1‑150 (MQ=255)
                                                                                                                                  cagcaccagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCAc                     >  2:1127111/1‑150 (MQ=255)
                                                                                                                                      accagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCATGGCGGCTGCATAGCACCAGATGTCAGTGTTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGTCCAGTGCCGCCACCTGc                 <  2:263409/150‑1 (MQ=255)
                                                                                                                                      accagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTggg                                                                    <  2:665327/99‑1 (MQ=255)
                                                                                                                                      accagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTggg                                                                    >  1:665327/1‑99 (MQ=255)
                                                                                                                                        cagcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGccc               <  2:324613/150‑1 (MQ=255)
                                                                                                                                          gcaTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGc                                            >  2:204204/1‑119 (MQ=255)
                                                                                                                                           caTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAg                                                                                                               >  1:1931753/1‑51 (MQ=255)
                                                                                                                                           caTCGCCCCGCCCGCGGCGCGGCTCAGGGCGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGGTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTg            >  2:1445444/1‑150 (MQ=255)
                                                                                                                                            aTCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGc           <  1:513793/150‑1 (MQ=255)
                                                                                                                                             gcgccgcGGGCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTgg                                                       <  2:667679/104‑1 (MQ=255)
                                                                                                                                             gcgccgcGGGCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTgg                                                       >  1:667679/2‑105 (MQ=255)
                                                                                                                                              cgccgcGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTc         <  2:1831078/150‑1 (MQ=255)
                                                                                                                                              cgccgcGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCAGACAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTc         <  1:941452/150‑1 (MQ=255)
                                                                                                                                                     cccGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCTTCCAGTGCCGCCACCTGCCCTTGCTCGAGCTTc  <  1:645999/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCAGGCAACAGTTCCGGGAAAGCCGAGGTGACGAAGCCGATCTTGCCGTGGCGTTCAACGGCATCGCAGTACAGCCCCAGCAGCAGGTCGCATTCCTTGCGTTCGCCGGTAAACGCCGCCGGGTAGTCCAGCACCAGCATCGCCGCGCCCGCGGCGCCGCTCAGGGTGTTGTCGAGCATGCGTTGCAGGGCGGCTTCATCGCCCCAGATGGCAGTGGTGAAGTCCAGTGGGTTGACCAGGTTGGCATAGCCCGGCAACACCTCGGCCAGTGCCGCCACCTGCCCTTGCTCGAGCTTC  >  NC_002947/4248798‑4249096

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