Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,823 C→G 5.9% H89D (CAC→GAC)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,8230CG 5.9% 310.6 / 6.5 97H89D (CAC→GAC) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base C (47/44);  new base G (2/4);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAGGACTTTGCCGGCTGCTG  >  NC_002947/980678‑980971
                                                                                                                                                   |                                                                                                                                                    
cGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACg                                                                                                                                                   >  7:242030/1‑149 (MQ=255)
   cGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGAcc                                                                                                                                                <  3:184932/149‑1 (MQ=255)
    gCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGAccc                                                                                                                                               >  2:271121/1‑149 (MQ=255)
    gCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGAcc                                                                                                                                                >  6:162689/1‑148 (MQ=255)
     ccTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCt                                                                                                                                              >  1:68251/1‑149 (MQ=255)
         cGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGa                                                                                                                                            <  8:261441/147‑1 (MQ=255)
         cGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACt                                                                                                                                          >  6:254415/1‑149 (MQ=255)
           aGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTc                                                                                                                                        <  3:45099/149‑1 (MQ=255)
              gTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctg                                                                                                                                     >  2:115056/1‑149 (MQ=255)
               ttatcagcGTATTTCGCTGTAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAAGTCAGGGTCGTGGACCGTGACCTCctg                                                                                                                                     <  2:304701/144‑1 (MQ=255)
               tGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAAGTCAGGGTCGTGGACCGTGACCTCctgg                                                                                                                                    <  4:98558/149‑2 (MQ=255)
                gAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgct                                                                                                                                   >  7:289428/1‑149 (MQ=255)
                 agcagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgct                                                                                                                                   <  8:289428/148‑1 (MQ=255)
                   cagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgctg                                                                                                                                  >  3:269799/1‑147 (MQ=255)
                   cagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgca                                                                                                                                <  1:82632/149‑1 (MQ=255)
                    agcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcag                                                                                                                               <  4:269799/149‑1 (MQ=255)
                      cGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagca                                                                                                                             <  5:254415/149‑1 (MQ=255)
                             cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAg                                                                                                                      >  7:127558/1‑149 (MQ=255)
                              gcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGt                                                                                                                     >  7:12515/1‑149 (MQ=255)
                                  gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGGAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                                                                                                  <  1:115056/148‑1 (MQ=255)
                                  gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGCCTTCCTGCTGCAGCAGATGCAGTGGCg                                                                                                                 >  7:311374/1‑149 (MQ=255)
                                  gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                                                                                                  <  5:191582/148‑1 (MQ=255)
                                  gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCg                                                                                                                 >  3:264918/1‑149 (MQ=255)
                                  gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTAGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                                                                                              >  2:46899/1‑136 (MQ=255)
                                  gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTAGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                                                                                              <  1:46899/136‑1 (MQ=255)
                                  gAAAAATCCGCCGCCCTCAACGAAGCTTACCAGCCGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGTAT‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                                                                                                  <  1:267168/149‑1 (MQ=255)
                                   aaaaaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAGCTGCTGGCCATCAGCGGCGACG‑AA‑GTGCCCCAGGAGGTGACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGt                                                                                                                <  1:244102/149‑1 (MQ=255)
                                      aaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                                                                                             <  5:145904/149‑1 (MQ=255)
                                        tCCGCGGCCCTCAACGAGGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCCCGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAaga                                                                                                           >  5:342917/1‑149 (MQ=255)
                                        tCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                                                                                             >  5:107883/1‑147 (MQ=255)
                                        tCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGAGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGAGATCCTGCTGGAGCAGATGCAGTGGCGTGaa                                                                                                             >  7:138000/1‑147 (MQ=255)
                                         ccGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                                                                                          >  1:80130/1‑149 (MQ=255)
                                          cgcgGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AACTTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                                                                                          <  1:271121/149‑1 (MQ=255)
                                          cgcgGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                                                                                          <  5:162689/148‑1 (MQ=255)
                                                 cTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaaga                                                                                                  <  2:139542/149‑1 (MQ=255)
                                                  tCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaaga                                                                                                  >  4:201640/1‑148 (MQ=255)
                                                   cAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCt                                                                                                        >  4:223715/1‑141 (MQ=255)
                                                   cAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCt                                                                                                        <  3:223715/141‑1 (MQ=255)
                                                    aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACT‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACt                                                                                               <  8:301636/149‑1 (MQ=255)
                                                    aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAAc                                                                                                <  2:68251/148‑1 (MQ=255)
                                                      cgacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                                                                                              >  3:121144/1‑148 (MQ=255)
                                                         cgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAgg                                                                                          <  3:201640/149‑1 (MQ=255)
                                                            ttACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGAt                                                                                        >  2:99920/1‑148 (MQ=255)
                                                            ttACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGAt                                                                                        >  7:28386/1‑148 (MQ=255)
                                                               ccAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                                                                                              <  8:76618/107‑1 (MQ=255)
                                                               ccAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                                                                                              >  7:76618/1‑107 (MQ=255)
                                                               ccAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg                                                                                    >  4:274834/1‑149 (MQ=255)
                                                               ccAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg                                                                                    <  2:80130/149‑1 (MQ=255)
                                                                cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGc                                                                                   >  1:46378/1‑149 (MQ=255)
                                                                cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGc                                                                                   >  4:233883/1‑149 (MQ=255)
                                                                  gACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCg                                                                                 <  8:242030/149‑1 (MQ=255)
                                                                  gACGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg                                                                                    <  7:186352/146‑1 (MQ=255)
                                                                  gACGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg                                                                                    >  8:186352/1‑146 (MQ=255)
                                                                  gACGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGcc                                                                                  >  8:46471/1‑148 (MQ=255)
                                                                   aCGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa                                                                                >  1:296019/1‑149 (MQ=255)
                                                                      cTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCt                                                                             >  4:131488/1‑149 (MQ=255)
                                                                       tGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCtt                                                                            >  1:304768/1‑149 (MQ=255)
                                                                         cgcAGTGCGCCGCGTCGCGCGCGCTACCTGCTGCCCATCAGCGGCCACG‑AA‑GTGCCCCAGCAGGTCACGGTCCACGACCCTGACGTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGATCTGCAGGATGAAGCCGACCtt                                                                            >  3:273221/1‑147 (MQ=255)
                                                                          gcAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGAc                                                                         <  1:196914/149‑1 (MQ=255)
                                                                                cgccgcGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTTGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTc                                                                   >  3:241274/1‑149 (MQ=255)
                                                                                cgccgcGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTc                                                                   >  1:3065/1‑149 (MQ=255)
                                                                                 gccgcGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCg                                                                  >  4:38233/1‑149 (MQ=255)
                                                                                         gcgcgcgcTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATg                                                                                       <  8:104584/120‑1 (MQ=255)
                                                                                         gcgcgcgcTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATg                                                                                       >  7:104584/1‑120 (MQ=255)
                                                                                         gcgcgcgcTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGtt                                                           <  3:131488/148‑1 (MQ=255)
                                                                                           gcgcgcTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCaa                                                        <  1:99920/149‑1 (MQ=255)
                                                                                            cgcgcTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCaag                                                       <  4:264918/149‑1 (MQ=255)
                                                                                               gcTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTTGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCaagaag                                                    <  4:241274/149‑1 (MQ=255)
                                                                                                   ccTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAgcgc                                                 <  2:46378/148‑1 (MQ=255)
                                                                                                    ctgctgGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCt                                               <  2:304768/149‑1 (MQ=255)
                                                                                                       ctgGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGaa                                            <  3:274834/149‑1 (MQ=255)
                                                                                                           ccATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGgcc                                        <  4:121144/149‑1 (MQ=255)
                                                                                                            cATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGgccg                                       >  3:217542/1‑149 (MQ=255)
                                                                                                              tCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGaa                                            >  6:183514/1‑142 (MQ=255)
                                                                                                              tCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGaa                                            <  5:183514/142‑1 (MQ=255)
                                                                                                                     ccACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCaag                                                       >  8:314293/1‑124 (MQ=255)
                                                                                                                     ccACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCaag                                                       <  7:314293/124‑1 (MQ=255)
                                                                                                                      cACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGacac                             >  2:273646/1‑149 (MQ=255)
                                                                                                                        cG‑AA‑GTGCGCCAGTTGGTCACGGTCGACGACCCTGACTTCCTGCTGCAGCAGATGCAGTTGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGTAGCGCCTGTAGGCCGCCCAGGACACGc                           <  8:138000/149‑1 (MQ=255)
                                                                                                                        cG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGc                           <  8:276568/149‑1 (MQ=255)
                                                                                                                         g‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCt                          <  3:202633/149‑1 (MQ=255)
                                                                                                                         g‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCt                          >  4:273372/1‑149 (MQ=255)
                                                                                                                         g‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCt                          <  4:217542/149‑1 (MQ=255)
                                                                                                                         g‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCt                          <  8:127558/149‑1 (MQ=255)
                                                                                                                              gTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGaa                       >  3:211729/1‑149 (MQ=255)
                                                                                                                              ctgcagcAGGAAGTCAGGGTCGTGGACCGTGACCTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGaa                       >  8:102097/7‑149 (MQ=38)
                                                                                                                               tgcagcAGGAAGTCAGGGTCGTGGACCGTGACCTCCTGGGGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAAt                      >  2:20389/6‑149 (MQ=37)
                                                                                                                               tGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAAt                      <  3:233883/149‑1 (MQ=255)
                                                                                                                                 cagcAGGAAGTCAGGGTCGTGGACCGTGACCTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGa                    <  7:102097/146‑1 (MQ=38)
                                                                                                                                  cccAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCTCCTGAATGCCTCCCAGGACACGCTGAATGAg                   >  8:292150/1‑149 (MQ=255)
                                                                                                                                  cccAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAg                   <  1:273646/149‑1 (MQ=255)
                                                                                                                                     aggaggTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGATGGCCGCCCAGGACACGCTGAGTGacgcc                >  4:46973/1‑145 (MQ=255)
                                                                                                                                     aggaggTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAGGAc                >  8:100692/1‑149 (MQ=255)
                                                                                                                                             aCGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGa                                <  2:306891/125‑1 (MQ=255)
                                                                                                                                             aCGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGa                                >  1:306891/1‑125 (MQ=255)
                                                                                                                                             aCGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAGGACTTTGCCg         <  2:3065/148‑1 (MQ=255)
                                                                                                                                              cGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAGGACTTTGCCGgc       >  3:67622/1‑149 (MQ=255)
                                                                                                                                                   cACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGaa                       <  6:91837/128‑1 (MQ=255)
                                                                                                                                                   cACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGaa                       >  5:91837/1‑128 (MQ=255)
                                                                                                                                                   cACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAGGACTTTGCCGGctgctg  >  2:304123/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGTGTCGGCGTGTTCAAGAAGCGCCTGAAGGCCGCCCAGGACACGCTGAATGAGGACTTTGCCGGCTGCTG  >  NC_002947/980678‑980971

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