Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,060,489 T→C 6.5% intergenic (+125/+56) PP_5311 → / ← PP_5312 signal transduction protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,060,4890TC 6.5% 171.0 / 2.2 63intergenic (+125/+56)PP_5311/PP_5312signal transduction protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/33);  new base C (2/2);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

AGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCT  >  NC_002947/6060345‑6060636
                                                                                                                                                |                                                                                                                                                   
aGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGtt                                                                                                                                                 >  6:58698/1‑149 (MQ=255)
  ttCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTc                                                                                                                                               >  8:26689/1‑149 (MQ=255)
        tcACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCAGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCggg                                                                                                                                          >  6:225035/2‑148 (MQ=255)
                     cattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTg                                                                                                                                                   >  1:321298/1‑126 (MQ=255)
                     cattaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTg                                                                                                                                                   <  2:321298/126‑1 (MQ=255)
                       ttaaCACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCccgctccca                                                                                                                            <  1:11984‑M2/147‑10 (MQ=255)
                             caTAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGtt                                                                                                                                                 <  4:52001/120‑1 (MQ=255)
                             caTAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGtt                                                                                                                                                 >  3:52001/1‑120 (MQ=255)
                             caTAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTa                                                                                                                      >  1:71660/1‑147 (MQ=255)
                              aTAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCccgctcccacaaggtgat                                                                                                                   <  3:175504‑M2/149‑19 (MQ=255)
                               tAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGc                                                                                                                  <  2:11984/149‑1 (MQ=255)
                                  ccTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtc                                                                                                               <  3:330047/149‑1 (MQ=255)
                                   gatctggAGAAACTGTATTTCTGAATTTCTAACAGCTGGAGGGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtc                                                                                                               <  8:155339/143‑1 (MQ=21)
                                   cTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtca                                                                                                              >  8:205369/1‑149 (MQ=255)
                                    tGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTcgcg                                                                                                                                            <  3:25900/118‑1 (MQ=255)
                                    tGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTcgcg                                                                                                                                            >  4:25900/1‑118 (MQ=255)
                                                  tCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGc                                                                                               <  4:26890/149‑1 (MQ=255)
                                                       cTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAg                                                                                          >  7:100853/1‑149 (MQ=255)
                                                       cTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAg                                                                                          <  7:272300/149‑1 (MQ=255)
                                                        tCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGtca                                                                                                              <  7:300475/128‑1 (MQ=255)
                                                        tCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGc                                                                                         >  3:108984/1‑149 (MQ=255)
                                                          aaTTTCTACCAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCt                                                                                       <  3:255081/149‑1 (MQ=255)
                                                              tCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTttc                                                                                    <  4:169399/148‑1 (MQ=255)
                                                                tAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTttcttc                                                                                 >  8:256458/1‑149 (MQ=255)
                                                                    aGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTg                                                                             <  7:26689/149‑1 (MQ=255)
                                                                     gCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGGCGCTTGTTTCGCGGGCACGCCCGCGCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCCGCCTTTCTTCTTTgg                                                                            >  5:365493/1‑149 (MQ=255)
                                                                        ggCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGttt                                                                         >  1:145031/1‑149 (MQ=255)
                                                                          cgcgATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCa                                                                       <  4:108984/149‑1 (MQ=255)
                                                                           gcgATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAg                                                                      >  7:201804/1‑149 (MQ=255)
                                                                           gcgATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAg                                                                      >  3:279582/1‑149 (MQ=255)
                                                                              aTTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTa                                                                   >  7:7191/1‑149 (MQ=255)
                                                                                   ggACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCAGGTGCCATGCACCGCGTCGCTTGTTTCGCGGTCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGTTTTCAGGTATTtc                                                               <  5:225035/148‑1 (MQ=255)
                                                                                    gACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAg                                                                      >  1:330961/1‑140 (MQ=255)
                                                                                    gACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAg                                                                      <  2:330961/140‑1 (MQ=255)
                                                                                    gACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtca                                                              <  7:200704/148‑1 (MQ=255)
                                                                                    gACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtca                                                              <  5:60950/148‑1 (MQ=255)
                                                                                     aCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtcatc                                                            <  2:102465/149‑1 (MQ=255)
                                                                                     aCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtcatc                                                            <  5:58698/149‑1 (MQ=255)
                                                                                     aCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCGTGTTTAGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCGTTCGTCTTTGGTTTGAGGTATTtcaa                                                             >  4:348710/1‑147 (MQ=37)
                                                                                        gtgGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTACAGGTATTTCATGAAc                                                         >  3:241184/1‑149 (MQ=255)
                                                                                           ggAGGGGGCTGGCCCGCGAACACCGGAGAAGCCGGAGCCATGCTCCGCGTCGCGTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAAccc                                                       <  6:365493/148‑1 (MQ=37)
                                                                                            gAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCt                                                     <  2:71660/149‑1 (MQ=255)
                                                                                            gAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCt                                                     <  2:145031/149‑1 (MQ=255)
                                                                                                       cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCACGCccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactccct                                           <  2:293514‑M2/148‑91 (MQ=255)
                                                                                                         cgcgAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtcatc                                                            >  1:167592/1‑129 (MQ=255)
                                                                                                         cgcgAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTtcatc                                                            <  2:167592/129‑1 (MQ=255)
                                                                                                         cgcgAACACCGGCAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCACGCccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcat                                                   >  1:213822‑M2/1‑56 (MQ=255)
                                                                                                         cgcgAACACCGGCAAAGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCACGCccgctcccacagactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcat                                                   <  2:213822‑M2/138‑83 (MQ=255)
                                                                                                                acCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGcc                                 >  1:323295/1‑149 (MQ=255)
                                                                                                                  cGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCgg                               <  7:256458/149‑1 (MQ=255)
                                                                                                                  cGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCgg                               <  1:141149/149‑1 (MQ=255)
                                                                                                                       aaGCCGGTGCCATGCACCGCGGCGCCTGTTTCGCGGGCACGCccgctcccacacactgattgccattaagcgggcctatcgctacagtgacggcccgttttttgttgcctcgctaacgaatcatgactcccttagaacgctatcaagca                          >  1:293514‑M2/1‑42 (MQ=255)
                                                                                                                         gCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGccc                        >  8:10623/1‑149 (MQ=255)
                                                                                                                          ccGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCt                       <  8:7191/149‑1 (MQ=255)
                                                                                                                                ccATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGa                    >  7:105848/1‑146 (MQ=255)
                                                                                                                                ccATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATc                  <  2:323295/148‑1 (MQ=255)
                                                                                                                                ccATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCGTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGa                    <  8:105848/146‑1 (MQ=255)
                                                                                                                                 cATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTg                <  4:279582/149‑1 (MQ=255)
                                                                                                                                     cACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGAt             <  8:100853/148‑1 (MQ=255)
                                                                                                                                        cgcgTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGc         >  6:112781/1‑149 (MQ=255)
                                                                                                                                        cgcgTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGc         <  3:348710/149‑1 (MQ=255)
                                                                                                                                              gCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTcc   <  7:205369/149‑1 (MQ=255)
                                                                                                                                               cTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCt  >  4:75131/1‑149 (MQ=255)
                                                                                                                                                ttGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCt  >  5:233489/1‑148 (MQ=255)
                                                                                                                                                ttGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCt  >  1:19421/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGTTCCATGCACCGCATTAAGCATTAACACATAGCCTGCGGGAGCAACTGTCTTGCTCAATTTCTAACAGCTGGCGCGATTTTGGACTGTGGGAGCGGGCTTGCCCGCGAACACCGGCGAAGCCGGTGCCATGCACCGCGTCGCTTGTTTCGCGGGCACGCCCGCTCCCACAGGTAAAGCGTCATCATCCCCGTCGGGCTTCAGCCTTTCTTCTTTGGTTTCAGGTATTTCATCAACCCCTGGAACCACATCACCAGCGCCGGGTTGCCCTTGATCTGGATGTGCTTGTCCT  >  NC_002947/6060345‑6060636

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: