Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,741,940 G→C 6.0% D90E (GAC→GAG) ‡ PP_5039 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,741,9400GC 6.0% 258.1 / 7.1 84D90E (GAC→GAG) ‡PP_5039hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base C (2/3);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

AGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACA  >  NC_002947/5741791‑5742085
                                                                                                                                                     |                                                                                                                                                  
agcaTTCGATAACGCTCCACGCACAGGCTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGCCCCCCTGGGCACGACgg                                                                                                                                                    <  2:1013190/150‑1 (MQ=255)
  caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTc                                                                                                                                                  <  1:845279/150‑1 (MQ=255)
  caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTc                                                                                                                                                  <  1:155592/150‑1 (MQ=255)
  caTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTCCATGGCGCCAACCTGAGCACGACCGTc                                                                                                                                                  <  2:504468/150‑1 (MQ=255)
   aTTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCg                                                                                                                                                 >  2:902291/1‑150 (MQ=255)
    ttCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTATACAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGc                                                                                                                                                >  1:1452901/1‑150 (MQ=255)
       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                                                                                                                                        >  1:892605/1‑136 (MQ=255)
                          gTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAg                                                                                                                                        <  2:892605/136‑1 (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)
                                  aGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGc                                                                                                                                                >  1:1478169/1‑120 (MQ=255)
                                  aGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGc                                                                                                                                                <  2:1478169/120‑1 (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                                                                                                                                     >  2:1615552/1‑121 (MQ=255)
                                            ttCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGgcg                                                                                                                                     <  1:1615552/121‑1 (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:1280966/1‑150 (MQ=255)
                                                                                        gcagcaGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCgggg                                                           >  2:522324/1‑150 (MQ=255)
                                                                                           gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        <  2:1166953/150‑1 (MQ=255)
                                                                                           gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        <  2:274651/150‑1 (MQ=255)
                                                                                           gcagcTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGttt                                                        <  1:68905/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                                                                                                                            <  1:1519865/72‑1 (MQ=18)
                                                                                                      tCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGCGCGCCTTCGAGGGCGACCGGCgtg                                                                                                                            >  2:1519865/1‑72 (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:929307/150‑1 (MQ=255)
                                                                                                                              aTGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGc                     <  1:1001662/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)
                                                                                                                                                     |                                                                                                                                                  
AGCATTCGATAACGCTCCACGCACAGGTTCTGCCAGTCCAGCACTTCGGCATGCTCCAGTACCAGCAGCACGAACACGCGCTGGATCGGCAGCAGCTGGTAGTCCCAGCCCTTCTGCAGGCCCTGCATGGCGACCACCTGGGCACGCCGGTCGCCCTCGAAGGCGCGCGGCGTGTCG‑CGGTAGATCATGCGTGGCAGCTGGTCGAGCAAAATCAGCAGGCCCAGCCAGCCCTGCGGGGTTTGCTGCCACTCTTCGAGCCCACCCGCCAGGGCGTGCTCGACCAGGTCGCCAAACA  >  NC_002947/5741791‑5742085

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