Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,784,880 C→A 8.4% intergenic (+46/‑117) kefA → / → PP_5068 intermediate mechanosensitive channel protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,784,8800CA 8.4% 293.1 / 13.9 96intergenic (+46/‑117)kefA/PP_5068intermediate mechanosensitive channel protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (49/38);  new base A (4/4);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

TCGGTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCACCTTCGA  >  NC_002947/5784733‑5785025
                                                                                                                                                   |                                                                                                                                                 
tcGGTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCg                                                                                                                                                  >  2:104745/1‑149 (MQ=255)
  ggTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTc                                                                                                                                                   <  8:46982/146‑1 (MQ=255)
  ggTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCg                                                                                                                                                  <  4:46048/147‑1 (MQ=255)
  ggTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGc                                                                                                                                                 <  1:104745/148‑1 (MQ=255)
  ggTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAAGTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGcc                                                                                                                                                <  6:118532/149‑1 (MQ=255)
   gTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGccc                                                                                                                                               <  5:82027/149‑1 (MQ=255)
      cGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTc                                                                                                                                             <  1:298021/148‑1 (MQ=255)
      cGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTc                                                                                                                                             <  2:81982/148‑1 (MQ=255)
      cGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCg                                                                                                                                            >  3:24880/1‑149 (MQ=255)
        ccAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGtt                                                                                                                                          >  6:27132/1‑149 (MQ=255)
         cAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCGAAGCCCTCgcag                                                                                                                                         <  1:181370/149‑4 (MQ=255)
                    gtTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCGAAGCCCTCgca                                                                                                                                          >  5:366610/1‑135 (MQ=255)
                    gtTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCGAAGCCCTCgca                                                                                                                                          <  6:366610/137‑3 (MQ=255)
                     ttCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGATTCGCCCTCGTTCGCGGGTAAAccc                                                                                                                             >  8:322639/1‑149 (MQ=255)
                       ccTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCg                                                                                                                                                  >  6:108600/1‑126 (MQ=255)
                       ccTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCg                                                                                                                                                  <  5:108600/126‑1 (MQ=255)
                        cTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCt                                                                                                                          >  3:356954/1‑149 (MQ=255)
                          gatggatgCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTcc                                                                                                                        >  6:374972/1‑149 (MQ=255)
                          gatggatgCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCGAAGCCCTCgca                                                                                                                                          <  8:67082/131‑3 (MQ=255)
                          gatggatgCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCGAAGCCCTCgca                                                                                                                                          >  7:67082/1‑129 (MQ=255)
                            tggatgCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTAATGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCt                                                                                                                       >  1:51854/1‑148 (MQ=255)
                                tgCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATc                                                                                                                  <  4:356954/149‑1 (MQ=255)
                                 gCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCg                                                                                                                 >  3:98252/1‑149 (MQ=255)
                                  cGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGccc                                                                                                                                               >  1:16488/1‑118 (MQ=255)
                                  cGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGccc                                                                                                                                               <  2:16488/118‑1 (MQ=255)
                                  cGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCgg                                                                                                                <  4:24880/149‑1 (MQ=255)
                                   gAAGGGCAGCAATCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATTGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCggg                                                                                                               <  8:161890/149‑1 (MQ=255)
                                   gAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGccc                                                                                                                                               >  3:16594/1‑117 (MQ=255)
                                   gAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGccc                                                                                                                                               <  4:16594/117‑1 (MQ=255)
                                     atggaagcaagcagcaaTTGATTCCTATTGAACAGCCAAAATCTGATGTCACTGCTTCGGCCTGAATTGTCTCGAATACTTTATTGTGTGCATAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGgc                                                                                                              <  5:274427/143‑1 (MQ=255)
                                         cagcaagcagcaaTTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTgg                                                                                                         >  5:299382/1‑149 (MQ=255)
                                            caagcagcaaTTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCg                                                                                                      >  8:289695/1‑149 (MQ=255)
                                                 agcaaTTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCtgtg                                                                                                 <  5:27132/149‑1 (MQ=255)
                                                   caaTTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTggg                                                                                               >  3:366412/1‑149 (MQ=255)
                                                   caaTTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAAACCGATCCTAAATCGGGCGATGGCCGCTGTggg                                                                                               >  3:372258/1‑149 (MQ=255)
                                                        gATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTgg                                                                                                <  4:159446/143‑1 (MQ=255)
                                                        gATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTgg                                                                                                >  3:159446/1‑143 (MQ=255)
                                                         aTTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGgggcggg                                                                                         <  4:98252/149‑1 (MQ=255)
                                                         aTTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGgggcggg                                                                                         <  5:374972/149‑1 (MQ=255)
                                                         aTTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGgggcggg                                                                                         >  2:238557/1‑149 (MQ=255)
                                                          ttCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGgggcgggc                                                                                        >  5:57430/1‑149 (MQ=255)
                                                          ttCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGgggcggg                                                                                         <  4:15710/148‑1 (MQ=255)
                                                            ccGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCt                                                                                       <  4:355925/148‑1 (MQ=255)
                                                             cgatttaaCAGACAAAATCGGATTGCAATGCTTCTGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCATCGTTCGCGGGTAAACCCGCTCCTAAATAGGGCTCTGGCCGCTGTGGGGGCGGGCTc                                                                                      <  4:372258/142‑1 (MQ=255)
                                                             cGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTc                                                                                      <  4:366412/148‑1 (MQ=255)
                                                              gATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCAc                                                                                    >  1:357621/1‑149 (MQ=255)
                                                              gATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCAc                                                                                    <  2:357621/149‑1 (MQ=255)
                                                              gATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCAc                                                                                    <  3:164320/149‑1 (MQ=255)
                                                              gATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCACTGGCCGCTGTGGGGGCGGGCTCAc                                                                                    >  6:110927/1‑149 (MQ=255)
                                                                   aaCAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCcgc                                                                                >  8:1301/1‑148 (MQ=255)
                                                                      aGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGaaa                                                                            <  2:194033/149‑1 (MQ=255)
                                                                       gCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGCTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAg                                                                           <  5:381031/149‑1 (MQ=255)
                                                                         cAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCACTCGTTCGCGGGTAAACCCGATCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGgc                                                                         >  8:108079/1‑149 (MQ=255)
                                                                             aTCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGgccgc                                                                      >  7:18930/1‑148 (MQ=255)
                                                                                       aCTGCTTCGGCCTGAAGGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGccc                                                            >  6:45843/1‑148 (MQ=255)
                                                                                           cTTCGGCCTGAATTGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGACATCGTTCGCGGGTAAACCCGCTCCTAAAACGGGCGCAGGCCGCTGTGTGGGCTGGCTCACCCCCGAAACGCCGACAAACGGCCCCATTg                                                       >  8:344175/1‑149 (MQ=255)
                                                                                           cTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCAtt                                                        <  7:289695/148‑1 (MQ=255)
                                                                                           cTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTg                                                       <  1:238557/149‑1 (MQ=255)
                                                                                           cTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCCCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTg                                                       >  6:37610/1‑149 (MQ=255)
                                                                                            ttCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTgg                                                      >  3:387389/1‑149 (MQ=255)
                                                                                              cGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTg                                                       <  3:116184/146‑1 (MQ=255)
                                                                                              cGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTg                                                       >  4:116184/1‑146 (MQ=255)
                                                                                               ggCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCAc                                                   >  5:184404/1‑149 (MQ=255)
                                                                                                gCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCAcc                                                  >  6:43742/1‑149 (MQ=255)
                                                                                                  cTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCAtt                                                        >  2:180744/1‑141 (MQ=255)
                                                                                                  cTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCAtt                                                        <  1:180744/141‑1 (MQ=255)
                                                                                                    gAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGAtt                                              >  6:87072/1‑149 (MQ=255)
                                                                                                    gAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGAtt                                              >  4:129310/1‑149 (MQ=255)
                                                                                                    gAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGAtt                                              >  4:57822/1‑149 (MQ=255)
                                                                                                     aaTGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGACGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTc                                             >  3:280990/1‑149 (MQ=255)
                                                                                                       tGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTcc                                            >  4:137292/1‑148 (MQ=255)
                                                                                                         ggCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAg                                         <  5:43742/149‑1 (MQ=255)
                                                                                                         ggCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAg                                         <  3:318304/149‑1 (MQ=255)
                                                                                                            tCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAg                                      >  4:168560/1‑149 (MQ=255)
                                                                                                            tCGAATACTTTATTGTGTGCAAAATGGCTGCGAGAGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAg                                                                  >  2:38258/1‑121 (MQ=255)
                                                                                                            tCGAATACTTTATTGTGTGCAAAATGGCTGCGAGAGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAg                                                                  <  1:38258/121‑1 (MQ=255)
                                                                                                              gAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGTGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAg                                      <  5:98624/147‑1 (MQ=255)
                                                                                                              gAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAg                                      <  7:1301/147‑1 (MQ=255)
                                                                                                                    tttATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCg                              >  2:333274/1‑149 (MQ=255)
                                                                                                                      tATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCAAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGt                             <  4:356622/148‑1 (MQ=255)
                                                                                                                      tATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTACTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCCCGAAAGGCCGCCAAGCGGCCCCATTGGCACCCATTCCAAGGAGCACTCCCGt                             >  6:317467/1‑148 (MQ=255)
                                                                                                                          gtgtgCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGaaa                         <  8:18930/148‑1 (MQ=255)
                                                                                                                           tgtgCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGc                       >  8:192035/1‑149 (MQ=255)
                                                                                                                             tgCAAAATGGCTGCGAGGGCTTGGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGccc                     >  7:58979/1‑149 (MQ=255)
                                                                                                                                 aaaTGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGa                 <  7:272407/149‑1 (MQ=255)
                                                                                                                                  aaTGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGAc                >  7:145332/1‑149 (MQ=255)
                                                                                                                                     gcgaacGAGGGCGAAGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCaa             >  2:223531/6‑149 (MQ=255)
                                                                                                                                      gCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCaa             <  7:322639/148‑1 (MQ=255)
                                                                                                                                       gaacGAGGGCGAAGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACt           <  1:223531/146‑1 (MQ=255)
                                                                                                                                        tGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTc          >  3:354997/1‑149 (MQ=255)
                                                                                                                                         gCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTcc         >  8:89388/1‑148 (MQ=255)
                                                                                                                                         acGAGGGCGAAGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACAGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCa         >  6:194072/2‑149 (MQ=255)
                                                                                                                                          cGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCAc        >  6:26281/1‑149 (MQ=255)
                                                                                                                                              ggCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCACCTTc    >  2:222614/1‑149 (MQ=255)
                                                                                                                                              ggCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCACCTTc    >  3:99931/1‑149 (MQ=255)
                                                                                                                                                 ttCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCACCTTCGa  >  8:113304/1‑148 (MQ=255)
                                                                                                                                                 ttCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCACCTTCGa  >  4:251733/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCGGTACGCCAGGTAGAGGTGTTCCTGATGGATGCGAAGGGCAGCAAGCAGCAATTGATTCCGATGGAACAGCCAAAATCGGATGGCACTGCTTCGGCCTGAATGGGCTCGAATACTTTATTGTGTGCAAAATGGCTGCGAGGGCTTCGCCCTCGTTCGCGGGTAAACCCGCTCCTAAATCGGGCGCTGGCCGCTGTGGGGGCGGGCTCACCCGCGAAAGGCCGCCAAGCGGCCCCATTGGCACCGATTCCAAGGAGCAGTCCCGTGAAAGCCCTCGACCAACTCACCTTCGA  >  NC_002947/5784733‑5785025

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: