Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,074 T→A 6.2% Q48L (CAG→CTG)  PP_2201 ← alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,0740TA 6.2% 451.5 / 13.4 142Q48L (CAG→CTG) PP_2201alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (67/66);  new base A (5/4);  total (72/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

GCGCGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATGGGGGC  >  NC_002947/2507925‑2508222
                                                                                                                                                     |                                                                                                                                                    
gcgcGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCt                                                                                                                                                      >  1:2181280/1‑150 (MQ=255)
    gATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGcc                                                                                                                                                  >  1:287127/1‑150 (MQ=255)
      tGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCt                                                                                                                                                <  1:2045996/150‑1 (MQ=255)
       gACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTg                                                                                                                                               <  2:806431/150‑1 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTgg                                                                                                                                                    >  2:2217124/1‑142 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTgg                                                                                                                                                    <  1:2217124/142‑1 (MQ=255)
          ccGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGa                                                                                                                                            <  1:897077/150‑1 (MQ=255)
                  cGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGaa                                                                                                                                      >  1:687244/1‑148 (MQ=255)
                  cGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTAGAGTGCTTGGTCTTGCCTTGGGCGGGCACGTCACTGCAGACCATCTGTAGGACCTGGCCCTGGGACTTGAAGt                                                                                                                                    >  1:1252230/1‑150 (MQ=255)
                   ggCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTc                                                                                                                                   >  2:895907/1‑150 (MQ=255)
                      cTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGaa                                                                                                                                <  1:1114171/150‑1 (MQ=255)
                          cTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGc                                                                                                                            >  2:1540241/1‑150 (MQ=255)
                          cTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGc                                                                                                                            <  2:1310206/150‑1 (MQ=255)
                              gggCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGa                                                                                                                        <  1:578800/150‑1 (MQ=255)
                                gCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAgg                                                                                                                      <  2:2314094/150‑1 (MQ=255)
                                            gtTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTgg                                                                                                                                              >  1:347523/1‑114 (MQ=255)
                                            gtTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTgg                                                                                                                                              <  2:347523/114‑1 (MQ=255)
                                               tCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTg                                                                                                       >  2:928183/1‑150 (MQ=255)
                                                 ccAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGaa                                                                                                     <  2:1252230/150‑1 (MQ=255)
                                                 ccAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCATTGGGCGGGCAAGTCCATGTAGCCCATCGGCAGGTTCTGGCCTTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGaa                                                                                                     <  2:1009288/150‑1 (MQ=255)
                                                   aGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGc                                                                                                   >  1:1361533/1‑150 (MQ=255)
                                                     gTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCt                                                                                                 >  2:2337771/1‑150 (MQ=255)
                                                       ggcggcGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTc                                                                                               <  1:1051448/150‑1 (MQ=255)
                                                        gcggcgCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCggg                                                                                                             <  2:1510113/135‑1 (MQ=255)
                                                        gcggcgCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCggg                                                                                                             >  1:1510113/1‑135 (MQ=255)
                                                        gcggcgCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCg                                                                                              >  2:121361/1‑150 (MQ=255)
                                                         cggcgCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGa                                                                                             <  2:1776276/150‑1 (MQ=255)
                                                          ggcgCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTCCAGGGCCAGGCCCTgcagat                                                                                                                                          >  1:258056/1‑99 (MQ=255)
                                                             gcAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTg                                                                                         <  1:39057/150‑1 (MQ=255)
                                                              cAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                        <  2:1141229/150‑1 (MQ=255)
                                                               aGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTggg                                                                                       >  1:1742163/1‑150 (MQ=255)
                                                                gAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgggg                                                                                      <  1:1859606/150‑1 (MQ=255)
                                                                 aaGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTggggg                                                                                     <  1:1585003/150‑1 (MQ=255)
                                                                    ttCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCg                                                                                  >  2:89160/1‑150 (MQ=255)
                                                                       ttGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgtag                                                                               <  2:525952/150‑1 (MQ=255)
                                                                           cGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgtaggtag                                                                           >  2:1858930/1‑150 (MQ=255)
                                                                             tGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGc                                                                         <  2:1576308/150‑1 (MQ=255)
                                                                             tGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGc                                                                         <  1:1540241/150‑1 (MQ=255)
                                                                                aTCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                        >  1:1805806/1‑132 (MQ=255)
                                                                                aTCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTgg                                                                                        <  2:1805806/132‑1 (MQ=255)
                                                                                    gCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCgg                                                                  <  1:895907/150‑1 (MQ=255)
                                                                                     caccacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGc                                                                 >  2:1912682/1‑150 (MQ=255)
                                                                                       ccacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAg                                                               >  1:847972/1‑150 (MQ=255)
                                                                                        cacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCgg                                                                                                              >  2:1959021/1‑102 (MQ=255)
                                                                                        cacGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCgg                                                                                                              <  1:1959021/102‑1 (MQ=255)
                                                                                        cacGTTGTGGCCGGTGGCCTGGCCTTGGGCGGGCCCGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGc                                                              >  2:735059/1‑150 (MQ=255)
                                                                                            ttGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATg                                                          <  2:1587820/150‑1 (MQ=255)
                                                                                              gtgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGc                                                                         <  1:1367889/133‑1 (MQ=255)
                                                                                              gtgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGc                                                                         >  2:1367889/1‑133 (MQ=255)
                                                                                               tgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAgggg                                                                                                                    <  1:935043/89‑1 (MQ=255)
                                                                                               tgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGcc                                                       >  1:1583144/1‑150 (MQ=255)
                                                                                               tgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGcc                                                       >  2:1843302/1‑150 (MQ=255)
                                                                                               tgGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGcc                                                       >  1:178001/1‑150 (MQ=255)
                                                                                                 gCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGa                                                     >  1:843287/1‑150 (MQ=255)
                                                                                                   cGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATg                                                          <  2:333396/143‑1 (MQ=255)
                                                                                                    gTTGGCCTGGCCTTGGGCGAGCACGTCCATGTAGCCCATCTGCAGGACCTGGACCTGGGACTTGAAGTCGAAATGCTTGAGGGGGAGCGGGAAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAAct                                                  >  2:1838342/1‑150 (MQ=255)
                                                                                                      tggcctggccTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgc                                                <  2:287127/150‑1 (MQ=255)
                                                                                                      tggcctggccTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTCGGCCGGGGACTTTAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgc                                                >  1:60581/1‑150 (MQ=255)
                                                                                                        gcctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgctg                                              <  1:1843036/150‑1 (MQ=255)
                                                                                                        gcctggcctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgctg                                              <  2:1722492/150‑1 (MQ=255)
                                                                                                              cctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAgg                                                                                                                      <  2:1268304/72‑1 (MQ=255)
                                                                                                              cctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAgg                                                                                                                      >  1:1268304/1‑72 (MQ=255)
                                                                                                              cctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAagg                                        <  1:121361/150‑1 (MQ=255)
                                                                                                              cctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGATGAAGCCTTCGAGCAGTGGGCGGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAagg                                        >  1:2343359/1‑150 (MQ=255)
                                                                                                               ctTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAaggc                                       >  1:667363/1‑150 (MQ=255)
                                                                                                                 tGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAaggcag                                     >  2:1461356/1‑150 (MQ=255)
                                                                                                                       ggcACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAgg                               <  2:843287/150‑1 (MQ=255)
                                                                                                                          aCGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCa                            >  2:1876559/1‑150 (MQ=255)
                                                                                                                          aCGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCa                            >  1:942877/1‑150 (MQ=255)
                                                                                                                           cGTCCATTTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGCGTAGCTGAAGCCGTCGCGCTGGGGGACGTACGTAGGCGACACGGAAGCGATGGCCCACCTGCTGGCAAGGCAGGCTAgcgcta                           <  2:31533/150‑3 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAg                                                                                                          <  2:483350/70‑1 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAg                                                                                                          >  1:483350/1‑70 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAg                          >  2:1755897/1‑150 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAg                          <  1:1755897/150‑1 (MQ=255)
                                                                                                                            gTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAg                          >  1:1965129/1‑150 (MQ=255)
                                                                                                                               cATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCggg                                                                                                             <  1:991413/64‑1 (MQ=255)
                                                                                                                                aTGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCgtaggta                                                                            <  1:2107077/96‑1 (MQ=255)
                                                                                                                                  gTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGa                                                            >  1:1384442/1‑110 (MQ=255)
                                                                                                                                  gTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGa                                                            <  2:1384442/110‑1 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGTCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAActgctg                                              >  2:2024883/1‑119 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGcc                                                       <  2:1035335/110‑1 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGcc                                                       >  1:1035335/1‑110 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc               >  2:931132/1‑150 (MQ=255)
                                                                                                                                       ccATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc               >  1:799827/1‑150 (MQ=255)
                                                                                                                                        cATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgcg              >  2:2227564/1‑150 (MQ=255)
                                                                                                                                          tCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgca            <  2:942877/150‑1 (MQ=255)
                                                                                                                                          gtcccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc               >  2:155763/5‑147 (MQ=255)
                                                                                                                                          gtcccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgc               <  1:155763/143‑1 (MQ=255)
                                                                                                                                          gtcccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgca            >  1:481902/5‑150 (MQ=255)
                                                                                                                                          gtcccAGGGCCAGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCGGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTgcgcgag            >  1:325052/5‑148 (MQ=255)
                                                                                                                                          gtcccAGGGCCAGGCCCTGAGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAagg                                        >  2:2120876/5‑122 (MQ=255)
                                                                                                                                            tGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCgcagc          >  1:2334529/1‑150 (MQ=255)
                                                                                                                                              cAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCAt        <  1:1876559/150‑1 (MQ=255)
                                                                                                                                               aGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCAta       >  2:2177538/1‑149 (MQ=255)
                                                                                                                                               aGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATg       >  2:935114/1‑150 (MQ=255)
                                                                                                                                               aGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGATTCGCGCAGCATg       >  1:1099661/1‑150 (MQ=255)
                                                                                                                                                  gCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATgggg    >  1:1463543/1‑150 (MQ=255)
                                                                                                                                                   ccTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATggggg   <  2:2334529/150‑1 (MQ=255)
                                                                                                                                                    cTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATgggg    <  1:2066261/148‑1 (MQ=255)
                                                                                                                                                    cTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATgggg    >  2:2066261/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCGATGACCCGGTAGCCGGCCTTGCTCAGGGCGTCGATGGTGGTTTCCCAGGTGGCGGCGCAGAAGTTCTTGCCGTGCATCAGCACCACGTTGTGGCCGTTGGCCTGGCCTTGGGCGGGCACGTCCATGTAGCCCATCTGCAGGGCCTGGCCCTGGGACTTGAAGTCGAAATGCTTGAGGGGGTGCGGGTAGCTGAAGCCTTCGAGCTGGGGGCCGTAGGTAGGCGACTCGGCAGCGATGGCCGAACTGCTGGCAAGGCAGGCCAGGGCAAGGGCGGTTGCGCGCAGCATGGGGGC  >  NC_002947/2507925‑2508222

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: