Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 873,425 T→C 8.2% L20P (CTT→CCT)  PP_5457 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947873,4250TC 8.2% 236.4 / 9.9 85L20P (CTT→CCT) PP_5457hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base C (5/2);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAATGAT  >  NC_002947/873280‑873570
                                                                                                                                                 |                                                                                                                                                  
ccAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGc                                                                                                                                                 >  1:397031/1‑149 (MQ=255)
 cAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGcc                                                                                                                                                <  5:113614/149‑1 (MQ=255)
 cAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGc                                                                                                                                                 >  1:401445/1‑148 (MQ=255)
    gtgtGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGc                                                                                                                                             <  6:102542/149‑1 (MQ=255)
    gtgtGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTGGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCa                                                                                                                                               <  6:332732/147‑1 (MQ=255)
         cTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAgg                                                                                                                                         <  5:202104/148‑1 (MQ=255)
           gTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                      >  1:249328/1‑149 (MQ=255)
            tCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGa                                                                                                                                      <  3:69849/148‑1 (MQ=255)
                ccTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                                 <  5:368957/149‑1 (MQ=255)
                ccTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGAt                                                                                                                                 <  8:265703/149‑1 (MQ=255)
                    tggtgTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTac                                                                                                                             <  2:198133/149‑1 (MQ=255)
                    tggtgTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTa                                                                                                                              <  2:43323/148‑1 (MQ=255)
                         tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTg                                                                                                                        >  8:272990/1‑149 (MQ=255)
                         tgCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTg                                                                                                                        >  3:354802/1‑149 (MQ=255)
                                  cAGGTGCATGGTGACGAATGGTTGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                               <  8:226617/149‑1 (MQ=255)
                                        cATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCtt                                                                                                                                                   >  3:344154/1‑107 (MQ=255)
                                          tGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGt                                                                                                       <  2:249328/149‑1 (MQ=255)
                                          tGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGt                                                                                                       <  2:202049/149‑1 (MQ=255)
                                              gACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                               <  3:165889/137‑1 (MQ=255)
                                              gACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGt                                                                                                               >  4:165889/1‑137 (MQ=255)
                                                 gAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgcc                                                                                                >  8:140641/1‑149 (MQ=255)
                                                  aaTGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgccg                                                                                               <  6:295998/149‑1 (MQ=255)
                                                    tGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTgccgcc                                                                                             >  4:11395/1‑149 (MQ=255)
                                                      gATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCgg                                                                                           >  8:335340/1‑149 (MQ=255)
                                                           ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGc                                                                                      >  4:26830/1‑149 (MQ=255)
                                                           ttGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGc                                                                                      >  2:31926/1‑149 (MQ=255)
                                                              aaGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCTGTGCCGCCGGCGAGCTg                                                                                    <  4:354802/148‑1 (MQ=255)
                                                                 tCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGATCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCt                                                                                >  1:292210/1‑149 (MQ=255)
                                                                    gccgcGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATATACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑Gc                                                                            <  4:46763/149‑1 (MQ=255)
                                                                     ccgcGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑g                                                                             >  6:49321/1‑147 (MQ=255)
                                                                       gcgAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAgg                                                                         <  7:272990/149‑1 (MQ=255)
                                                                        cgAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAggg                                                                        <  8:302568/149‑1 (MQ=255)
                                                                         gAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGc                                                                       <  2:397031/149‑1 (MQ=255)
                                                                            gacgacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCctg                                                                    >  4:81898/1‑149 (MQ=255)
                                                                             acgacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCctgc                                                                   >  8:341051/1‑149 (MQ=255)
                                                                              cgacACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCctgc                                                                   <  7:56262/148‑1 (MQ=255)
                                                                                   ccTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGc                                                             <  7:341051/149‑1 (MQ=255)
                                                                                    cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGc                                                             >  4:18456/1‑148 (MQ=255)
                                                                                    cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGc                                                             <  3:18456/148‑1 (MQ=255)
                                                                                    cTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCt                                                            >  3:1838/1‑149 (MQ=255)
                                                                                       aTGCAGCGTATCTACCAACAGTAAAGCCTTCTTGAGGCGCAAATGCTGTCAGGCATGCTTGCCAGCGAGGGGATCGATTTACACGTGGTCGTCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCCA‑GCAGGCCCTGCAAGGGcagc                                                          <  3:33460/148‑4 (MQ=255)
                                                                                         gcggCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCtcgtcg                                                       >  7:312948/4‑149 (MQ=255)
                                                                                         gCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCtcgtcg                                                       >  6:87225/1‑149 (MQ=255)
                                                                                           aGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTg                                                     >  5:391442/1‑149 (MQ=255)
                                                                                            gCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTgg                                                    >  2:287352/1‑149 (MQ=255)
                                                                                            gCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTgg                                                    <  3:26830/149‑1 (MQ=255)
                                                                                               tATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATCTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCtg                                                 <  3:11395/149‑1 (MQ=255)
                                                                                                aTCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCtga                                                >  3:291693/1‑149 (MQ=255)
                                                                                                aTCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTCGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCtga                                                >  1:127510/1‑149 (MQ=255)
                                                                                                 tCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTTCGCAGGGCCTGCTGGGGCTCGTCGTGGCtg                                                 <  1:237148/148‑1 (MQ=255)
                                                                                                    aCGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAg                                            >  5:270534/1‑149 (MQ=255)
                                                                                                     cGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCtgatga                                             <  1:31926/147‑1 (MQ=255)
                                                                                                       aaCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAg                                         >  7:254899/1‑149 (MQ=255)
                                                                                                       aaCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGATCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAg                                         <  2:292210/149‑1 (MQ=255)
                                                                                                        aCCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAgg                                        >  5:36103/1‑149 (MQ=255)
                                                                                                          cGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTgg                                                    >  3:68787/1‑135 (MQ=255)
                                                                                                          cGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTgg                                                    <  4:68787/135‑1 (MQ=255)
                                                                                                            gAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCTTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGa                                    <  6:270534/149‑1 (MQ=255)
                                                                                                            gAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCg                                     >  2:5181/1‑148 (MQ=255)
                                                                                                              aaGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCtgat                                               <  5:223538/136‑1 (MQ=255)
                                                                                                              aaGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCtgat                                               >  6:223538/1‑136 (MQ=255)
                                                                                                               aGCCTGCTTGAGGCGCATATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCgg                                                                                           >  5:368606/1‑92 (MQ=255)
                                                                                                               aGCCTGCTTGAGGCGCATATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCgg                                                                                           <  6:368606/92‑1 (MQ=255)
                                                                                                                        gAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTgg                                                    >  7:176588/1‑121 (MQ=255)
                                                                                                                        gAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTgg                                                    <  8:176588/121‑1 (MQ=255)
                                                                                                                         aGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAg                       >  4:82473/1‑149 (MQ=255)
                                                                                                                             gcAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGa                   >  1:14533/1‑149 (MQ=255)
                                                                                                                              cAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGAt                  <  7:335340/149‑1 (MQ=255)
                                                                                                                              cAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGAt                  <  5:87225/149‑1 (MQ=255)
                                                                                                                               aaaTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGc                                                                                          >  6:85486/1‑77 (MQ=255)
                                                                                                                               aaaTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGc                                                                                          <  5:85486/77‑1 (MQ=255)
                                                                                                                               aaaTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc                 >  5:321567/1‑149 (MQ=255)
                                                                                                                                ctcgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc                 >  5:44518/4‑148 (MQ=255)
                                                                                                                                aaTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCCTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc                 >  3:211944/1‑148 (MQ=255)
                                                                                                                                aaTGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCATGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATc                 >  3:206534/1‑148 (MQ=255)
                                                                                                                                  cgCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCGTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCCCACGGCAGCTGCTCGAt              >  7:257025/2‑149 (MQ=255)
                                                                                                                                   gCTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATg             >  4:303477/1‑149 (MQ=255)
                                                                                                                                    cTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGc            <  4:291693/149‑1 (MQ=255)
                                                                                                                                    cTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTgggg                                                              >  4:332052/1‑99 (MQ=255)
                                                                                                                                    cTGGCAAGCATGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTgggg                                                              <  3:332052/99‑1 (MQ=255)
                                                                                                                                       gcaggcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGt          <  3:81898/148‑1 (MQ=255)
                                                                                                                                         aggcaAGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTAc        >  5:251680/1‑148 (MQ=255)
                                                                                                                                           gcaTGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAat     >  6:188038/1‑149 (MQ=255)
                                                                                                                                           gcaTGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAat     >  5:262039/1‑149 (MQ=255)
                                                                                                                                           gcaTGCCTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAat     <  6:44518/149‑1 (MQ=255)
                                                                                                                                              tGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAatgat  <  6:36103/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCAGGTGTGCTGTCAGCCTGTGGTGTGCATTCTGCAGGTGCATGGTGACGAATGGATGCTTGAAGTCCGCCGCGAGGACGACACCTGATGCAGCGTATCTACGAACCGGAAAGCCTGCTTGAGGCGCAAATGCTGGCAGGCATGCTTGCCAGCGAGGGGATCGATGTACACGTGGTCGGCCGTGACCTTGTCGGTGCCGCCGGCGAGCTGCCCTT‑GCAGGGCCTGCTGGGGCTCGTCGTGGCTGATGAGCAGGCCGAATGCGCACGGCAGCTGATCGATGCGTACAATGAT  >  NC_002947/873280‑873570

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