Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 346,120 T→C 9.1% N576D (AAC→GAC) ‡ PP_0287 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947346,1200TC 9.1% 280.9 / 14.6 99N576D (AAC→GAC) ‡PP_0287hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/48);  new base C (7/2);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 9.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCG  >  NC_002947/345974‑346266
                                                                                                                                                  |                                                                                                                                                  
gCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGa                                                                                                                                                  >  8:48137/1‑149 (MQ=255)
gCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGa                                                                                                                                                  <  5:296961/149‑1 (MQ=255)
  cATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATg                                                                                                                                                >  3:209529/1‑149 (MQ=255)
   aatGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATgg                                                                                                                                               >  8:181922/3‑149 (MQ=255)
    ttGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGc                                                                                                                                              <  1:263011/149‑1 (MQ=255)
       cGTACCCCCCGGCTGCCGCGCCGCCCCTGGCAGGGGGGGTCCTTGCCGCCTGGTTCGCTGCTCACCGTCATGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  7:232269/149‑1 (MQ=255)
       cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGTGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  1:253503/149‑1 (MQ=255)
       cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  7:181922/149‑1 (MQ=255)
       cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGaa                                                                                                                                           <  3:208039/149‑1 (MQ=255)
        gTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAg                                                                                                                                          >  8:240790/1‑149 (MQ=255)
         tACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGt                                                                                                                                         >  7:4752/1‑149 (MQ=255)
              ccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcg                                                                                                                                    >  4:333585/1‑149 (MQ=255)
                 ggCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGGGAAGTTTGCGCTg                                                                                                                                 >  7:267491/1‑149 (MQ=255)
                 ggCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTg                                                                                                                                 <  3:343020/149‑1 (MQ=255)
                    tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCg                                                                                                                              >  6:258137/1‑149 (MQ=255)
                    tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCg                                                                                                                              >  6:242886/1‑149 (MQ=255)
                    tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCg                                                                                                                              >  5:309518/1‑149 (MQ=255)
                      ccGCGCAGCTCCTGTAAGGGGGTGTCCTTGCGGCGTGGTTCGCTTCTCAGCGTCTTGCGATTGAGCGTGTCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCAGGCCATCGTTGATGGCGAAGTTTGCGCTGCCGtt                                                                                                                            <  1:262212/149‑1 (MQ=255)
                       cgcgCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCGGCGTCTTGCGATTGAGCGTGGCGATGGATTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCAACGTTGATGGCGAAGTTTGCGGTGCCGtt                                                                                                                            >  5:108821/1‑148 (MQ=255)
                         cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                         <  6:309518/149‑1 (MQ=255)
                         cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                         <  1:373316/149‑1 (MQ=255)
                         cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAg                                                                                                                         <  3:250930/149‑1 (MQ=255)
                          gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        <  8:4752/149‑1 (MQ=255)
                          gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        <  7:240790/149‑1 (MQ=255)
                               tCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcgc                                                                                                                                   >  3:254697/1‑133 (MQ=255)
                               tCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTgcgc                                                                                                                                   <  4:254697/133‑1 (MQ=255)
                                 cTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGa                                                                                                                 >  6:320438/1‑149 (MQ=255)
                                  tGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                >  7:344873/1‑149 (MQ=255)
                                   ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                >  6:24788/1‑148 (MQ=255)
                                    gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                <  6:264138/147‑1 (MQ=255)
                                    gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGAc                                                                                                                <  8:344873/147‑1 (MQ=255)
                                      aGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                            <  5:258137/149‑1 (MQ=255)
                                        ggggTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTg                                                                                                           <  5:252897/148‑1 (MQ=255)
                                              tCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTg                                                                                                    >  7:126545/1‑149 (MQ=255)
                                                cTTGCCGTGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGc                                                                                                  >  4:85906/1‑149 (MQ=255)
                                                 ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCt                                                                                                                                  <  3:60730/116‑1 (MQ=255)
                                                 ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCt                                                                                                                                  >  4:60730/1‑116 (MQ=255)
                                                 ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGc                                                                                                                                >  1:243836/1‑118 (MQ=255)
                                                 ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGc                                                                                                                                <  2:243836/118‑1 (MQ=255)
                                                 ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                 >  5:180602/1‑149 (MQ=255)
                                                   gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                                 <  6:103368/147‑1 (MQ=255)
                                                      gcgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        >  2:315764/1‑121 (MQ=255)
                                                       cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        <  4:45065/120‑1 (MQ=255)
                                                       cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAgg                                                                                                                        >  3:45065/1‑120 (MQ=255)
                                                           gTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGTCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTg                                                                                       <  5:42953/149‑1 (MQ=255)
                                                            ttCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTgg                                                                                      >  6:62976/1‑149 (MQ=255)
                                                               gctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGTCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCt                                                                                    <  4:135318/148‑1 (MQ=255)
                                                                ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  <  5:242886/149‑1 (MQ=255)
                                                                ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  >  7:165169/1‑149 (MQ=255)
                                                                ctgctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  <  8:126545/149‑1 (MQ=255)
                                                                   ctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  >  1:241117/1‑146 (MQ=255)
                                                                    tCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                                  >  1:75996/1‑145 (MQ=255)
                                                                      aGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGAt                                                                             <  1:317902/148‑1 (MQ=255)
                                                                       gCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCa                                                                           <  8:165169/149‑1 (MQ=255)
                                                                          tCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTg                                                                        <  4:294227/149‑1 (MQ=255)
                                                                           cTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTgg                                                                       >  1:164339/1‑149 (MQ=255)
                                                                              gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGCCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTgg                                                                       <  4:209529/146‑1 (MQ=255)
                                                                              gCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtc                                                                     <  7:252817/148‑1 (MQ=255)
                                                                               cGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtca                                                                    <  5:81379/148‑1 (MQ=255)
                                                                                gATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                                  <  1:13786/149‑1 (MQ=255)
                                                                                     aGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggtcag                                                             <  3:333585/149‑1 (MQ=255)
                                                                                       cGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagca                                                           <  2:59788/149‑1 (MQ=255)
                                                                                        gTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTCACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagcag                                                          <  5:320438/149‑1 (MQ=255)
                                                                                         tGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCt                                                                                    <  5:332139/122‑1 (MQ=255)
                                                                                         tGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCt                                                                                    >  6:332139/1‑122 (MQ=255)
                                                                                                       cacaACTGCTGTTCCAGGTTGGCATAGACCAGCACGCCATCGTTGATGGCGAAGATTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcggcgg                                           >  4:255415/1‑149 (MQ=255)
                                                                                                          aaCTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGgcggcg                                         >  8:116426/1‑148 (MQ=255)
                                                                                                             tgctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTg                                     >  6:2095/1‑149 (MQ=255)
                                                                                                               ctgTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTg                                     <  2:75996/147‑1 (MQ=255)
                                                                                                                  ttCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCt                                >  2:4262/1‑149 (MQ=255)
                                                                                                                    ccAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCtt                               <  7:48137/148‑1 (MQ=255)
                                                                                                                        gTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCtt                          >  3:179128/1‑149 (MQ=255)
                                                                                                                           ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                        <  2:241117/148‑1 (MQ=255)
                                                                                                                           ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                        <  6:180602/148‑1 (MQ=255)
                                                                                                                           ggCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCg                        <  5:62976/148‑1 (MQ=255)
                                                                                                                              aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTc                     >  4:276984/1‑148 (MQ=255)
                                                                                                                              aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTc                     <  3:276984/148‑1 (MQ=255)
                                                                                                                              aTTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTACAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCt                    >  2:284836/1‑149 (MQ=255)
                                                                                                                                ctacgaGCACGCCATCGTTGCTGGCGAAGTTTGCGCTGCCGTTCGGGGTGGCGACCAGTGCTTTGTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCtt                   <  8:63066/143‑1 (MQ=255)
                                                                                                                                  aCCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                >  6:108406/1‑149 (MQ=255)
                                                                                                                                   ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCGGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCCGCGGCGTCTGTTCCTTGCCTTCGGTCtt                   >  7:63066/1‑145 (MQ=255)
                                                                                                                                   ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGa                                                  >  5:87216/1‑114 (MQ=255)
                                                                                                                                   ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGa                                                  <  6:87216/114‑1 (MQ=255)
                                                                                                                                   ccAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat                >  1:226701/1‑148 (MQ=255)
                                                                                                                                   aaacttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGGCGtgat                >  4:263116/7‑148 (MQ=255)
                                                                                                                                    aacttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTtgatg               >  6:128337/6‑148 (MQ=255)
                                                                                                                                     aGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa             <  7:116426/149‑1 (MQ=255)
                                                                                                                                     aGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa             <  5:108406/149‑1 (MQ=255)
                                                                                                                                      cttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTCCCTTGCCTTCGGTCTTGATGAAg            >  3:273920/4‑149 (MQ=255)
                                                                                                                                       ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggt                                                                <  6:142767/94‑1 (MQ=255)
                                                                                                                                       ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagggt                                                                >  5:142767/3‑96 (MQ=255)
                                                                                                                                       ttcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCTCTTTACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgt           >  6:6598/3‑149 (MQ=255)
                                                                                                                                        tcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTg                                                                        >  3:283479/2‑87 (MQ=255)
                                                                                                                                        tcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTg                                                                        <  4:283479/86‑1 (MQ=255)
                                                                                                                                        tcGCCATCAACGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGACCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTATTGATGAAgtg          >  6:110987/2‑149 (MQ=255)
                                                                                                                                         cGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgtgt         <  5:24788/149‑1 (MQ=255)
                                                                                                                                          gCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt        >  2:106373/1‑149 (MQ=255)
                                                                                                                                           ccATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt        >  7:97282/1‑148 (MQ=255)
                                                                                                                                                gTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCg  >  8:306521/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGCCATCGTTGATGGCGAAGTTTGCGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGGCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTTCCACCG  >  NC_002947/345974‑346266

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: