Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,872 C→A 5.5% Q94H (CAG→CAT)  PP_5039 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,8720CA 5.5% 235.8 / 3.0 72Q94H (CAG→CAT) PP_5039hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CGTCCAGCAGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGT  >  NC_002947/5741727‑5742020
                                                                                                                                                 |                                                                                                                                                     
cGTCCAGCAGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTg                                                                                                                                                     <  7:79334/148‑1 (MQ=255)
  tCCAGCAGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTggg                                                                                                                                                   >  8:10910/1‑148 (MQ=255)
    cagcagCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCAc                                                                                                                                                >  6:98741/1‑149 (MQ=255)
      gcagcaTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              <  1:183271/149‑1 (MQ=255)
      gcagcaTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              <  8:174331/149‑1 (MQ=255)
      gcagcaTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              <  6:281998/149‑1 (MQ=255)
      gcagcaTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              <  3:63384/149‑1 (MQ=255)
      gcagcaTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAACAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              >  3:255286/1‑149 (MQ=255)
          caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGt                                                                                                                                          <  1:69070/149‑1 (MQ=255)
               gATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCggagatag                                                                                                                                     >  5:282556/1‑143 (MQ=255)
               gATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCg                                                                                                                                            <  6:282556/142‑1 (MQ=255)
                aTAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCt                                                                                                                                    <  7:149184/149‑1 (MQ=255)
                   aCGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGa                                                                                                                                 <  3:181222/149‑1 (MQ=255)
                    cGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAACAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGaa                                                                                                                                <  4:255286/149‑1 (MQ=255)
                          aCGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              <  8:72017/129‑1 (MQ=255)
                          aCGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGc                                                                                                                                              >  7:72017/1‑129 (MQ=255)
                            gCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGcgcgcg                                                                                                                        <  5:136702/149‑1 (MQ=255)
                            gCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGcgcgcg                                                                                                                        >  7:159564/1‑149 (MQ=255)
                              acaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcggc                                                                                                                      >  7:300147/1‑149 (MQ=255)
                                  gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCgtgt                                                                                                                  >  1:75052/1‑149 (MQ=255)
                                        ccAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                             >  1:83551/1‑148 (MQ=255)
                                         cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCgg                                                                                                             >  5:92471/1‑147 (MQ=255)
                                         cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTa                                                                                                           <  8:159564/149‑1 (MQ=255)
                                                gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATg                                                                                                    >  2:28859/1‑149 (MQ=255)
                                                gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACACCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATg                                                                                                    >  3:316630/1‑149 (MQ=255)
                                                gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACATGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATg                                                                                                    >  4:231585/1‑149 (MQ=255)
                                                gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACATGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATg                                                                                                    >  4:227294/1‑149 (MQ=255)
                                                 cACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                               >  7:305520/1‑121 (MQ=255)
                                                 cACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                               <  8:305520/121‑1 (MQ=255)
                                                 cACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   >  1:303948/1‑149 (MQ=255)
                                                 cACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   >  1:308679/1‑149 (MQ=255)
                                                    ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   <  6:183623/146‑1 (MQ=255)
                                                    ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   >  5:183623/1‑146 (MQ=255)
                                                       ggCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCTGTAGATCATGc                                                                                                   <  5:191591/143‑1 (MQ=255)
                                                       ggCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTCCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   >  6:191591/1‑143 (MQ=255)
                                                            gCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   >  3:220569/1‑138 (MQ=255)
                                                            gCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGc                                                                                                   <  4:220569/138‑1 (MQ=255)
                                                            gCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTgg                                                                                        >  8:19063/1‑149 (MQ=255)
                                                             cTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTgg                                                                                        <  5:98741/148‑1 (MQ=255)
                                                                cAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTc                                                                                      >  7:121758/1‑147 (MQ=255)
                                                                  gTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAg                                                                                   >  4:172505/1‑148 (MQ=255)
                                                                      cagcagCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaa                                                                              >  1:233696/1‑149 (MQ=255)
                                                                      cagcagCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaa                                                                              <  2:238164/149‑1 (MQ=255)
                                                                      cagcagCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaa                                                                              >  1:238164/1‑149 (MQ=255)
                                                                      cagcagCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaa                                                                              <  2:233696/149‑1 (MQ=255)
                                                                           gcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACATGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATc‑agc                                                                        <  3:227294/149‑1 (MQ=255)
                                                                           gcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACATGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATc‑agc                                                                        <  3:231585/149‑1 (MQ=255)
                                                                              cGAACCCGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAgg                                                                     <  8:264794/149‑1 (MQ=255)
                                                                               gAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCAt                                                                                                     <  7:15881/117‑1 (MQ=255)
                                                                               gAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGCCCCTGCATGGCGACCACCTGGGCACGCACCTCGCCCTCGAACGCGCGCGGCGTGTCGCGGTAGATCAt                                                                                                     >  8:15881/1‑117 (MQ=255)
                                                                                     gcgcTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCCGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCccagc                                                               >  3:67018/1‑148 (MQ=255)
                                                                                      cgcTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCGTGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCccagcca                                                             <  1:254606/149‑1 (MQ=255)
                                                                                           gATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCagccagcc                                                          <  6:92471/147‑1 (MQ=255)
                                                                                                gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCggg                                                   <  2:303948/149‑1 (MQ=255)
                                                                                                gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCggg                                                   <  7:10910/149‑1 (MQ=255)
                                                                                                gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCggg                                                   >  4:176784/1‑149 (MQ=255)
                                                                                                gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCGCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCggg                                                   <  2:308679/149‑1 (MQ=255)
                                                                                                  agcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGt                                                 <  3:172505/149‑1 (MQ=255)
                                                                                                          gTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                            <  1:28859/147‑1 (MQ=255)
                                                                                                             gTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCAc                                       <  2:75052/148‑1 (MQ=255)
                                                                                                                  aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCg                                                     <  8:22867/129‑1 (MQ=255)
                                                                                                                  aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCg                                                     >  7:22867/1‑129 (MQ=255)
                                                                                                                     ccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCg                                 >  8:225571/1‑146 (MQ=255)
                                                                                                                     ccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCg                                 <  7:225571/146‑1 (MQ=255)
                                                                                                                       ttCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGccc                            >  3:235002/1‑149 (MQ=255)
                                                                                                                       ttCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTACCACTCTTCGAGcc                             >  2:179801/1‑148 (MQ=255)
                                                                                                                         cTGCAGGCCCTGCATGGCGACCACCTGGGCACACCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGTTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTTCCACTCTTCGAGccc                            <  4:316630/147‑1 (MQ=255)
                                                                                                                            cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCg                       >  6:306906/1‑149 (MQ=255)
                                                                                                                            cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCg                       >  4:109482/1‑149 (MQ=255)
                                                                                                                                 ccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAgg                  >  3:294741/1‑149 (MQ=255)
                                                                                                                                       tGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc            >  3:314137/1‑149 (MQ=255)
                                                                                                                                       tGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc            >  1:59134/1‑149 (MQ=255)
                                                                                                                                             ccaccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCgcc       >  6:316428/1‑146 (MQ=255)
                                                                                                                                               accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTTGCAGCTGGTCGAGCAAAATC‑AGCAGGCccagcca                                                             <  2:134113/92‑1 (MQ=255)
                                                                                                                                               accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCccagcca                                                             >  1:134113/1‑92 (MQ=255)
                                                                                                                                               accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAg    <  4:235002/149‑1 (MQ=255)
                                                                                                                                                 cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGt  >  3:242213/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CGTCCAGCAGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCGCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATC‑AGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGT  >  NC_002947/5741727‑5742020

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: ‑