Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,134,025 T→A 8.6% Q350L (CAG→CTG)  PP_3638 ← acyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,134,0250TA 8.6% 209.6 / 10.3 70Q350L (CAG→CTG) PP_3638acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base A (4/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

CACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGT  >  NC_002947/4133877‑4134172
                                                                                                                                                    |                                                                                                                                                   
caccacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGcctgcctg                                                                                                                                                    <  1:1647111/150‑1 (MQ=255)
 accacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTgg                                                                                                                                                   <  2:12083/150‑1 (MQ=255)
   cacTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCg                                                                                                                                                 >  2:1597873/1‑150 (MQ=255)
        gggCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCa                                                                                                                                            <  1:737447/150‑1 (MQ=255)
          gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCACCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                                          >  2:291209/1‑150 (MQ=255)
          gCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGt                                                                                                                                          >  2:723647/1‑150 (MQ=255)
              aGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGACGTCACGGGCCTCCAGCAACCCCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCgg                                                                                                                                      <  1:539009/150‑1 (MQ=255)
                       aCGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAAtt                                                                                                                             >  2:1054013/1‑150 (MQ=255)
                          aCCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCa                                                                                                                          <  2:895955/150‑1 (MQ=255)
                           tccgatcTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCTCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCAc                                                                                                                         <  2:1041450/144‑1 (MQ=255)
                             tCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGTCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCa                                                                                                                       >  1:57633/1‑150 (MQ=255)
                               aCCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGTCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                                     <  2:57633/150‑1 (MQ=255)
                               aCCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgg                                                                                                                     >  1:649672/1‑150 (MQ=255)
                                 cTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCt                                                                                                                   >  2:196293/1‑150 (MQ=255)
                                  tGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAGCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTg                                                                                                                  >  1:279951/1‑150 (MQ=255)
                                    gTGCCAGCTGCAGGTCGAACTGCACAAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAgatcgga                                                                                                                >  1:1041450/1‑144 (MQ=255)
                                     tGCCAGCTGCAGGTCGCACTGCCCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                               <  1:1426634/150‑1 (MQ=255)
                                         aGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                           >  1:1373089/1‑150 (MQ=255)
                                         aGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                           <  1:1343463/150‑1 (MQ=255)
                                         aGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGa                                                                                                           <  2:1286053/150‑1 (MQ=255)
                                             gCAGGTCGAACCGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCATTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGAt                                                                                                       <  1:1523232/150‑1 (MQ=255)
                                                 gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcctgcct                                                                                                                                                     >  2:503250/1‑100 (MQ=255)
                                                 gTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACgcctgcct                                                                                                                                                     <  1:503250/100‑1 (MQ=255)
                                                      aCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                   <  1:1583373/145‑1 (MQ=255)
                                                      aCTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGc                                                                                                   >  2:1583373/1‑145 (MQ=255)
                                                       cTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAg                                                                                             <  1:311402/150‑1 (MQ=255)
                                                       cTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTCGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAg                                                                                             <  2:885301/150‑1 (MQ=255)
                                                          cgcCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGcc                                                                                          <  2:1373089/150‑1 (MQ=255)
                                                            ccAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGCAGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGa                                                                                        >  2:714199/1‑150 (MQ=255)
                                                             cAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCAGGCAGGCGTTGCa                                                                                                                                            <  1:1017582/97‑1 (MQ=255)
                                                             cAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCAGGCAGGCGTTGCa                                                                                                                                            >  2:1017582/1‑97 (MQ=255)
                                                               gCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTc                                                                                     >  2:1362087/1‑150 (MQ=255)
                                                                         ccAGCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcgg                                                                           <  2:279951/150‑1 (MQ=255)
                                                                           gaccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                         >  2:1096489/2‑150 (MQ=255)
                                                                           aaCCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcg                                                                         >  2:646457/1‑150 (MQ=255)
                                                                             ccaccCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCggg                                                                       >  2:869171/1‑150 (MQ=255)
                                                                                 ccTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCaaa                                                                   <  2:649672/150‑1 (MQ=255)
                                                                                     aCGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCt                                                               >  2:1492957/1‑150 (MQ=255)
                                                                                      cgcgGGTAGCCCAAGCTGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCtt                                                              <  2:772763/150‑1 (MQ=255)
                                                                                                ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGaa                                                                                              >  2:1430359/1‑108 (MQ=255)
                                                                                                ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGaa                                                                                              <  1:1430359/108‑1 (MQ=255)
                                                                                                ccAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTc                                                    >  1:1645911/1‑150 (MQ=255)
                                                                                                  aaGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGc                                                  >  2:637420/1‑150 (MQ=255)
                                                                                                    gCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                               <  2:180348/87‑1 (MQ=255)
                                                                                                    gCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGc                                                                                                               >  1:180348/1‑87 (MQ=255)
                                                                                                          gTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGACTGCCTGGCGCCGTAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTaa                                          >  1:1265851/1‑150 (MQ=255)
                                                                                                             aCGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCAGGGCGGGACGAGCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGcc                                       >  1:197252/1‑150 (MQ=255)
                                                                                                                  ccTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTgcc                                  <  1:1492957/150‑1 (MQ=255)
                                                                                                                     ccAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                                   >  1:1465963/1‑130 (MQ=255)
                                                                                                                     ccAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCg                                                   <  2:1465963/130‑1 (MQ=255)
                                                                                                                      cAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCa                              >  1:1639881/1‑150 (MQ=255)
                                                                                                                       aGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTTGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCAc                             <  1:1597873/150‑1 (MQ=255)
                                                                                                                        gCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCCGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCAGACTGCTGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACt                            >  1:837646/1‑150 (MQ=255)
                                                                                                                         cAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTg                           <  1:723647/150‑1 (MQ=255)
                                                                                                                             acCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                               <  2:127108/126‑1 (MQ=255)
                                                                                                                             acCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCaga                                               >  1:127108/1‑126 (MQ=255)
                                                                                                                              ccGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                                     <  1:899974/119‑1 (MQ=255)
                                                                                                                              ccGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTAtt                                                     >  2:899974/1‑119 (MQ=255)
                                                                                                                                  cgaactGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGc                                                           <  2:122134/104‑1 (MQ=38)
                                                                                                                                  cgaactGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGc                                                           >  1:122134/6‑109 (MQ=38)
                                                                                                                                    tCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCt                                      >  2:939504/1‑128 (MQ=255)
                                                                                                                                    tCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCt                                      <  1:939504/128‑1 (MQ=255)
                                                                                                                                     actGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTcc               >  2:1417625/3‑150 (MQ=255)
                                                                                                                                      aTGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                        >  2:1493956/1‑92 (MQ=255)
                                                                                                                                      aTGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATcggcgg                                                                        <  1:1493956/92‑1 (MQ=255)
                                                                                                                                       tGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg             <  1:646457/150‑1 (MQ=255)
                                                                                                                                       tGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAACTACACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAg             <  1:500740/150‑1 (MQ=255)
                                                                                                                                       tGCAACGCCAGGCAGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCACCACTGGCCTGCAACTCCAg             >  1:1082782/1‑150 (MQ=255)
                                                                                                                                        gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCTTTGCCGCCACTGGCCTGCAACTCCAGt            >  1:414709/1‑150 (MQ=255)
                                                                                                                                        gCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGt            <  2:414709/150‑1 (MQ=255)
                                                                                                                                          aaCGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCtt                                     <  2:1417158/123‑1 (MQ=255)
                                                                                                                                          aaCGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCtt                                     >  1:1417158/1‑123 (MQ=255)
                                                                                                                                            cGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTa        >  1:606820/1‑150 (MQ=255)
                                                                                                                                             gcctgcctGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTAc       <  1:1362087/150‑1 (MQ=255)
                                                                                                                                                  cctgGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGt  >  2:1369652/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACCACTGGGGCCAAGGATGCTCACGACCTCACCTGGTGCCAGCTGCAGGTCGAACTGCGCCAGCACCTGCTGCCAACCACCCTCACGCGGGTAGCCCAAGCGGATGTCACGGGCCTCCAGCAACACCGCGCTCATGCAACGCCTGCCTGGCGTTGCAGTTCGGCGCGCAATTGCACCAGGCTGGGCGTGACGATCGGCACGAAGGCCGACTCCCGCCATCGGCGGGCAAACCCTTCGCCGTATTCGCAGAGGTAAGCCTTGCCGCCACTGGCCTGCAACTCCAGTTGTACCGCGT  >  NC_002947/4133877‑4134172

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: