Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,492,062 T→A 5.7% L481Q (CTG→CAG)  PP_2185 → formate dehydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,492,0620TA 5.7% 335.3 / 7.2 104L481Q (CTG→CAG) PP_2185formate dehydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (55/43);  new base A (3/3);  total (58/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 108 total aligned reads displayed.

CAAGGGGCGCGGCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGCCCGGCTGT  >  NC_002947/2491914‑2492211
                                                                                                                                                     |                                                                                                                                                     
cAAGGGGCGGGGCTGATCGGGATCGGCCCACGGCGTATCGACCTGGGCGACTCAGCGCACGCCCGCGCCGGCCTGCACCTGCAACTGCGCCCGGGCACCCACGTGGCCATACTCAACGCCCCGGCCCACGTGATCATCACCAAGGGCctg                                                                                                                                                      >  2:89047/1‑150 (MQ=255)
cAAGGGGCGCGGCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCctg                                                                                                                                                      <  1:910484/150‑1 (MQ=255)
cAAGGGGCGCGGCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCctg                                                                                                                                                      >  1:2157714/1‑150 (MQ=255)
     ggCGCGGCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGc                                                                                                                                                 >  1:1997393/1‑150 (MQ=255)
      gcgcGGCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGcc                                                                                                                                                <  2:1245138/150‑1 (MQ=255)
          ggCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCaaa                                                                                                                                            <  1:2114882/150‑1 (MQ=255)
          ggCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCAGCAGGccctcg                                                                                                                                            <  2:45440/150‑4 (MQ=255)
              gatcgtgatcgACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCtt                                                                                                                                        >  2:2215125/1‑150 (MQ=255)
              gatcgtgatcgACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCtt                                                                                                                                        >  2:2215124/1‑150 (MQ=255)
               atcgtgatcgACCCACGGCGTATCGACCTGGTCGCCTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCAAACGGTTGTTTCCGACTTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTc                                                                                                                                       <  1:2225055/150‑1 (MQ=255)
                tcgtgatcgACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCa                                                                                                                                      >  1:518064/1‑150 (MQ=255)
                     atcgAGCCACGGCGTATCGACCTGGTCGCCTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGGCCGGGCACCAACGTGGCCAGGCTCAACCCCCTGGCCAACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGAt                                                                                                                                 <  2:867212/150‑1 (MQ=255)
                     atcgACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGAt                                                                                                                                 >  1:868148/1‑150 (MQ=255)
                         aCCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTCGCCCAAAGCTTCATCGATGacc                                                                                                                             >  2:1037782/1‑150 (MQ=255)
                             aCGGCGTATCGACCTGGTCGAGTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTg                                                                                                                         >  1:1633650/1‑150 (MQ=255)
                             aCGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTg                                                                                                                         >  1:444878/1‑150 (MQ=255)
                              cGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGc                                                                                                                        <  2:1333751/150‑1 (MQ=255)
                                  gTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGaaa                                                                                                                    >  1:190522/1‑150 (MQ=255)
                                       gACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGga                                                                                                               <  2:190522/150‑1 (MQ=255)
                                       gACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGga                                                                                                               <  1:158416/150‑1 (MQ=255)
                                       gACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTCACACGACGTCCGGGCACGAACGTGGCCATGCACAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGga                                                                                                               <  1:542212/150‑1 (MQ=255)
                                           tGGTCGACTCGCCGCTCTCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCCCATTGCCCTGAACCACGAACGCCCTACCCCACGTGAACGACAACGAGGGCATGCTGTTACAAATCGTCATCGATGCCCGCTGCGAAACCGGAGACt                                                                                                           <  1:1151574/150‑1 (MQ=255)
                                                gACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCg                                                                                                                            <  2:2271417/128‑1 (MQ=255)
                                                gACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCg                                                                                                                            >  1:2271417/1‑128 (MQ=255)
                                                 aCTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGccc                                                                                                     >  1:524069/1‑150 (MQ=255)
                                                   tCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGc                                                                                                                                                 >  1:2308345/1‑104 (MQ=255)
                                                   tCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGc                                                                                                                                                 <  2:2308345/104‑1 (MQ=255)
                                                     gccgcACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGGCTTCGCCCGTTg                                                                                                 >  1:1522578/1‑150 (MQ=255)
                                                        gcACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCACGTTGgcg                                                                                              >  1:1182261/1‑150 (MQ=255)
                                                          aCGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGcgcg                                                                                            >  1:679054/1‑150 (MQ=255)
                                                              ccGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCAGCAGGccctcgg                                                                                                                                           <  1:1273877/99‑5 (MQ=255)
                                                              ccGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCAGCAGGccctcgg                                                                                                                                           >  2:1273877/1‑95 (MQ=255)
                                                               cgcgCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTc                                                                                       <  2:1323076/150‑1 (MQ=255)
                                                                 cgcCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGt                                                                                     <  2:1962935/150‑1 (MQ=255)
                                                                   ccGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCa                                                                                   >  2:1962787/1‑150 (MQ=255)
                                                                   ccGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCa                                                                                   <  2:1085324/150‑1 (MQ=255)
                                                                    cGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAg                                                                                  <  2:626141/150‑1 (MQ=255)
                                                                     gACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCTTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGc                                                                                 <  1:1143816/150‑1 (MQ=255)
                                                                       cctgcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCt                                                                               <  2:444878/150‑1 (MQ=255)
                                                                        ctgcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTg                                                                              >  2:735421/1‑150 (MQ=255)
                                                                        ctgcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTg                                                                              >  2:1896327/1‑150 (MQ=255)
                                                                        ctgcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCCTTGTCGCGACTTCGTCAGCCTg                                                                              >  2:558873/1‑150 (MQ=255)
                                                                         tgcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGc                                                                             >  1:2314359/1‑150 (MQ=255)
                                                                          gcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGcc                                                                                >  2:1100100/1‑146 (MQ=255)
                                                                          gcacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGcc                                                                                <  1:1100100/146‑1 (MQ=255)
                                                                           cacctgcaACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCa                                                                           >  1:2092418/1‑150 (MQ=255)
                                                                                  aaCTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTg                                                                                                 <  2:468422/121‑1 (MQ=255)
                                                                                          ccGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAg                                                            <  1:2215124/150‑1 (MQ=255)
                                                                                          ccGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAg                                                            >  2:69750/1‑150 (MQ=255)
                                                                                          ccGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAg                                                            <  1:2215125/150‑1 (MQ=255)
                                                                                               caacAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGACCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCt                                                       <  2:2297339/147‑1 (MQ=255)
                                                                                                aCCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTc                                                      >  2:1210712/1‑150 (MQ=255)
                                                                                                  cAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCGCCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTAGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCgg                                                    >  1:2297561/1‑150 (MQ=255)
                                                                                                  cAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACATCGCCCTTTTGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCgg                                                    >  2:1873326/1‑150 (MQ=255)
                                                                                                    aCGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACAGCTTCAGCCTGCCAGACAACGCCCCCGAGGAGATCggaa                                                  >  1:634689/1‑148 (MQ=255)
                                                                                                     cGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGccc                                                 <  2:868148/150‑1 (MQ=255)
                                                                                                       tGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCgg                                               <  1:754805/150‑1 (MQ=255)
                                                                                                           cATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGAt                                                                                                                                 >  1:1064787/1‑64 (MQ=255)
                                                                                                           cATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGAt                                                                                                                                 <  2:1064787/64‑1 (MQ=255)
                                                                                                           cATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTtg                                           >  2:2252609/1‑150 (MQ=255)
                                                                                                           cATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTtg                                           <  2:2157714/150‑1 (MQ=255)
                                                                                                            aTGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTtgt                                          <  2:2143092/150‑1 (MQ=255)
                                                                                                             ctgttccaCGCCCTGCCCCACGTGCACGTCACCGAGCGCCTCCTGGCCCAAAGCTCCCTCGATGCCCGCTGCGAAACCGGAGACGTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCACACAACGCCCCCGAGGTGCTCGGCCCGGTTtgt                                          <  1:1794457/143‑1 (MQ=255)
                                                                                                               cTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGc                                       <  2:1624186/150‑1 (MQ=255)
                                                                                                               cTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGc                                       <  2:524069/150‑1 (MQ=255)
                                                                                                               cTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGc                                       <  1:1896327/150‑1 (MQ=255)
                                                                                                                tCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCg                                      >  2:1783663/1‑150 (MQ=255)
                                                                                                                tCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCg                                      >  2:51420/1‑150 (MQ=255)
                                                                                                                   aaCGCGCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTg                                    <  2:1075106/150‑1 (MQ=255)
                                                                                                                    aCGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGc                                                                                 <  1:1001238/104‑1 (MQ=255)
                                                                                                                       cccTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGt                                     <  1:2117899/145‑1 (MQ=255)
                                                                                                                        ccTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCgg                                               >  1:2083086/1‑134 (MQ=255)
                                                                                                                        ccTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCgg                                               <  2:2083086/134‑1 (MQ=255)
                                                                                                                        ccTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCccgc                               >  1:984052/1‑150 (MQ=255)
                                                                                                                        ccTGGCCCACGTGATCGGCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAACATGCCAGCCTGCACGCGGTAGCTGCGCGGCCCGGTTTTTGGCGTGCccgc                               >  2:1230988/1‑150 (MQ=255)
                                                                                                                           ggCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGATCGGCCCGGTTTGTGGCGTGCCCGCCGa                            <  1:1873326/150‑1 (MQ=255)
                                                                                                                            gCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAg                           >  1:1279679/1‑150 (MQ=255)
                                                                                                                            gCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAg                           >  1:1279210/1‑150 (MQ=255)
                                                                                                                              ccACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGga                                                                                                               >  1:1534809/1‑64 (MQ=255)
                                                                                                                              ccACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCa                         <  1:89047/150‑1 (MQ=255)
                                                                                                                               cACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGa                                                                                           <  1:668319/83‑1 (MQ=255)
                                                                                                                               cACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGa                                                                                           >  2:668319/1‑83 (MQ=255)
                                                                                                                                  gTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCccgccg                             >  2:1796527/1‑142 (MQ=255)
                                                                                                                                  gTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTCCGCCCGTTGGCGTGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCccgccg                             <  1:1796527/142‑1 (MQ=255)
                                                                                                                                    gATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAGCCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGCTTGTGGCGTGCCCGCCGAGCAGATCcg                   >  2:1393487/1‑150 (MQ=255)
                                                                                                                                       cGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCgcgc                <  1:1345181/150‑1 (MQ=255)
                                                                                                                                        gTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTACTTCCCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGTGTGCCAGCCGAGCAGATCCGCGCt               >  1:2332266/1‑150 (MQ=255)
                                                                                                                                         tCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTg              >  2:2274935/1‑150 (MQ=255)
                                                                                                                                         acACCGCGAGCCTGCTGGCCCAACGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTg              >  2:2088595/2‑150 (MQ=255)
                                                                                                                                            ctttgggCCAGCAGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCgg           >  2:786999/5‑150 (MQ=255)
                                                                                                                                            ccGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGTAATTCGCGGCCACGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCgg           <  2:63395/150‑1 (MQ=255)
                                                                                                                                            ccGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCgg           >  2:1665044/1‑150 (MQ=255)
                                                                                                                                            ccGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCgg           <  2:2297561/150‑1 (MQ=255)
                                                                                                                                             tttgggCCAGCAGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCgg           >  1:1882164/4‑149 (MQ=255)
                                                                                                                                             tttgggCCAGCAGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCgg           <  2:1882164/146‑1 (MQ=255)
                                                                                                                                             cGAGGGCCTGCTGGCCCAAAGCTTCATCGCTTCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGc          >  1:63395/1‑150 (MQ=255)
                                                                                                                                                 ggCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGCCCgg      >  2:1257903/1‑150 (MQ=255)
                                                                                                                                                  gCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGCCCGGc     >  1:231129/1‑150 (MQ=255)
                                                                                                                                                   ccTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGCCCGGCt    >  1:427973/1‑150 (MQ=255)
                                                                                                                                                     tgctgGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGCCCGGCTGt  >  1:22078/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGGGGCGCGGCTGATCGTGATCGACCCACGGCGTATCGACCTGGTCGACTCGCCGCACGCCCGCGCCGACCTGCACCTGCAACTGCGCCCGGGCACCAACGTGGCCATGCTCAACGCCCTGGCCCACGTGATCGTCACCGAGGGCCTGCTGGCCCAAAGCTTCATCGATGCCCGCTGCGAAACCGGAGACTTCGCCCGTTGGCGCGACTTCGTCAGCCTGCCAGACAACGCCCCCGAGGTGCTCGGCCCGGTTTGTGGCGTGCCCGCCGAGCAGATCCGCGCTGCGGCCCGGCTGT  >  NC_002947/2491914‑2492211

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