Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,946 C→T 6.5% E88E (GAG→GAA)  PP_5039 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,9460CT 6.5% 226.5 / 4.7 77E88E (GAG→GAA) PP_5039hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base T (2/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

GATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGG  >  NC_002947/5741798‑5742094
                                                                                                                                                    |                                                                                                                                                     
gATAACGCTCCACGCACAGGTTCTGCGGCGCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCt                                                                                                                                                      >  1:50324/1‑150 (MQ=255)
gATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCCCGCCGGTCGCCCt                                                                                                                                                      >  1:1908941/1‑150 (MQ=255)
     cGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 >  1:710333/1‑150 (MQ=255)
        tCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                              <  2:1063186/150‑1 (MQ=255)
          cACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcg                                                                                                                                            >  1:1166953/1‑150 (MQ=255)
          cACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcg                                                                                                                                            >  1:193841/1‑150 (MQ=255)
           aCGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcgcgc                                                                                                                                           <  1:1606690/150‑1 (MQ=255)
             gCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCgcgg                                                                                                                                         >  2:1666962/1‑150 (MQ=255)
              cacaGGTTCTGTCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                                        <  2:710307/150‑7 (MQ=255)
              cacaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGTATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGgcgaccggc                                                                                                                                        <  2:163840/150‑7 (MQ=255)
               acaGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGgcgaccggcg                                                                                                                                       >  1:836278/1‑143 (MQ=255)
                   gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 <  2:892605/136‑1 (MQ=255)
                   gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                                 >  1:892605/1‑136 (MQ=255)
                   gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTACTCGTATGCCGTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTc                                                                                                                                   >  1:1772373/1‑150 (MQ=255)
                      cTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTcg‑cg                                                                                                                               >  2:1519927/1‑150 (MQ=255)
                          cAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAg                                                                                                                           >  1:1200002/1‑150 (MQ=255)
                             tCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATc                                                                                                                        >  1:542535/1‑150 (MQ=255)
                             cctaGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATc                                                                                                                        >  2:505077/4‑150 (MQ=255)
                               cAGCACTACGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCCTCTGCAGGCCCTGCATGGCGACCACCTGTGCACGCCGGTCCCCCTCGAAGGCGCGCGGCGTGTCG‑CGGCAGATccc                                                                                                                      >  2:90827/1‑148 (MQ=255)
                                atcattgcGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATg                                                                                                                     <  2:367075/143‑1 (MQ=255)
                                 gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGc                                                                                                                    <  2:193841/150‑1 (MQ=255)
                                 gCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGc                                                                                                                    <  2:542535/150‑1 (MQ=255)
                                     ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                              <  1:1615552/121‑1 (MQ=255)
                                     ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                              >  2:1615552/1‑121 (MQ=255)
                                       cGGCATGCTCCAGTACCAGCAGCACGCCCACGCGCTGGCACGGCAGCAGCTCGTAGTCCCAGCCCTTCTGCCGCTGCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCa                                                                                                              <  1:1067265/150‑1 (MQ=255)
                                         gCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGc                                                                                                            >  1:1811049/1‑150 (MQ=255)
                                           attcTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTg                                                                                                          >  1:274651/4‑150 (MQ=255)
                                           aTGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTg                                                                                                          <  2:1488370/150‑1 (MQ=255)
                                            tGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTgg                                                                                                         <  2:1440719/150‑1 (MQ=255)
                                              cTCCAGTACCAGCAGCACGAACACGCACAGGATCGTCAGCAGCTAATAGTCCCAGCCCTGCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGCAGCCGCCCGGCGTGTCG‑CGCTAGATCATGCCTGGCAGCTGGTc                                                                                                       <  2:1552445/150‑1 (MQ=255)
                                                ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                     <  1:1844185/150‑1 (MQ=255)
                                                ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                     <  1:1844158/150‑1 (MQ=255)
                                                ccAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATCGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                     <  1:233370/150‑1 (MQ=255)
                                                       cagcagCACGAACACGCGCTGGATCGGCAGCAGCTAGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGGAAGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCaaaaa                                                                                              >  2:61714/1‑149 (MQ=255)
                                                          cagcaCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATcag                                                                                           <  1:395913/150‑1 (MQ=255)
                                                              aCGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAgg                                                                                       <  1:1197836/150‑1 (MQ=255)
                                                                    acGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCcaagc                                                                                 <  1:902291/150‑1 (MQ=255)
                                                                            gATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTg                                                                         >  1:166024/1‑150 (MQ=255)
                                                                                ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                     <  1:1519927/150‑1 (MQ=255)
                                                                                ggCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCggg                                                                     <  1:248282/150‑1 (MQ=255)
                                                                                 gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                    >  2:522324/1‑150 (MQ=255)
                                                                                 gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                                    >  2:1280966/1‑150 (MQ=255)
                                                                                    gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                                 <  1:68905/150‑1 (MQ=255)
                                                                                    gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                                 <  2:274651/150‑1 (MQ=255)
                                                                                    gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                                 <  2:1166953/150‑1 (MQ=255)
                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAGGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                              >  1:1100788/1‑150 (MQ=255)
                                                                                       gcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgat                                                              <  2:935835/150‑3 (MQ=255)
                                                                                        cTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgctg                                                             >  2:1001662/1‑150 (MQ=255)
                                                                                               tCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGACCGGCgtg                                                                                                                                     >  2:1519865/1‑72 (MQ=18)
                                                                                               tCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGACCGGCgtg                                                                                                                                     <  1:1519865/72‑1 (MQ=18)
                                                                                                 ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                                    >  1:953992/1‑150 (MQ=255)
                                                                                                 ccAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTc                                                    >  1:672591/1‑150 (MQ=255)
                                                                                                   aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                     <  2:1510477/99‑1 (MQ=255)
                                                                                                   aGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGa                                                                                                     >  1:1510477/1‑99 (MQ=255)
                                                                                                     cccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGc                                                <  2:672591/150‑1 (MQ=255)
                                                                                                      ccTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGcc                                               >  1:1083037/1‑150 (MQ=255)
                                                                                                        ttCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCa                                             <  2:1772373/150‑1 (MQ=255)
                                                                                                           tGCAGGCCCTGCATGGCGACCTCCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCAccc                                          <  1:1513856/150‑1 (MQ=255)
                                                                                                            gCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCg                                         >  1:1694029/1‑150 (MQ=255)
                                                                                                             cAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGc                                        >  2:951147/1‑150 (MQ=255)
                                                                                                                gCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAg                                     >  2:1737263/1‑150 (MQ=255)
                                                                                                                 cccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGAG‑TAGTCCAGGCGTCGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAgg                                    <  2:50324/150‑1 (MQ=255)
                                                                                                                  ccTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAggg                                   >  2:138713/1‑150 (MQ=255)
                                                                                                                   cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                            <  2:1633754/92‑1 (MQ=255)
                                                                                                                   cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATca                                                                                            >  1:1633754/1‑92 (MQ=255)
                                                                                                                   cTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGc                                  <  1:138713/150‑1 (MQ=255)
                                                                                                                      cATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGccc                                                                           >  2:321309/1‑106 (MQ=255)
                                                                                                                      cATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGccc                                                                           <  1:321309/106‑1 (MQ=255)
                                                                                                                       aTGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc                              <  1:1001662/150‑1 (MQ=255)
                                                                                                                       aTGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc                              <  1:929307/150‑1 (MQ=255)
                                                                                                                               caccTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAg                      <  1:1737263/150‑1 (MQ=255)
                                                                                                                                accTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg                     >  1:123564/1‑150 (MQ=255)
                                                                                                                                accTGGGCACGCCGGTCGCCCTCGAAGGCGCGAGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAgg                     >  1:439834/1‑150 (MQ=255)
                                                                                                                                  cTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTc                   <  2:300604/150‑1 (MQ=255)
                                                                                                                                          gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                              <  2:1443282/99‑1 (MQ=255)
                                                                                                                                          gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTtgct                                                              >  1:1443282/1‑99 (MQ=255)
                                                                                                                                          gCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACa           >  1:240759/1‑150 (MQ=255)
                                                                                                                                                   ccTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATggg  >  2:533197/1‑150 (MQ=255)
                                                                                                                                                   ccTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATggg  >  2:342536/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACAGCGCATGGG  >  NC_002947/5741798‑5742094

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