Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,310,451 T→A 5.9% D363V (GAC→GTC)  PP_1144 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,310,4510TA 5.9% 324.5 / 10.9 119D363V (GAC→GTC) PP_1144membrane protein
Reads supporting (aligned to +/- strand):  ref base T (56/55);  new base A (5/2);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

Only 100 of 121 total aligned reads displayed.

GTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCG  >  NC_002947/1310307‑1310600
                                                                                                                                                |                                                                                                                                                     
gtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTgg                                                                                                                                                  <  1:514108/150‑1 (MQ=255)
  ggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGtt                                                                                                                                                <  1:1561854/150‑1 (MQ=255)
    tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGc                                                                                                                                              >  1:1808849/1‑150 (MQ=255)
    tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGc                                                                                                                                              >  2:19634/1‑150 (MQ=255)
    tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGc                                                                                                                                              >  1:1945280/1‑150 (MQ=255)
      aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                            >  2:356164/1‑150 (MQ=255)
      aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                            >  1:851224/1‑150 (MQ=255)
      aTGACAACGCGCGGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg                                                                                                                                            <  2:1358341/150‑1 (MQ=255)
         aCCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTa                                                                                                                                         <  1:246249/150‑1 (MQ=255)
            aTGCGCAGCCCGCACAGGTCGGTGGGGTGTTTAAGAAGTTGCATGATGCGGTGGGACAGGAGGCGGGCTTCGTCGCGGCTGTTTAGGCTGGTGACCACACAGAACTCGTCGCCACCACAGCGCGCAACTAGGGCGGGGTTGCGGGTGGCg                                                                                                                                      >  2:1028266/1‑150 (MQ=255)
                gccgccCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCt                                                                                                                                  >  1:896029/1‑150 (MQ=255)
                 ccgcccgcCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCtt                                                                                                                                 >  2:1950528/1‑150 (MQ=255)
                        cccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTgg                                                                                                                          <  2:1268123/150‑1 (MQ=255)
                        cccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGATCAGGCTGGTGACCACACAGGACTCGTCGCCCCCAAAGCGCGCAACCCGGTCGGGGTTGCGGGTGGCGGCCTTGGTGTgg                                                                                                                          >  2:2148571/1‑150 (MQ=255)
                              tCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCg                                                                                                                    <  2:2044312/150‑1 (MQ=255)
                                 aTGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATc                                                                                                                 >  1:94665/1‑150 (MQ=255)
                                 aTGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATc                                                                                                                 <  1:1930475/150‑1 (MQ=255)
                                  tGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGGTGGCGaaga                                                                                                                >  2:961066/1‑147 (MQ=255)
                                    gggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAcc                                                                                                              <  1:19634/150‑1 (MQ=255)
                                     ggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAccg                                                                                                             >  1:2430156/1‑149 (MQ=255)
                                      gTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCtt                                                                                                            <  1:127497/150‑1 (MQ=255)
                                        cTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGa                                                                                                          <  1:321301/150‑1 (MQ=255)
                                         tttCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGag                                                                                                         >  2:2296114/1‑150 (MQ=255)
                                         tttCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGag                                                                                                         >  1:1473644/1‑150 (MQ=255)
                                             ccccGTTCCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                     <  1:1028266/148‑1 (MQ=255)
                                             aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                     <  2:278610/150‑1 (MQ=255)
                                             aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGAGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                     <  2:2421971/150‑1 (MQ=255)
                                              tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgc                                                                                                                                            <  2:351781/110‑2 (MQ=255)
                                              tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgc                                                                                                                                            >  1:351781/1‑109 (MQ=255)
                                               ccGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTc                                                                                                   <  2:1640422/150‑1 (MQ=255)
                                                 gTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTCGCTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGt                                                                                                 <  1:61820/150‑1 (MQ=255)
                                                    gCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCAc                                                                                              >  1:1134072/1‑150 (MQ=255)
                                                     cATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCAcc                                                                                             <  2:2175000/150‑1 (MQ=255)
                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                            >  2:2120167/1‑150 (MQ=255)
                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                            >  2:2217508/1‑150 (MQ=255)
                                                      atgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCg                                                                                            <  2:1945280/150‑1 (MQ=255)
                                                       tgatgCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCgg                                                                                           >  1:2176236/1‑150 (MQ=255)
                                                         atgCGCTGGGCCAGGTGGCGGGCTTGGTCTCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGcg                                                                                         >  1:1638665/1‑150 (MQ=255)
                                                         atgCGCTGGGCCAGGTGGCGGGCTTATTCTGGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGcg                                                                                         >  1:1638598/1‑150 (MQ=255)
                                                              cTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTg                                                                                    <  2:851224/150‑1 (MQ=255)
                                                                   ccAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGc                                                                                                                                    <  2:1265881/97‑1 (MQ=255)
                                                                   ccAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCaa                                                                               >  1:2003792/1‑150 (MQ=255)
                                                                     aGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAgg                                                                             <  2:1134072/150‑1 (MQ=255)
                                                                             ggCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGtt                                                                     <  2:1554057/150‑1 (MQ=255)
                                                                                    tcgNGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGc                                                              >  2:1504927/1‑150 (MQ=255)
                                                                                    tcgCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCt                                                                                                             >  1:2084891/1‑103 (MQ=255)
                                                                                    tcgCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCt                                                                                                             <  2:2084891/103‑1 (MQ=255)
                                                                                    tcgCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGc                                                              <  1:138884/150‑1 (MQ=255)
                                                                                      gcgGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCtt                                                            <  1:2207196/150‑1 (MQ=255)
                                                                                         gCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaa                                                         <  2:94665/150‑1 (MQ=255)
                                                                                         gCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTAGTTGCGGGTAGCCGCCTTGATGTAGTTGGCGATCTCCTTGAGCTGCTCGGCAGCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaa                                                         <  2:605378/150‑1 (MQ=255)
                                                                                            gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGAGCCTCAGGCTTGCGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATg                                                      <  2:1638598/150‑1 (MQ=255)
                                                                                             ttCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATgg                                                     >  1:176049/1‑150 (MQ=255)
                                                                                                aGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCg                                                  >  1:2451963/1‑150 (MQ=255)
                                                                                                   cTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATg                                               >  2:297689/1‑150 (MQ=255)
                                                                                                   cTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATg                                               <  2:1060184/150‑1 (MQ=255)
                                                                                                       tGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGa                                           <  1:297689/150‑1 (MQ=255)
                                                                                                        gACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAg                                          <  1:2346491/150‑1 (MQ=255)
                                                                                                          ccACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGa                                        <  2:425232/150‑1 (MQ=255)
                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGa                                                                   >  2:463552/1‑122 (MQ=255)
                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGa                                                                   <  1:463552/122‑1 (MQ=255)
                                                                                                           cacacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGaa                                       >  2:1962356/1‑150 (MQ=255)
                                                                                                            acacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGGTACGCTTGAAATGGTCGATGTCGAGGAAc                                      >  1:389305/1‑150 (MQ=255)
                                                                                                            acacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGCGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGGCGTTGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                      <  1:1022887/150‑1 (MQ=255)
                                                                                                              acaGAACTCGTCGCCCCCAAGGCGCGCGACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACAt                                    <  2:389305/150‑1 (MQ=255)
                                                                                                              acaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACAt                                    >  1:1943091/1‑150 (MQ=255)
                                                                                                              acaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACAt                                    >  1:38221/1‑150 (MQ=255)
                                                                                                                aGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGa                                  >  2:892682/1‑150 (MQ=255)
                                                                                                                      cgtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCa                            <  2:2003792/150‑1 (MQ=255)
                                                                                                                          gCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                     <  1:1016481/73‑1 (MQ=255)
                                                                                                                          gCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagcagc                                                                                                     >  2:1016481/1‑73 (MQ=255)
                                                                                                                           ccaccaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCgg                       >  1:1533646/1‑150 (MQ=255)
                                                                                                                            caccaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAAc                                      <  2:1684255/134‑1 (MQ=255)
                                                                                                                             accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAcc                                                                                                              >  2:1496959/1‑61 (MQ=255)
                                                                                                                             accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAcc                                                                                                              <  1:1496959/61‑1 (MQ=255)
                                                                                                                             accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt                     >  2:183496/1‑150 (MQ=255)
                                                                                                                               caAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCg                                            <  1:1376287/125‑1 (MQ=255)
                                                                                                                               caAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCg                                            >  2:1376287/1‑125 (MQ=255)
                                                                                                                               caAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtttt                   <  1:183496/150‑1 (MQ=255)
                                                                                                                                aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc                  <  1:588550/150‑1 (MQ=255)
                                                                                                                                aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc                  >  1:1829872/1‑150 (MQ=255)
                                                                                                                                  aGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc                <  2:1943091/150‑1 (MQ=255)
                                                                                                                                    cccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTg              >  2:2338077/3‑150 (MQ=255)
                                                                                                                                     ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTgg             >  2:204019/2‑150 (MQ=255)
                                                                                                                                      cgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                    >  2:964173/1‑110 (MQ=255)
                                                                                                                                      cgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGt                                                    <  1:964173/110‑1 (MQ=255)
                                                                                                                                       gcAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCt                                                             <  1:1273544/100‑1 (MQ=255)
                                                                                                                                       gcAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCt                                                             >  2:1273544/1‑100 (MQ=255)
                                                                                                                                        cAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCt          >  2:185767/1‑150 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGAt                                                                  >  2:119887/1‑93 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGAt                                                                  <  1:119887/93‑1 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAg                                          <  1:930957/117‑1 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAg                                          >  2:930957/1‑117 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTg         >  1:1801984/1‑150 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTg         >  2:1815275/1‑150 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGAAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGc                                                              <  2:257292/97‑1 (MQ=255)
                                                                                                                                         aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGAAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGc                                                              >  1:257292/1‑97 (MQ=255)
                                                                                                                                           ccAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGAt       <  1:2148571/150‑1 (MQ=255)
                                                                                                                                                tCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCg  <  2:38221/150‑1 (MQ=255)
                                                                                                                                                tCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCg  >  2:51741/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCG  >  NC_002947/1310307‑1310600

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: