Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,532,328 A→T 6.1% L285L (CTT→CTA) ‡ PP_3121 ← transcriptional regulator, LysR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,532,3280AT 6.1% 199.4 / 2.1 65L285L (CTT→CTA) ‡PP_3121transcriptional regulator, LysR family
Reads supporting (aligned to +/- strand):  ref base A (28/33);  new base T (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGAC  >  NC_002947/3532184‑3532474
                                                                                                                                                |                                                                                                                                                  
gCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAg                                                                                                                                                  >  4:369188/1‑147 (MQ=255)
gCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAg                                                                                                                                                  <  3:369188/147‑1 (MQ=255)
 cAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAg                                                                                                                                               >  7:358201/1‑149 (MQ=255)
        gcgACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTGCTTGAAGTGCttct                                                                                                                                        <  8:323053/149‑3 (MQ=255)
         cgACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTtgtg                                                                                                                                       >  4:142471/1‑149 (MQ=255)
           aCTTAGGCAGCCTTGCGTCCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGtgct                                                                                                                                     <  4:253694/149‑1 (MQ=255)
                ggCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCgg                                                                                                                                 <  2:376503/148‑1 (MQ=255)
                 gCAGCCTTGCGACCTGTGTCTTCGCACATGTACCAGGTGCCTTGGCCGATGGCTGCAACCCTCGTGCATGCCAGATTTACGTAAAGCATGTCACCTCCGGGATACGGGTATTGTAAAGAAGCACTTCAAGCAGTCCTTGTGCTGCGGTg                                                                                                                               >  6:100369/1‑149 (MQ=255)
                  cAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTgg                                                                                                                              <  7:112295/149‑1 (MQ=255)
                     ccTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGaaa                                                                                                                           <  1:198515/149‑1 (MQ=255)
                         gCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCtt                                                                                                                       >  3:163940/1‑149 (MQ=255)
                         gCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCt                                                                                                                        <  4:202382/148‑1 (MQ=255)
                          cGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCtt                                                                                                                       <  6:283041/148‑1 (MQ=255)
                          cGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCtt                                                                                                                       <  2:162535/148‑1 (MQ=255)
                           gACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTtgt                                                                                                                     <  8:358201/149‑1 (MQ=255)
                             ccTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGc                                                                                                                   <  3:158493/149‑1 (MQ=255)
                                gggTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCt                                                                                                                  <  2:193484/147‑1 (MQ=255)
                                gggTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCt                                                                                                                  >  1:193484/1‑147 (MQ=255)
                                   tCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCt                                                                                                              >  5:76973/1‑148 (MQ=255)
                                   tCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTg                                                                                                             <  3:142471/149‑1 (MQ=255)
                                   tCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCGTGTGCTGCGGTGGAAAGCTTGTGCTCGCTg                                                                                                             >  4:71728/1‑149 (MQ=255)
                                        gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTtg                                                                                                                      >  7:47787/1‑135 (MQ=255)
                                        gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTtg                                                                                                                      <  8:47787/135‑1 (MQ=255)
                                        gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAggg                                                                                                         >  5:322987/1‑148 (MQ=255)
                                        gCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGACTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAggg                                                                                                         >  2:17945/1‑148 (MQ=255)
                                         cccATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGc                                                                                                        >  5:7695/1‑148 (MQ=255)
                                            aTGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATc                                                                                                    <  4:132033/149‑1 (MQ=255)
                                                 ccAGCTGCCTTGGCCGATGGCTGCAACGCTCCTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGc                                                                                                <  4:27061/148‑1 (MQ=255)
                                                     gTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGc                                                                                                <  1:131480/144‑1 (MQ=255)
                                                     gTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGc                                                                                                >  2:131480/1‑144 (MQ=255)
                                                        ccTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATgcgc                                                                                        <  6:7695/149‑1 (MQ=255)
                                                                 aTGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACt                                                                               <  8:231111/149‑1 (MQ=255)
                                                                     cTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTgg                                                                           >  2:87050/1‑149 (MQ=255)
                                                                       gCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTggg                                                                          <  6:113831/148‑1 (MQ=255)
                                                                       gCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTggg                                                                          >  5:113831/1‑148 (MQ=255)
                                                                         aaCGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTc                                                                       <  6:76973/149‑1 (MQ=255)
                                                                            gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTc                                                                                   <  6:15808/134‑1 (MQ=255)
                                                                            gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTc                                                                                   >  5:15808/1‑134 (MQ=255)
                                                                            gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACtgt                                                                             <  4:198344/140‑1 (MQ=255)
                                                                            gCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACtgt                                                                             >  3:198344/1‑140 (MQ=255)
                                                                                 tgcctgccAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAgg                                                               <  1:87050/149‑1 (MQ=255)
                                                                                   cctgccAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGTAAAGCTTGTGCTCGCTGAGGGCAATCACGCCTATGCGCCGCTCGACTGTGGGCTCGACCAAGGTg                                                             <  5:100369/149‑1 (MQ=255)
                                                                                     tgccAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGAc                                                           >  7:24773/1‑149 (MQ=255)
                                                                                            tttACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCgg                                                    <  4:163940/149‑1 (MQ=255)
                                                                                                    aCGCATGTCACCTCCCGGATCCGGGTATTGGAAAGAACCACTTCAACCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAg                                            <  6:322987/149‑1 (MQ=255)
                                                                                                               cTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCagat                                 >  1:100604/1‑146 (MQ=255)
                                                                                                               cTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCt                                 >  3:158442/1‑149 (MQ=255)
                                                                                                                     gATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGa                           <  1:145634/149‑1 (MQ=255)
                                                                                                                                 ggAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcg               >  1:69853/1‑149 (MQ=255)
                                                                                                                                  gAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcgg              >  4:169863/1‑149 (MQ=255)
                                                                                                                                  gAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcgg              >  2:36868/1‑149 (MQ=255)
                                                                                                                                  gAAAGAAGCACTTCAAGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCAGATGGTTGATGCACAGGGCcgg              >  2:160324/1‑149 (MQ=255)
                                                                                                                                    aagaagCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCCTCTGGTTGAAGCACAGGGCcggca            >  6:264528/1‑149 (MQ=255)
                                                                                                                                     agaagCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTTCATCTGGTTGATGCACAGGGCcggcac           <  1:17945/149‑1 (MQ=255)
                                                                                                                                      caagCACTGCTTGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGAt                          <  7:90775/132‑1 (MQ=255)
                                                                                                                                      caagCACTGCTTGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGAt                          >  8:90775/2‑133 (MQ=255)
                                                                                                                                       aagCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGACGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCTGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcggcacgg         <  1:160324/149‑1 (MQ=255)
                                                                                                                                       aagCACTGCTTGAAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCcggcacgg         >  1:233043/1‑149 (MQ=255)
                                                                                                                                             ttCAAGCAGTGCTTGTGCTGCGGTGTAAAGCTTGTGCTCTCTGATGGCAATCACGCCGATTCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGtt                             <  5:230934/123‑1 (MQ=255)
                                                                                                                                             ttCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCTGGCCCCCAGCTCCTGCATCttgat                             >  6:230934/1‑119 (MQ=255)
                                                                                                                                             ttCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCt                                          <  4:87852/110‑1 (MQ=255)
                                                                                                                                             ttCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCt                                          >  3:87852/1‑110 (MQ=255)
                                                                                                                                               cAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCAcc                                               <  4:349308/103‑1 (MQ=255)
                                                                                                                                               cAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCAcc                                               >  3:349308/1‑103 (MQ=255)
                                                                                                                                               cAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGAc  <  4:158442/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCAGCGCCGCGACTTCGGCAGCCTTGCGACCTGGGTCTTCGCCCATGTACCAGGTGCCTTGGCCGATGGCTGCAACGCTCGTGCCTGCCAGATTTACGTAACGCATGTCACCTCCGGGATCCGGGTATTGGAAAGAAGCACTTCAAGCAGTGCTTGTGCTGCGGTGGAAAGCTTGTGCTCGCTGAGGGCAATCACGCCGATGCGCCGCTCGACTGTGGGCTCGACCAAGGTGACGCAGCGGGCACCCAGCTCCTGCATCTGGTTGATGCACAGGGCCGGCACGGCACTGAC  >  NC_002947/3532184‑3532474

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: