Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 672,887 G→T 8.7% A211A (GCC→GCA PP_0574 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947672,8870GT 8.7% 196.2 / 7.3 65A211A (GCC→GCAPP_0574LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (31/28);  new base T (4/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

ACTCAATGGGTCGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCACGGGCT  >  NC_002947/672741‑673031
                                                                                                                                                   |                                                                                                                                                  
aCTCAATGGGTCGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTAGCAGG‑‑c                                                                                                                                                >  7:189461/1‑149 (MQ=255)
          tCGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTaa                                                                                                                                      >  8:348866/1‑149 (MQ=255)
           cGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAAt                                                                                                                                     >  2:64802/1‑149 (MQ=255)
            gTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAAt                                                                                                                                     <  2:293185/148‑1 (MQ=255)
                aCCGCTTGTGCACGGGTGACGTTCTGGGTGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGAGTCAAACGCAATACCTGAGACACATAGAGCTTCACTGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGt                                                                                                                                 <  6:154657/148‑1 (MQ=255)
                aCCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTc                                                                                                                                >  6:135291/1‑149 (MQ=255)
                 ccGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTct                                                                                                                               <  7:330889/149‑1 (MQ=255)
                  cGCTGGTGCACGGGTGACGAGCTGGGAGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGt                                                                                                                                 >  1:3666/1‑147 (MQ=255)
                  cGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGt                                                                                                                                 <  2:3666/146‑1 (MQ=255)
                    cTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTcttgc                                                                                                                            <  4:331687/149‑1 (MQ=255)
                        tGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTcttgccctt                                                                                                                        >  3:54546/1‑149 (MQ=255)
                        tGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTCCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTcttgccctt                                                                                                                        <  4:211916/149‑1 (MQ=255)
                             gggTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTa                                                                                                                   >  1:134466/1‑149 (MQ=255)
                                  aCGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAgg                                                                                                              >  4:162411/1‑149 (MQ=255)
                                                      cAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGc                                                                                           >  5:326937/1‑149 (MQ=255)
                                                          gCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTgg                                                                                       <  3:162411/149‑1 (MQ=255)
                                                           ccAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGga                                                                                      <  7:348866/149‑1 (MQ=255)
                                                             aGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGgaga                                                                                    <  7:174933/149‑1 (MQ=255)
                                                               cTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTAc                                                                                                                  >  5:270902/1‑117 (MQ=255)
                                                               cTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTAc                                                                                                                  <  6:270902/117‑1 (MQ=255)
                                                                  gggtcGGGTGTGCATATGCGGGAACCGCAATACCTGAGACACAAAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAg                                                                                 <  8:309456/143‑1 (MQ=255)
                                                                       ggTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTtctc                                                                           >  6:266781/1‑148 (MQ=255)
                                                                             gCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACgct                                                                    >  8:108838/1‑149 (MQ=255)
                                                                               atatGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACgctgc                                                                  >  6:232310/1‑149 (MQ=255)
                                                                                tatGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACgctgct                                                                 <  4:54546/149‑1 (MQ=255)
                                                                                tatGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGGGGAGACAGTTTCTAGACCgctgct                                                                 >  4:15409/1‑149 (MQ=255)
                                                                                    cGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGc                                                              <  7:108838/148‑1 (MQ=255)
                                                                                     gTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGtt                                                                               >  4:122505/1‑130 (MQ=255)
                                                                                     gTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGtt                                                                               <  3:122505/130‑1 (MQ=255)
                                                                                     gTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAg                                                            <  3:99521/149‑1 (MQ=255)
                                                                                      tCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGAGTGGAGACAGTTTCTAGAACGCTGCTGGCAGc                                                           >  8:153684/1‑149 (MQ=255)
                                                                                        aaaCGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGGGTCGATCTGGTAATTGTTCTTGCCATTGCCTACAAGGCGCGCCACTTCAAGCTGGCGTGGAGACATTTTCTCGAACGATCCTgg                                                               >  1:108887/1‑145 (MQ=255)
                                                                                        aaaCGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCtc                                                         <  1:183474/149‑1 (MQ=255)
                                                                                              aaTACCTGATACAGATAGAGCTTCTCGGTGTTTTCGGTAATGCCCAGTTCGCTTG‑‑CGATCTGGTAATTTTTCTTTCCATTGCCTACCAGGCGCGCCACTTCCATCTTGCGTGGTGTCAGTTTCTCGAACGCTGCTTGCAGCtagctttc                                                   <  8:189461/149‑1 (MQ=255)
                                                                                              aaTACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGga                                                                                      >  3:139231/1‑114 (MQ=255)
                                                                                              aaTACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGga                                                                                      <  4:139231/114‑1 (MQ=255)
                                                                                                   cTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTcgctttcgcttt                                              <  5:266781/149‑1 (MQ=255)
                                                                                                        acacATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTtcttct                                         <  5:232310/149‑1 (MQ=255)
                                                                                                          acaTAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTg                                        >  7:47324/1‑148 (MQ=255)
                                                                                                          acaTAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGc                                       <  2:134466/149‑1 (MQ=255)
                                                                                                            aTAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTATTTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCtcgctttcg                                                  <  5:192390/136‑1 (MQ=255)
                                                                                                            aTAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCtcgctttcg                                                  >  6:192390/1‑136 (MQ=255)
                                                                                                            aTAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCAt                                     >  5:366284/1‑149 (MQ=255)
                                                                                                              agagCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCa                                   >  2:211010/1‑149 (MQ=255)
                                                                                                                  cTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTgg                               >  6:296195/1‑149 (MQ=255)
                                                                                                                          tgtTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCgg                       >  5:382439/1‑149 (MQ=255)
                                                                                                                            ttttCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCggtg                     >  7:174120/1‑149 (MQ=255)
                                                                                                                              ttCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGc                                       >  3:183531/1‑129 (MQ=255)
                                                                                                                              ttCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGc                                       <  4:183531/129‑1 (MQ=255)
                                                                                                                              ttCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCggtgg                    <  8:174120/148‑1 (MQ=255)
                                                                                                                                cGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGATTCTCGAACGCTGCTGGCAGCtcgctttc                                                   >  4:224386/1‑115 (MQ=255)
                                                                                                                                cGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGATTCTCGAACGCTGCTGGCAGCtcgctttc                                                   <  3:224386/115‑1 (MQ=255)
                                                                                                                                 ggTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGc                <  2:108887/149‑1 (MQ=255)
                                                                                                                                   tAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCATTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCt              >  2:202386/1‑149 (MQ=255)
                                                                                                                                   caattaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCggtgg                    <  4:323948/137‑1 (MQ=255)
                                                                                                                                   caattaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCggtgg                    >  3:323948/7‑143 (MQ=255)
                                                                                                                                   caattaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCt              >  1:184071/7‑149 (MQ=255)
                                                                                                                                     attaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTg            >  6:26770/5‑149 (MQ=255)
                                                                                                                                     attaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTg            >  6:108857/5‑149 (MQ=255)
                                                                                                                                      tGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTgg           >  5:154911/1‑149 (MQ=255)
                                                                                                                                       gCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGc          >  6:70734/1‑149 (MQ=255)
                                                                                                                                        cccAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGc          <  6:154911/148‑1 (MQ=255)
                                                                                                                                        cccAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGc          <  5:70734/148‑1 (MQ=255)
                                                                                                                                         ccAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCa        >  6:36999/1‑149 (MQ=255)
                                                                                                                                                gCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCACGGGCt  <  5:36999/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACTCAATGGGTCGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCACGGGCT  >  NC_002947/672741‑673031

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: