Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,333,762:1 +A 7.5% coding (82/1155 nt) PP_1159 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,333,7621.A 7.5% 412.6 / 27.0 106coding (82/1155 nt)PP_1159membrane protein
Reads supporting (aligned to +/- strand):  ref base . (48/50);  new base A (4/4);  total (52/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 107 total aligned reads displayed.

CGACCTCTACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTTGCTGCTGGGGGCGTG  >  NC_002947/1333616‑1333910
                                                                                                                                                   |                                                                                                                                                    
cGACCTCTACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtg                                                                                                                                                   >  1:1321927/1‑150 (MQ=255)
cGACCTCTACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtg                                                                                                                                                   <  2:1623695/150‑1 (MQ=255)
 gACCTCTACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtgg                                                                                                                                                  <  1:1986570/150‑1 (MQ=255)
  aCCTCTACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggc                                                                                                                                                 <  2:1067793/150‑1 (MQ=255)
     tctACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATtggcca                                                                                                                                               >  2:1932353/1‑150 (MQ=255)
     tctACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccat                                                                                                                                              >  1:89430/1‑150 (MQ=255)
     tctACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccat                                                                                                                                              <  2:2224425/150‑1 (MQ=255)
      ctACGATGTCGTCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccatg                                                                                                                                             <  1:371389/150‑1 (MQ=255)
       tACGATGTCCCCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccatgg                                                                                                                                            <  2:538481/150‑1 (MQ=255)
       tACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCCAGACTtggccatgg                                                                                                                                            >  2:1377701/1‑150 (MQ=255)
       tACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccatgg                                                                                                                                            >  1:1360500/1‑150 (MQ=255)
       tACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccatgg                                                                                                                                            >  1:1366510/1‑150 (MQ=255)
       tACGATGTCACCCTTGGCTATCGCTTCCCTTCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTtggccatgg                                                                                                                                            >  1:65121/1‑150 (MQ=255)
             gTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCaagg                                                                                                                                      <  1:968156/150‑3 (MQ=255)
              tCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATTGGCCATGGCCaagg                                                                                                                                      <  1:1932353/150‑3 (MQ=255)
              tCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGATAGCGTTGCTTTTCTTAAGTATTCTTGTGTTTCTTGCTGCAAGCCTTGGCCATGGCCAATgg                                                                                                                                     >  2:2201287/1‑150 (MQ=255)
                  ccTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAAt                                                                                                                                       <  1:2037836/144‑1 (MQ=255)
                  ccTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAAt                                                                                                                                       >  2:2037836/1‑144 (MQ=255)
                  ccTTGGCTATCGCTTGCATGCACCTATGCCACGAGGATTTGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTgg                                                                                                                                 >  1:1771155/1‑150 (MQ=255)
                       gCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTAt                                                                                                                            <  1:1449697/150‑1 (MQ=255)
                       gCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTAGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTAt                                                                                                                            <  1:1377701/150‑1 (MQ=255)
                        cTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCTATGGCCAATGGCTGGATTATg                                                                                                                           >  1:1545874/1‑150 (MQ=255)
                            cGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGcaca                                                                                                                       <  2:1360500/150‑1 (MQ=255)
                             gCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACag                                                                                                                      >  2:718386/1‑150 (MQ=255)
                                   atgcaCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCt                                                                                                                                   <  1:876118/131‑1 (MQ=255)
                                   atgcaCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCt                                                                                                                                   >  2:876118/1‑131 (MQ=255)
                                   atgcaCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACagag                                                                                                                    <  2:1913912/146‑1 (MQ=255)
                                   atgcaCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACagag                                                                                                                    >  1:1913912/1‑146 (MQ=255)
                                   atgcaCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTg                                                                                                                >  2:2284964/1‑150 (MQ=255)
                                    tgcaCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTgg                                                                                                               >  2:2195755/1‑150 (MQ=255)
                                       aCCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGa                                                                                                            >  1:1059027/1‑150 (MQ=255)
                                              ccACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTtga                                                                                                     >  2:1418108/1‑150 (MQ=255)
                                              ccACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTtga                                                                                                     >  1:985260/1‑150 (MQ=255)
                                               cACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTtgat                                                                                                    <  1:2291199/150‑1 (MQ=255)
                                                  gTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACagag                                                                                                                    >  2:1520289/1‑131 (MQ=255)
                                                  gTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACagag                                                                                                                    <  1:1520289/131‑1 (MQ=255)
                                                    tGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGcc                                                                                               <  2:915951/150‑1 (MQ=255)
                                                                       acGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCAGAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTgg                                                                            >  1:1944815/1‑150 (MQ=255)
                                                                        cGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGa                                                                           <  2:1373904/150‑1 (MQ=255)
                                                                        cGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGa                                                                           <  2:985260/150‑1 (MQ=255)
                                                                        cGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGa                                                                           <  1:2195755/150‑1 (MQ=255)
                                                                         gCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTtga                                                                                                     >  1:1756067/1‑123 (MQ=255)
                                                                         gCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTtga                                                                                                     <  2:1756067/123‑1 (MQ=255)
                                                                           tcttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCtt                                                                                                       >  1:113924/1‑119 (MQ=255)
                                                                           tcttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCtt                                                                                                       <  2:113924/119‑1 (MQ=255)
                                                                             ttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGa                                                                                                                                <  1:346360/92‑1 (MQ=255)
                                                                             ttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGa                                                                                                                                >  2:346360/1‑92 (MQ=255)
                                                                             ttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACagag                                                                                                                    >  2:338239/1‑104 (MQ=255)
                                                                             ttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACagag                                                                                                                    <  1:338239/104‑1 (MQ=255)
                                                                             ttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCCGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGca                                                                                                                         >  1:1666121/1‑99 (MQ=255)
                                                                             ttcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCCGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGca                                                                                                                         <  2:1666121/99‑1 (MQ=255)
                                                                              tcttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGAttgttg                                                                     <  1:2201287/150‑1 (MQ=255)
                                                                               cttctCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTgg                                                                    <  2:1545874/150‑1 (MQ=255)
                                                                                   tcACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATTGGCCATGGCCaaggc                                                                                                                                     <  2:1907535/82‑4 (MQ=39)
                                                                                   tcACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATTGGCCATGGCCaaggc                                                                                                                                     >  1:1907535/1‑79 (MQ=39)
                                                                                   tcACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGat                                                                >  1:2279575/1‑150 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATtggccatggccagggc                                                                                                                                     >  2:502438/1‑76 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATtggccatggccagggc                                                                                                                                     <  1:502438/80‑5 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGAttgtt                                                                      <  2:948668/142‑1 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGAttgtt                                                                      >  1:948668/1‑142 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGatat                                                              <  1:1788822/150‑1 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGGTGGTGatat                                                              >  1:2111954/1‑150 (MQ=255)
                                                                                     aCGACAGCAGTATTGGCTTGCGTTACTTGTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGatat                                                              <  1:1840/150‑1 (MQ=255)
                                                                                      cGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGc                                                                                                >  2:2136881/1‑115 (MQ=255)
                                                                                      cGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGc                                                                                                <  1:2136881/115‑1 (MQ=255)
                                                                                       gACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATtggccat                                                                                                                                              >  2:622491/1‑69 (MQ=255)
                                                                                       gACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCATtggccat                                                                                                                                              <  1:622491/69‑1 (MQ=255)
                                                                                        aCAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAg                                                           <  1:1418108/150‑1 (MQ=255)
                                                                                            cagTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGAATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGaa                                                       <  1:99026/150‑1 (MQ=255)
                                                                                                aTTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTgg                                                   <  1:1777345/150‑1 (MQ=255)
                                                                                                 ttGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGc                                                  >  1:2247025/1‑150 (MQ=255)
                                                                                                     cTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGAAATGGTTAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTc                                              <  2:65121/150‑1 (MQ=255)
                                                                                                      ttGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGa                                                                                                            >  1:979554/1‑87 (MQ=255)
                                                                                                      ttGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGa                                                                                                            <  2:979554/87‑1 (MQ=255)
                                                                                                            tGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGcc                                                                                       <  2:1670301/102‑1 (MQ=255)
                                                                                                            tGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGcc                                                                                       >  1:1670301/1‑102 (MQ=255)
                                                                                                              aTTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCAc                                     <  1:1866709/150‑1 (MQ=255)
                                                                                                                       tgcctgcTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGat                                                                >  1:1504824/1‑114 (MQ=255)
                                                                                                                       tgcctgcTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGat                                                                <  2:1504824/114‑1 (MQ=255)
                                                                                                                       tgcctgcTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGccc                            >  2:609540/1‑150 (MQ=255)
                                                                                                                           tgctgGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCgcgc                        >  2:107066/1‑150 (MQ=255)
                                                                                                                             ctgGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGt                                                                                                                 <  1:1780858/59‑1 (MQ=255)
                                                                                                                             ctgGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGt                                                                                                                 >  2:1780858/1‑59 (MQ=255)
                                                                                                                             ctgGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTg                                                    >  1:1631150/1‑120 (MQ=255)
                                                                                                                             ctgGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTg                                                    <  2:1631150/120‑1 (MQ=255)
                                                                                                                                 tgtTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCt                  <  2:2247025/150‑1 (MQ=255)
                                                                                                                                    ttCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTtgc               <  1:2284964/150‑1 (MQ=255)
                                                                                                                                      ctcGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTg                                                                  >  2:2225047/1‑97 (MQ=255)
                                                                                                                                      ctcGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTg                                                                  <  1:2225047/97‑1 (MQ=255)
                                                                                                                                      ctcGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTATACAAGCACGAACTGCCCGCGCTGGGCTtgatg             >  2:267650/1‑150 (MQ=255)
                                                                                                                                       tcGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCttt                                               >  2:682178/1‑115 (MQ=255)
                                                                                                                                       tcGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCttt                                               <  1:682178/115‑1 (MQ=255)
                                                                                                                                       tcGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTtgctgc            >  1:1747614/1‑150 (MQ=255)
                                                                                                                                         gCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGca                                                                                                                         <  1:1689233/39‑1 (MQ=255)
                                                                                                                                         gCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGca                                                                                                                         >  2:1689233/1‑39 (MQ=255)
                                                                                                                                         gCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTTGctgctg          <  2:1747614/150‑1 (MQ=255)
                                                                                                                                             cAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGAt                                                                          >  1:9275/1‑82 (MQ=255)
                                                                                                                                             cAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGAt                                                                          <  2:9275/82‑1 (MQ=255)
                                                                                                                                              aaGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTTGCTGCTGGGGGc     >  1:238581/1‑150 (MQ=255)
                                                                                                                                                 ccTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTTGCTGCTGGGGGCGTg  >  2:548241/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGACCTCTACGATGTCACCCTTGGCTATCGCTTGCATGCACCTATGCCACGTTGAGGAGTTGCGAATGAACACGCTCTTCTTCTCACGACAGCAGTATTGGCTTGCGTTGATTTTCGGCTGCCTGCTGGTGTTTCTCGCTGCAAGCCTTGGCCATGGCCAATGGCTGGATTATGCACAGAGGGTGGCGACGCTTGATGAGCCGCTCAGCCGCTTGCGCTGGATTGTTGGTGATATCAGTGAAGTGGCTTTCTACAAGCACGAACTGCCCGCGCTGGGCTTGCTGCTGGGGGCGTG  >  NC_002947/1333616‑1333910

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: