Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 672,892 T→A 9.2% I210F (ATC→TTC)  PP_0574 ← LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947672,8920TA 9.2% 180.5 / 8.2 65I210F (ATC→TTC) PP_0574LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base A (4/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

TCGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCACGGGCTTT  >  NC_002947/672751‑673033
                                                                                                                                                |                                                                                                                                             
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                                                     <  3:224386/115‑1 (MQ=255)
                                                                                                                      cGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGATTCTCGAACGCTGCTGGCAGCtcgctttc                                                     >  4:224386/1‑115 (MQ=255)
                                                                                                                       ggTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGc                  <  2:108887/149‑1 (MQ=255)
                                                                                                                         tAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCATTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCt                >  2:202386/1‑149 (MQ=255)
                                                                                                                         caattaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCggtgg                      >  3:323948/7‑143 (MQ=255)
                                                                                                                         caattaccAGATCGCCTG‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCggtgg                      <  4:323948/137‑1 (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)
                                                                                                                                       cctg‑‑CGAACTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCACGGGCttt  <  5:26770/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TCGTGAACCGCTGGTGCACGGGTGACGAGCTGGGGGTCAGCGCCAGAGCCAGCTGGGTGCGGTTGTGCATATGCGTCAAACGCAATACCTGAGACACATAGAGCTTCACGGTGTTTTCGGTAATGCCCAGTTCGCAGG‑‑CGATCTGGTAATTGGTCTTGCCCTTGCCTACCAGGCGCGCCACTTCCAGCTGGCGTGGAGACAGTTTCTCGAACGCTGCTGGCAGCTCGCTTTCGCTTTCTTCTGCATCACTGGCACGCCGGTGGGTGCCTTGGCCACGGGCTTT  >  NC_002947/672751‑673033

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: