Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,207,599 G→A 5.4% G344R (GGG→AGG)  PP_1950 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,207,5990GA 5.4% 270.7 / 2.6 91G344R (GGG→AGG) PP_1950hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/40);  new base A (2/3);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CAGCTTAAAGCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATCTCTGCCATTACGCCTGG  >  NC_002947/2207452‑2207740
                                                                                                                                                     |                                                                                                                                                   
cAGCTTAAAGCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGATGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTgg                                                                                                                                                     <  7:10544/148‑1 (MQ=255)
 atcTTAAAGCATGGGATCGCATAGCGTGCTTCCTTGGCGCAGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTggg                                                                                                                                                    <  1:306203/146‑1 (MQ=255)
    ttAAAGCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATACTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCAATCTTGCCTGGGGc                                                                                                                                                  >  1:162522/1‑147 (MQ=255)
      aaaGCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCa‑‑‑a                                                                                                                                             <  8:252733/147‑1 (MQ=255)
      aaaGCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGc                                                                                                                                           >  4:426058/1‑149 (MQ=255)
         gCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCGTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGATAAG‑GC‑GTGATGGGGCATACATCTTTCGGTTTCGGTCTCCGTTCTTGCGTGGGGCA‑‑‑AGCGGt                                                                                                                                        >  4:271731/1‑149 (MQ=255)
         gCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCACCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCCCAGGCA‑‑‑Agaatg                                                                                                                                        <  1:112968/149‑5 (MQ=255)
          cAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCg                                                                                                                                         >  6:216739/1‑147 (MQ=255)
          cAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATGCTTGCCCCAGGCA‑‑‑Agaatgg                                                                                                                                       <  5:33583/149‑6 (MQ=255)
           aGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAGGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCCCAGGCA‑‑‑Agaa                                                                                                                                          >  7:191104/1‑143 (MQ=255)
           aGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAGGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCCCAGGCA‑‑‑Agaa                                                                                                                                          <  8:191104/145‑3 (MQ=255)
             gggATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGcc                                                                                                                                    >  6:46689/1‑149 (MQ=255)
                aTCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑Gt                                                                                                                              <  1:214917/149‑1 (MQ=255)
                    cATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATg                                                                                                                           >  5:312214/1‑148 (MQ=255)
                     aTAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGa                                                                                                                         <  3:94088/149‑1 (MQ=255)
                      tAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAg                                                                                                                        <  3:194016/149‑1 (MQ=255)
                       aGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGt                                                                                                                       >  7:355717/1‑149 (MQ=255)
                        gCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAGCTCGGTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGa                                                                                                                         <  1:353353/148‑1 (MQ=255)
                        gCGTGCTGCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGtt                                                                                                                      <  3:269547/149‑1 (MQ=255)
                         cGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGACATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTAGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGACC‑‑‑GTATGGAGtt                                                                                                                      >  3:100820/1‑148 (MQ=255)
                           tgattCCTTGGCTCTGCATGACTTGATCCTATGTAGTTGTATCGCCATGACGAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCaa                                                                                                                   <  7:188788/146‑1 (MQ=255)
                                ccTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCaa                                                                                                              >  1:55213/1‑149 (MQ=255)
                                 cTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAAt                                                                                                             >  6:76247/1‑149 (MQ=255)
                                  ttGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAAt                                                                                                             >  7:195453/1‑148 (MQ=255)
                                    ggCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATc                                                                                                            <  4:265900/147‑1 (MQ=255)
                                       gcTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGtt                                                                                                       >  4:6814/1‑149 (MQ=255)
                                       gcTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGtt                                                                                                       >  4:732/1‑149 (MQ=255)
                                        cTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTg                                                                                                      >  8:25446/1‑149 (MQ=255)
                                           ccGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGcc                                                                                                   >  2:132219/1‑149 (MQ=255)
                                           ccGGACGTGATCCTATGAAGTGGGAGATCCTTGAGGAGTTGGATTTCACGCTTAATG‑GC‑GTGATTGGACTTACATCTTTCTGTTTCGGTCTCCATTCTTGACTGGGGCT‑‑‑AGCCGTATCCC‑‑‑GTATTGAGTTGATAGCAATCATGTTggac                                                                                                   <  2:162522/149‑3 (MQ=255)
                                             ggACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCtt                                                                                                 <  5:76247/149‑1 (MQ=255)
                                             ggACGTGAACCTATGAAGTGGGATCGCCATGACCAGTTCGATATAACGCGAAAAG‑TA‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCAATTCATGACTTAGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCAtt                                                                                                 >  6:296752/1‑149 (MQ=255)
                                               aCGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTgg                                                                                               <  4:186988/149‑1 (MQ=255)
                                                   gATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCGTGTTGGCCTTGGGTa                                                                                            >  7:319941/1‑148 (MQ=255)
                                                   gATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTa                                                                                            >  7:108660/1‑148 (MQ=255)
                                                    aTCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTaaa                                                                                          <  1:411332/149‑1 (MQ=255)
                                                       cTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAc                                                                                         >  8:336482/1‑147 (MQ=255)
                                                         aTGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAcgagt                                                                                     >  1:311982/1‑149 (MQ=255)
                                                          tGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAcgagt                                                                                     <  3:732/148‑1 (MQ=255)
                                                          tGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAcgagt                                                                                     <  3:6814/148‑1 (MQ=255)
                                                            aaGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAcgagtcga                                                                                  <  3:426058/149‑1 (MQ=255)
                                                             aGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAcgagtcgag                                                                                 >  4:157801/1‑149 (MQ=255)
                                                              gTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAAcgagtcgagt                                                                                <  3:264369/149‑1 (MQ=255)
                                                               tGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAgtg                                                                               <  3:271731/149‑1 (MQ=255)
                                                               tGGGATCGCCATGACCAGTTCGATATCACGCGAAAAC‑GC‑GTGATGGGACATACATGTTTCGGTTTAGGTCTCCATTATTTCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTAAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAgtg                                                                               >  1:70879/1‑149 (MQ=255)
                                                                 ggATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTa                                                                             >  1:39624/1‑149 (MQ=255)
                                                                 ggATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTa                                                                             >  3:47281/1‑149 (MQ=255)
                                                                    tCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTAtt                                                                           >  7:261605/1‑148 (MQ=255)
                                                                      gCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGa                                                                         <  5:46689/148‑1 (MQ=255)
                                                                       ccATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTgg                                                                                               <  8:285476/125‑1 (MQ=255)
                                                                       ccATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGATTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGt                                                                       >  3:283072/1‑149 (MQ=255)
                                                                       ccATGACCACTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTgg                                                                                               >  7:285476/1‑125 (MQ=255)
                                                                        cATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGtt                                                                      <  2:311982/149‑1 (MQ=255)
                                                                             ccAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAg                                                                 <  8:195453/149‑1 (MQ=255)
                                                                               aGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGAGAAAAGAGGCGTAGCCCGAAGAATGTAGTTAAAAACAAACATGTTTGAATTTGATAAAAGAGTCGAGTGTAttttattt                                                                     >  1:353939/1‑143 (MQ=255)
                                                                                  tCGATATCACGCGAAATG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCa                                                             <  8:355717/148‑1 (MQ=255)
                                                                                    gaaatCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGc                                                           <  7:192519/145‑1 (MQ=255)
                                                                                    gATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACTAGTCGAGTGTATTGAGTTGGTAGGTCAGc                                                           <  6:312214/148‑1 (MQ=255)
                                                                                    gATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGc                                                           <  1:132219/148‑1 (MQ=255)
                                                                                     atatCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTTGCTAAACGAGTCGAGTGTATTGAGTTGGt                                                                   <  1:112410/139‑1 (MQ=255)
                                                                                     atatCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGt                                                                   >  2:112410/1‑139 (MQ=255)
                                                                                     atatCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCt                                                         >  7:141212/1‑149 (MQ=255)
                                                                                      tatCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTg                                                        <  2:70879/149‑1 (MQ=255)
                                                                                        tCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGTAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGt                                                       <  1:321761/148‑1 (MQ=255)
                                                                                          aCGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑CGACGTAGCAC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCATGt                                                       >  5:286293/1‑146 (MQ=255)
                                                                                             cgAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGaaa                                                 <  5:216739/149‑1 (MQ=255)
                                                                                              gAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGaaa                                                 >  6:396718/1‑148 (MQ=255)
                                                                                                aaaG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGaaa                                                 <  3:93773/146‑1 (MQ=255)
                                                                                                aaaG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGaaa                                                 >  4:93773/1‑146 (MQ=255)
                                                                                                     gC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATcaac                                          <  2:39624/149‑1 (MQ=255)
                                                                                                     gC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATcaa                                           >  3:228314/1‑148 (MQ=255)
                                                                                                      c‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATcaac                                          >  5:190109/1‑148 (MQ=255)
                                                                                                        gTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATcaacaa                                        >  4:53620/1‑149 (MQ=255)
                                                                                                            tGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTAGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCaa                                    <  2:55213/149‑1 (MQ=255)
                                                                                                              ggACATACATCTTTCGGTTTCGGGCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGc                                  <  2:281939/149‑1 (MQ=255)
                                                                                                                   tacatCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTg                             >  8:358855/1‑149 (MQ=255)
                                                                                                                   tacatCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTg                             >  6:157865/1‑149 (MQ=255)
                                                                                                                    acatCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACa                                               >  2:194967/1‑130 (MQ=255)
                                                                                                                         tttcggtttcggTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGt                                                                       >  3:336878/1‑101 (MQ=255)
                                                                                                                         tttcggtttcggTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGt                                                                       <  4:336878/101‑1 (MQ=255)
                                                                                                                         tttcggtttcggTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTgg                        >  2:265845/1‑148 (MQ=255)
                                                                                                                           tcggtttcggtCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTg                         <  5:396718/145‑1 (MQ=255)
                                                                                                                             ggtttcggtCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAAt                    <  7:25446/148‑1 (MQ=255)
                                                                                                                               tttcggtCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAAtctc                 >  7:97127/1‑149 (MQ=255)
                                                                                                                                 tcggtCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATctct                <  8:108660/148‑1 (MQ=255)
                                                                                                                                     tctcCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATCTCTGCCAt           >  8:199825/1‑149 (MQ=255)
                                                                                                                                          aTTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCa                                                             >  2:160805/1‑94 (MQ=255)
                                                                                                                                          aTTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCa                                                             <  1:160805/94‑1 (MQ=255)
                                                                                                                                           ttCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAAt                                       <  5:428684/115‑1 (MQ=255)
                                                                                                                                           ttCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAAt                                       >  6:428684/1‑115 (MQ=255)
                                                                                                                                           ttCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATCTCTGCCATTACGcc     >  5:236567/1‑149 (MQ=255)
                                                                                                                                             cTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATCTCTGCCATTACGCCt    <  4:283072/148‑1 (MQ=255)
                                                                                                                                              ttGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATCTCTGCCATTACGCCTgg  >  3:218766/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGCTTAAAGCAGGGGATCGCATAGCGTGCTTCCTTGGCGCTGCCGGACGTGATCCTATGAAGTGGGATCGCCATGACCAGTTCGATATCACGCGAAAAG‑GC‑GTGATGGGACATACATCTTTCGGTTTCGGTCTCCATTCTTGCCTGGGGCA‑‑‑AGCCGTAGCCC‑‑‑GTATGGAGTTCAAAGCAATCATGTTGGCCTTGGCTAAACGAGTCGAGTGTATTGAGTTGGTAGGTCAGCCTGTCCGAAACATCAACAATCAAGCCAGTGGTTGGAAATCTCTGCCATTACGCCTGG  >  NC_002947/2207452‑2207740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑