Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 813,203 C→A 9.2% intergenic (+110/+21) PP_0699 → / ← PP_t70 LysE family transporter/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947813,2030CA 9.2% 196.7 / 8.9 65intergenic (+110/+21)PP_0699/PP_t70LysE family transporter/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base A (4/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

GGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCA  >  NC_002947/813056‑813350
                                                                                                                                                   |                                                                                                                                                   
ggCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCg                                                                                                                                                    <  5:52804/149‑1 (MQ=255)
       ccTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGcc                                                                                                                                             >  6:380400/1‑149 (MQ=255)
        cTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGccc                                                                                                                                            <  1:239841/149‑1 (MQ=255)
                 gCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTg                                                                                                                                   >  3:272105/1‑149 (MQ=255)
                       aCCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                             <  5:380400/149‑1 (MQ=255)
                       aCCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtgg                                                                                                                              <  2:91051/148‑1 (MQ=255)
                        ccGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggt                                                                                                                             <  4:272105/148‑1 (MQ=255)
                         cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtgg                                                                                                                           <  8:188320/149‑1 (MQ=255)
                         cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtgg                                                                                                                           <  3:359898/149‑1 (MQ=255)
                         cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtgg                                                                                                                           >  8:124305/1‑149 (MQ=255)
                         cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGCCCttt                                                                                                                                         <  6:61031/135‑1 (MQ=255)
                         cGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCGCAAGGTAAGCCCttt                                                                                                                                         >  5:61031/1‑135 (MQ=255)
                             aTACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTac                                                                                                                       <  7:124305/149‑1 (MQ=255)
                              tACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtg                                                                                                                            >  1:432854/1‑143 (MQ=255)
                              tACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTtggtg                                                                                                                            <  2:432854/143‑1 (MQ=255)
                                cGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTAcaca                                                                                                                    <  2:342702/149‑1 (MQ=255)
                                   tGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAgg                                                                                                                  >  8:103315/1‑148 (MQ=255)
                                      ggCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACt                                                                                                              <  3:222388/149‑1 (MQ=255)
                                       gCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACtt                                                                                                             >  4:346242/1‑149 (MQ=255)
                                         gacgaTAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTATTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGCCTTg                                                                                                            >  1:412255/1‑148 (MQ=255)
                                              tAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAAcccc                                                                                                      >  5:155382/1‑149 (MQ=255)
                                                 cagcGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGa                                                                                                   <  7:103315/149‑1 (MQ=255)
                                                      ccGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCa                                                                                              <  1:199267/149‑1 (MQ=255)
                                                         aaaCCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCa                                                                                           <  5:27594/149‑1 (MQ=255)
                                                            ccTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTa                                                                                        <  6:380941/149‑1 (MQ=255)
                                                            ccTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTa                                                                                        >  8:232349/1‑149 (MQ=255)
                                                              tgtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTAt                                                                                       >  2:33483/1‑148 (MQ=255)
                                                               gtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATg                                                                                      >  8:324165/1‑148 (MQ=255)
                                                               gtgATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATg                                                                                      <  7:324165/148‑1 (MQ=255)
                                                                 aaTTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAat                                                                                   >  4:162406/2‑149 (MQ=255)
                                                                  aTTCGGCGCCCGGCTGCTGCGGGATTGGTGGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAat                                                                                   <  5:228606/148‑1 (MQ=255)
                                                                   ttCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGAGAc                                                                                                               <  5:247014/119‑1 (MQ=255)
                                                                   ttCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGAGAc                                                                                                               >  6:247014/1‑119 (MQ=255)
                                                                       gcgcCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctatg                                                                             >  5:161236/1‑149 (MQ=255)
                                                                       gcgcCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAatgctatg                                                                             >  7:260612/1‑149 (MQ=255)
                                                                           ccGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGctct                                                                         <  7:212114/149‑1 (MQ=255)
                                                                                   tgcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                                 >  7:81973/1‑149 (MQ=255)
                                                                                    gcGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                                 >  5:407725/1‑148 (MQ=255)
                                                                                     cGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                                 <  6:155382/147‑1 (MQ=255)
                                                                                         aTTTGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTa                                                           <  1:370730/149‑1 (MQ=255)
                                                                                             gttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGAc                                                                                                               <  8:90938/93‑1 (MQ=255)
                                                                                             gttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                                 <  8:296009/139‑1 (MQ=255)
                                                                                             gttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACACCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACt                                                                 >  7:296009/1‑139 (MQ=255)
                                                                                              ttggCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGACTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAg                                                      >  8:7004/1‑149 (MQ=255)
                                                                                                       gCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGgcg                                             >  2:342966/1‑149 (MQ=255)
                                                                                                                  tttttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATtcg                                  >  4:245344/1‑148 (MQ=255)
                                                                                                                     ttCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTctct                                <  6:407725/148‑1 (MQ=255)
                                                                                                                           aaacaaGACAAAGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGc                         >  8:312385/5‑149 (MQ=38)
                                                                                                                            gCAAAGAAAAAGGGCTTACCTTGCGGTAAGCGCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCGATGTAGCCGATGGCGCGCAATATTCTCTTGAACGcc                        >  2:249349/1‑149 (MQ=255)
                                                                                                                              aaaGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGACCCCCGGACCCCAGCATTATGAATGGTATGCTCTAACCAACTGAGCTATTTAGCCGCTGGCGCGCATTATTCTCTTGAACGcccc                      >  3:308500/1‑149 (MQ=255)
                                                                                                                              aaaGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGAGCGCATTATTCGCTTGACCGcccc                      >  7:120089/1‑149 (MQ=255)
                                                                                                                              aaaGAAAAAGGGCTTACATTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGAACCCAGCATTATGAATGCTATGCTCTAACCAAATGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGcccc                      >  8:376491/1‑149 (MQ=255)
                                                                                                                                aGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCct                    <  3:346242/149‑1 (MQ=255)
                                                                                                                                 gacaaaGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGa                                                  >  7:376762/4‑118 (MQ=38)
                                                                                                                                 gacaaaGGGCTTACCGCAAGGTAAGCCCTTTTTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGa                                                  <  8:376762/115‑1 (MQ=38)
                                                                                                                                    aaaGGGCTTACCGCAAGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGt                >  8:348558/1‑149 (MQ=255)
                                                                                                                                      aGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTc               >  1:235098/1‑148 (MQ=255)
                                                                                                                                        ggCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCTTCTGTCAac            <  8:81973/149‑1 (MQ=255)
                                                                                                                                         gCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAaca           >  6:213810/1‑149 (MQ=255)
                                                                                                                                          cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCtt                               <  1:187861/128‑1 (MQ=255)
                                                                                                                                          cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCtt                               >  2:187861/1‑128 (MQ=255)
                                                                                                                                          cTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAacac          >  3:66419/1‑149 (MQ=255)
                                                                                                                                               cTTGCGGTAAGCCCTTTGTCTTGTTTGGTGACTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCAtt     <  7:7004/149‑1 (MQ=255)
                                                                                                                                                  gCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGTTTGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTc   <  7:232349/148‑1 (MQ=255)
                                                                                                                                                  gCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTc   <  1:33483/148‑1 (MQ=255)
                                                                                                                                                   cGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCa  <  2:235098/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCTTGGCCTGGCTGTGGCGGTGACCGGAATACGTTGAGGCGACGATAGCAGCGCCGAAACCTGTGATTCGGCGCCCGGCTGCTGCGGGATTGGTTGGCGCGGGCATCAGCAACTTTTTCGCAGGCAAAGAAAAAGGGCTTACCTTGCGGTAAGCCCTTTGTCTTGTTTGGTGGCTACACAGGGACTTGAACCCCGGACCCCAGCATTATGAATGCTATGCTCTAACCAACTGAGCTATGTAGCCGATGGCGCGCATTATTCTCTTGAACGCCCCCTCTGTCAACACTCATTTCA  >  NC_002947/813056‑813350

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: ‑