Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,727,168 T→A 6.3% L113Q (CTG→CAG)  PP_3292 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,727,1680TA 6.3% 200.2 / 3.0 65L113Q (CTG→CAG) PP_3292membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGC  >  NC_002947/3727019‑3727317
                                                                                                                                                      |                                                                                                                                                     
aGGCTGGGCGGTGCCGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCct                                                                                                                                                       <  1:1156435/150‑1 (MQ=255)
    tGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCctgctg                                                                                                                                                   <  2:1011858/150‑1 (MQ=255)
     gggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTgg                                                                                                                                                  <  1:1288787/150‑1 (MQ=255)
      ggCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTggg                                                                                                                                                 <  1:1190511/150‑1 (MQ=255)
         ggTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGcac                                                                                                                                              >  2:879914/1‑150 (MQ=255)
          ttGCAGCCGGGGGTGAAGCGCAGCGTTTGCAAAGCGAATTTCACCGGCCAGGAGACAGGCGCGCATGCTAGACTGGGGTTCCTTCGCTCAGGGCCTCCGGAATCGGGCAGCCTGCCGGCAGTGCTGTTGTGCCTG‑GCCCTGCGGGGcaca                                                                                                                                             <  2:669592/149‑1 (MQ=255)
           tGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAg                                                                                                                                            <  2:280634/150‑1 (MQ=255)
                   ggCGTGAAGCGCACCGTTGGTAAAGCCATTTTCAGCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAcaaca                                                                                                                                    <  2:263058/150‑5 (MQ=255)
                           gcgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggctggc                                                                                                                            >  2:1603808/1‑150 (MQ=255)
                            cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                           <  1:534579/150‑1 (MQ=255)
                              cACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTAGGGAATCCAACATCCTCCTGGTGGTTGTGTTGGTTCTG‑GGTCTGCTGGGCACAGGCACCTTTTGGCTGGCAgg                                                                                                                         <  2:1006705/150‑1 (MQ=255)
                               aCCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAggg                                                                                                                        <  1:764549/150‑1 (MQ=255)
                                  gTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGgcgc                                                                                                                     >  1:688133/1‑150 (MQ=255)
                                        aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCtttt                                                                                                                                   >  1:727322/1‑130 (MQ=255)
                                        aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCtttt                                                                                                                                   <  2:727322/130‑1 (MQ=255)
                                        aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTg                                                                                                               >  1:1614152/1‑150 (MQ=255)
                                         aaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGa                                                                                                              <  2:688133/150‑1 (MQ=255)
                                                tttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGAt                                                                                                       >  2:261690/1‑150 (MQ=255)
                                                    aaCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTg                                                                                                   <  2:1599530/150‑1 (MQ=255)
                                                     aCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTgg                                                                                                  <  1:261690/150‑1 (MQ=255)
                                                             gagaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGc                                                                                                 <  1:186569/143‑1 (MQ=255)
                                                             gagaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGc                                                                                                 >  2:186569/1‑143 (MQ=255)
                                                             gagaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGc                                                                                          <  1:1595804/150‑1 (MQ=255)
                                                                aCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAg                                                                                       <  2:299867/150‑1 (MQ=255)
                                                                aCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAg                                                                                       <  2:1639118/150‑1 (MQ=255)
                                                                   ggCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggc                                                                                                                                <  2:181927/106‑1 (MQ=255)
                                                                   ggCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggc                                                                                                                                >  1:181927/1‑106 (MQ=255)
                                                                    gCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGTAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCa                                                                                   >  2:886867/1‑150 (MQ=255)
                                                                    gCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCa                                                                                   >  1:395776/1‑150 (MQ=255)
                                                                     cgatctggCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCaa                                                                                  <  1:130619/144‑1 (MQ=255)
                                                                         cAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTg                                                                              <  1:1035188/150‑1 (MQ=255)
                                                                            gCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGGGGCGGAGATCGTGGCGAACTGCGAGGCCAAcagatcg                                                                           >  2:130619/1‑144 (MQ=255)
                                                                                    ggCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGt                                                                   >  1:489109/1‑150 (MQ=255)
                                                                                        ttCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGat                                                               >  1:649189/1‑150 (MQ=255)
                                                                                           cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTCTAGTTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                            >  1:298393/1‑150 (MQ=255)
                                                                                            ttCGCTCAGTGTCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGTGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAc                                                           <  2:653666/150‑1 (MQ=255)
                                                                                            ttCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggctg                                                                                                                              >  1:32985/1‑83 (MQ=255)
                                                                                            ttCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggctg                                                                                                                              <  2:32985/83‑1 (MQ=255)
                                                                                            ttCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAc                                                           >  1:189371/1‑150 (MQ=255)
                                                                                            ttCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAc                                                           >  2:654079/1‑150 (MQ=255)
                                                                                               gCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCctgct                                                                                                                                                    <  2:1275252/58‑1 (MQ=255)
                                                                                               gCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCctgct                                                                                                                                                    >  1:1275252/1‑58 (MQ=255)
                                                                                                 tCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                                      <  2:988867/150‑1 (MQ=255)
                                                                                                 tCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                                      >  1:577679/1‑150 (MQ=255)
                                                                                                       gCTAGGGAATCGAACATGCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaa                                                >  2:1594702/1‑150 (MQ=255)
                                                                                                           gggAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAAc                                            >  1:300339/1‑150 (MQ=255)
                                                                                                                  gAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                                     >  2:939591/1‑150 (MQ=255)
                                                                                                                  gAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCCGCTCGGCACAGGCACATTTTGGCTGACAGGGCGGCGAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                                     >  2:800487/1‑150 (MQ=255)
                                                                                                                      aTCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCtt                                 >  1:1116125/1‑150 (MQ=255)
                                                                                                                          tgctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCa                                         <  1:417041/138‑1 (MQ=255)
                                                                                                                          tgctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCa                                         >  2:417041/1‑138 (MQ=255)
                                                                                                                           gctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCgg                            <  2:649189/150‑1 (MQ=255)
                                                                                                                            ctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTg                                                                                                   <  2:1372498/78‑1 (MQ=255)
                                                                                                                            ctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTg                                                                                                   >  1:1372498/1‑78 (MQ=255)
                                                                                                                                tggtgCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAAc                       <  2:395776/150‑1 (MQ=255)
                                                                                                                                   ccaaaaggtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCgg                    >  1:582669/8‑150 (MQ=255)
                                                                                                                                    gCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgt                   <  2:189371/150‑1 (MQ=255)
                                                                                                                                      tgttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGTAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt                 <  1:886867/150‑1 (MQ=255)
                                                                                                                                        ttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGAt                                                   <  2:1272353/114‑1 (MQ=255)
                                                                                                                                        ttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGAt                                                   >  1:1272353/1‑114 (MQ=255)
                                                                                                                                          gtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAAcc                                                                                <  2:996004/83‑1 (MQ=37)
                                                                                                                                          gtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAAcc                                                                                >  1:996004/1‑83 (MQ=37)
                                                                                                                                          gtgCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCAc             <  2:300339/150‑1 (MQ=255)
                                                                                                                                              cTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcg                                      <  2:1318939/121‑1 (MQ=255)
                                                                                                                                              cTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcg                                      >  1:1318939/1‑121 (MQ=255)
                                                                                                                                              cTG‑GCCCTGCTGGGCACAGGAACCTGTTGGCGGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGc         <  2:1614152/150‑1 (MQ=255)
                                                                                                                                               tG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGcc        >  1:898343/1‑150 (MQ=255)
                                                                                                                                                g‑gCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCt       <  2:489109/150‑1 (MQ=255)
                                                                                                                                                    ccTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgc    >  1:473306/1‑150 (MQ=255)
                                                                                                                                                      tgctgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt                 <  2:1595108/135‑1 (MQ=255)
                                                                                                                                                      tgctgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt                 >  1:1595108/1‑135 (MQ=255)
                                                                                                                                                      tgctgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCATGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCCCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgcgc  >  2:1064065/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGC  >  NC_002947/3727019‑3727317

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