Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,620,134 G→A 5.3% L88F (CTT→TTT)  PP_3190 ← ornithine cyclodeaminase/mu‑crystallin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,620,1340GA 5.3% 220.2 / 2.2 75L88F (CTT→TTT) PP_3190ornithine cyclodeaminase/mu‑crystallin family protein
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base A (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

ATGCCAATGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGC  >  NC_002947/3619998‑3620279
                                                                                                                                        |                                                                                                                                                 
aTGCCAATGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCcatca                                                                                                                                      <  2:85167/150‑1 (MQ=255)
 tGCCAATGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCcatcat                                                                                                                                     <  2:1326472/150‑1 (MQ=255)
  gCCAATGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCAtt                                                                                                                                    <  1:685388/150‑1 (MQ=255)
    cAATGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCCGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGc                                                                                                                                  >  1:841334/1‑150 (MQ=255)
     aaTGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCt                                                                                                                                 >  1:1756980/1‑150 (MQ=255)
           cTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAg                                                                                                                           >  1:2088474/1‑150 (MQ=255)
                ccTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCt                                                                                                                      <  1:2309857/150‑1 (MQ=255)
                 cTGGGATGGCGCCAGCAGCCGGGCGGCCATACGCCCGGCCAACGCGGTGCGCCTGGCCGTCGGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCGtt                                                                                                                     >  1:996756/1‑150 (MQ=255)
                   gggATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGTGCCTGCAGCAGTGACAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTgg                                                                                                                   >  1:755490/1‑150 (MQ=255)
                       tGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               <  2:1756980/150‑1 (MQ=255)
                       tGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               <  1:310822/150‑1 (MQ=255)
                       tGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGACTATGGCAGTCAGCCAGTCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               <  1:502522/150‑1 (MQ=255)
                        ggCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCAAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACggg                                                                                                              <  2:1769048/150‑1 (MQ=255)
                           gcCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTg                                                                                                           <  2:2140823/150‑1 (MQ=255)
                           gcCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTg                                                                                                           >  1:1337362/1‑150 (MQ=255)
                           gcCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTg                                                                                                           >  2:2143121/1‑150 (MQ=255)
                                   ccGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGaa                                                                                                   <  1:2143121/150‑1 (MQ=255)
                                     gggCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGc                                                                                                 <  2:1337362/150‑1 (MQ=255)
                                         ggcgATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGt                                                                                             <  1:1166326/150‑1 (MQ=255)
                                                 gCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCt                                                                                     >  2:1340551/1‑150 (MQ=255)
                                                 gCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCt                                                                                     >  2:1361677/1‑150 (MQ=255)
                                                  gatctccAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCTGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGTCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCtt                                                                                    <  1:1894784/145‑1 (MQ=255)
                                                  cccGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCtt                                                                                    <  2:1642125/150‑1 (MQ=255)
                                                    cGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGa                                                                                  <  1:1428034/150‑1 (MQ=255)
                                                     ggCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGa                                                                                                                  <  2:2296721/117‑1 (MQ=255)
                                                     ggCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGa                                                                                                                  >  1:2296721/1‑117 (MQ=255)
                                                      gCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               <  2:146905/119‑1 (MQ=255)
                                                      gCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               >  1:146905/1‑119 (MQ=255)
                                                        cAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCCGGCCCTTGGACGGGCTGTCCTAGAACCCTGCCGAAAGCTTGagatc                                                                              >  2:1894784/1‑146 (MQ=255)
                                                         aaCGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGTAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTg                                                                             >  1:1737999/1‑150 (MQ=255)
                                                           cgcgGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGaa                                                                           <  1:1216158/150‑1 (MQ=255)
                                                            gcgGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAg                                                                          >  2:1636730/1‑150 (MQ=255)
                                                                   gcATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGc                                                                   >  1:2182998/1‑150 (MQ=255)
                                                                            tCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCgg                                                                                                                              <  2:1723562/82‑1 (MQ=255)
                                                                            tCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCgg                                                                                                                              >  1:1723562/1‑82 (MQ=255)
                                                                            tCAGCCAGCCTTGGTCCTGCAGCAGCGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGa                                                          >  2:727386/1‑150 (MQ=255)
                                                                               gccagccTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTc                                                                                            >  2:2247725/1‑113 (MQ=255)
                                                                               gccagccTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTc                                                                                            <  1:2247725/113‑1 (MQ=255)
                                                                               gccagccTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGt                                                       >  1:946432/1‑150 (MQ=255)
                                                                                 cagccTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGCGAAGGTCGGGCTGCCTTCGATGTAg                                                     <  2:784258/150‑1 (MQ=255)
                                                                                  agccTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAgg                                                    >  2:665651/1‑150 (MQ=255)
                                                                                     cTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGAc                                                                                                                 <  2:1810547/86‑1 (MQ=255)
                                                                                     cTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGAc                                                                                                                 >  1:1810547/1‑86 (MQ=255)
                                                                                            cTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAg                                                     >  1:1396105/1‑139 (MQ=255)
                                                                                            cTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAg                                                     <  2:1396105/139‑1 (MQ=255)
                                                                                             tGCAGCAGTGCCAGCGGCTGACCGGTGTGTGTAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACg                                         >  2:1439294/1‑150 (MQ=255)
                                                                                              gcagcaGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGaa                                                                           <  1:1703112/115‑1 (MQ=255)
                                                                                              gcagcaGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGaa                                                                           >  2:1703112/1‑115 (MQ=255)
                                                                                                agcagTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCaa                                      <  1:363071/150‑1 (MQ=255)
                                                                                                  cagTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               <  2:259547/75‑1 (MQ=255)
                                                                                                  cagTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACgg                                                                                                               >  1:259547/1‑75 (MQ=255)
                                                                                                      gCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCATGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCa                                <  1:1195715/150‑1 (MQ=255)
                                                                                                      gCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGCCGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCACCAGTCa                                >  1:2282590/1‑150 (MQ=255)
                                                                                                       ccAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCAc                               <  2:2104524/150‑1 (MQ=255)
                                                                                                        cAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCAcc                              >  2:632723/1‑150 (MQ=255)
                                                                                                         aGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAg                                                     >  2:2173389/1‑126 (MQ=255)
                                                                                                         aGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAg                                                     <  1:2173389/126‑1 (MQ=255)
                                                                                                            ggCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTc                          <  1:1381301/150‑1 (MQ=255)
                                                                                                             gCTGACCGGTGTGTGCAGACAACACAAACATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACa                                    <  2:1558994/139‑1 (MQ=255)
                                                                                                             gCTGACCGGTGTGTGCAGACAACACAAACATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACa                                    <  2:1558011/139‑1 (MQ=255)
                                                                                                             gCTGACCGGTGTGTGCAGACAACACAAACATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACa                                    >  1:1558011/1‑139 (MQ=255)
                                                                                                             gCTGACCGGTGTGTGCAGACAACACAAACATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACa                                    >  1:1558994/1‑139 (MQ=255)
                                                                                                                 aCCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGAc                                                                                                                 >  2:1654251/1‑58 (MQ=255)
                                                                                                                 aCCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGAc                                                                                                                 <  1:1654251/58‑1 (MQ=255)
                                                                                                                   cGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCt                   >  2:315673/1‑150 (MQ=255)
                                                                                                                          tgCAGACAGCACAAGAATCAAGGCATCATTGCGCGGCAGGCACTGTGACGGATTGTCGTAGAAGCCTGTCGAAAGCTTGACGGGGAAGGTCGGGCTGCCTTCGATGTAGGCCGACGTAACGCAACAGTAACCGTTCGCGCCCTTGAATGc            >  2:914533/1‑150 (MQ=255)
                                                                                                                          tgCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGc            >  1:217934/1‑150 (MQ=255)
                                                                                                                              gACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTc                          >  2:25709/1‑132 (MQ=255)
                                                                                                                              gACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTc                          <  1:25709/132‑1 (MQ=255)
                                                                                                                               acaacaCAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGc       >  2:1926919/1‑150 (MQ=255)
                                                                                                                                 aacaCAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGGCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCGGACTTCGCGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGccc     >  2:1310918/1‑150 (MQ=255)
                                                                                                                                  acacAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCt    <  2:755490/150‑1 (MQ=255)
                                                                                                                                    acaAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGc  >  1:360922/1‑150 (MQ=255)
                                                                                                                                      aaGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCCTTc                          <  2:2287095/124‑1 (MQ=255)
                                                                                                                                      aaGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCCTTc                          >  1:2287095/1‑124 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
ATGCCAATGGCCTTGACCTGGGATGGCGCCAGCAGCCGGGCGGCGATACGCCCGGCCAACGCGGTGCGCATGGCAGTCAGCCAGCCTTGGTCCTGCAGCAGTGCCAGCGGCTGACCGGTGTGTGCAGACAACACAAGCATCAAGCCATCATTGCTCGGCAGGCCCTTGGACGGGTTGTCGTAGAAGCCTGTCGAAAGCTTGACGGTGAAGGTCGGGCTGCCTTCGATGTAGGCCGACTTCACGCAACAGTCACCGTTCGCGCCCTTGAATGCGAAGCCCTGC  >  NC_002947/3619998‑3620279

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