Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,691,962 T→G 6.0% I20S (ATC→AGC)  PP_1488 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,691,9620TG 6.0% 196.5 / 4.3 67I20S (ATC→AGC) PP_1488methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base G (2/2);  total (33/34)
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.06e-01

GACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCAC  >  NC_002947/1691822‑1692100
                                                                                                                                              |                                                                                                                                          
gacAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCAtgctgc                                                                                                                                    <  6:311186/149‑1 (MQ=255)
   aaTAACTTCAATACTGTAAGTACTCTTGTGTTTATTTCGTCTATTTTTCGCTCTCGTATTGTATTGCTGAGTTCTGAATGTGAAGAAGTGGACCTTGTGCCAACTGATCCTGGCAAGTTTCGCCGTTATCATCGGCATCATgcttctta                                                                                                                                 <  8:325950/149‑3 (MQ=255)
         aTCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCg                                                                                                                           >  4:243491/1‑149 (MQ=255)
         aTCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCg                                                                                                                           >  7:275826/1‑149 (MQ=255)
                 cGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGgccgcc                                                                                                                   >  4:154677/1‑149 (MQ=255)
                   tAAGACCTATGGGATATACTTCGGCGATTTTACCCGCTCGGAGTGAATTACTGAGGCCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCTCGGTCATCATCTCCATCATGCTGCTGATGATCGTGgccggcaa                                                                                                                 >  6:201912/1‑145 (MQ=255)
                         ctctGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTcgcg                                                                                                           >  6:127310/1‑149 (MQ=255)
                            tGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTAAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCgg                                                                                                          <  5:211252/147‑1 (MQ=255)
                             gggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTg                                                                                                                         <  4:181667/131‑1 (MQ=255)
                             gggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTg                                                                                                                         >  3:181667/1‑131 (MQ=255)
                              ggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTgg                                                                                                      <  6:235489/149‑1 (MQ=255)
                              ggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTg                                                                                                       <  5:127310/148‑1 (MQ=255)
                               tcTATACTTCGTCGATTTTTGCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAATAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCAGCCTACTCGCGGTTGGt                                                                                                     <  8:165940/148‑1 (MQ=255)
                               gCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCAGCAGCATGATGGCGATGATca                                                                                                                           >  2:324696/1‑126 (MQ=255)
                               gCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCAGCAGCATGATGGCGATGATca                                                                                                                           <  1:324696/127‑2 (MQ=255)
                                   taCTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGa                                                                                                   <  4:273469/147‑1 (MQ=255)
                                   taCTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGa                                                                                                   <  5:47494/147‑1 (MQ=255)
                                        cGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCg                                                                                             <  1:152194/148‑1 (MQ=255)
                                         ggCGATTTTTCCCTCTCGTTTTGAATTCCTGATGGATGAATGTGAAGAACTGGACCTTTCGCCAACGGATCCTTGCAAGTTTCGCCGTTATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGa                                                                                            <  2:295886/148‑1 (MQ=255)
                                         ggCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGaa                                                                                           <  3:154677/149‑1 (MQ=255)
                                            gATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCACCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAAcc                                                                                        >  4:184040/1‑149 (MQ=255)
                                             aTTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAAccg                                                                                       >  2:31870/1‑149 (MQ=255)
                                                   cccGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGaa                                                                                  <  3:243491/148‑1 (MQ=255)
                                                      gCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAgg                                                                              >  3:151249/1‑149 (MQ=255)
                                                            tgtgAATTCCTGAGGGCTTAATGTTAACAACTCTCCCTTGCGCCAAGGGATCCTGGCCTATTTTCGCCGTTATTCATCCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGt                                                                          <  1:280286/148‑1 (MQ=255)
                                                                aaTTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGAc                                                                    >  6:125054/1‑149 (MQ=255)
                                                                aaTTCCTGAGGGCGGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCACCATGCTGCTGATGATCGGGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGAc                                                                    <  3:184040/149‑1 (MQ=255)
                                                                    ccTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGAc                                                                    <  2:143362/145‑1 (MQ=255)
                                                                    ccTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGAc                                                                    >  1:143362/1‑145 (MQ=255)
                                                                     cTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGa                                                                 >  4:130624/1‑147 (MQ=255)
                                                                      tGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAg                                                              >  3:174606/1‑149 (MQ=255)
                                                                      tGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAg                                                              >  1:259865/1‑149 (MQ=255)
                                                                        aGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAg                                                              <  2:283356/147‑1 (MQ=255)
                                                                          ggCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTg                                                                                                    >  5:148853/1‑107 (MQ=255)
                                                                          ggCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTg                                                                                                    <  6:148853/107‑1 (MQ=255)
                                                                                  gtgAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGt                                                  >  6:220059/1‑149 (MQ=255)
                                                                                  gtgAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGt                                                  >  6:215062/1‑149 (MQ=255)
                                                                                    gaagaaCTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTc                                                 >  8:36104/1‑148 (MQ=255)
                                                                                    gaagaaCTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTc                                                 >  8:36289/1‑148 (MQ=255)
                                                                                     aagaaCTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTcta                                               >  8:242907/1‑149 (MQ=255)
                                                                                       gaaCTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTg                                                                                                       >  1:343361/1‑91 (MQ=255)
                                                                                       gaaCTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTg                                                                                                       <  2:343361/91‑1 (MQ=255)
                                                                                          cTGGACCTTGCGCCAACGGATCCTGGCTAGTTTCGCCGTGATCATCGCAATCATGCTGCTTATTATCGTGGCCGCCTACTCGCTGTTGGTGACTATCGAAACCGCCGAAGAGGCTGTGGGGACCGACAGCATTCCCGGCGTCTACTACa                                          <  5:201912/149‑1 (MQ=255)
                                                                                          cTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTctacta                                            >  8:177601/1‑147 (MQ=255)
                                                                                          cTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTctacta                                            >  8:182546/1‑147 (MQ=255)
                                                                                          cTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTctacta                                            >  7:211584/1‑147 (MQ=255)
                                                                                               ccTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCg                                     <  7:36104/149‑1 (MQ=255)
                                                                                               ccTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCg                                     <  7:242907/149‑1 (MQ=255)
                                                                                                cTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCga                                    <  5:125054/149‑1 (MQ=255)
                                                                                                     gcCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATc                               <  7:182546/149‑1 (MQ=255)
                                                                                                       cAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATc                               <  7:177601/147‑1 (MQ=255)
                                                                                                         aCGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCa                           >  8:113782/1‑149 (MQ=255)
                                                                                                                   gCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGt                 >  8:304824/1‑149 (MQ=255)
                                                                                                                   gCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGt                 <  5:215062/149‑1 (MQ=255)
                                                                                                                   gCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGt                 <  5:220059/149‑1 (MQ=255)
                                                                                                                         tttCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGGCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTgg                   <  1:84263/142‑1 (MQ=255)
                                                                                                                         tttCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTgg                   >  2:84263/1‑142 (MQ=255)
                                                                                                                         tttCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGa              >  4:43857/1‑147 (MQ=255)
                                                                                                                           tCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAg           <  7:304824/148‑1 (MQ=255)
                                                                                                                              ccGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTAc       <  8:211584/149‑1 (MQ=255)
                                                                                                                              ccGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTAc       >  6:124777/1‑149 (MQ=255)
                                                                                                                               cGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGGCGACAGCTACg      <  3:130624/149‑1 (MQ=255)
                                                                                                                                gTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGt     <  7:36289/149‑1 (MQ=255)
                                                                                                                                gTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGt     >  1:1628/1‑149 (MQ=255)
                                                                                                                                 cgATCATCAGCAGCATGATGGCGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGc                                                    <  8:209706/100‑1 (MQ=38)
                                                                                                                                 cgATCATCAGCAGCATGATGGCGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGc                                                    >  7:209706/2‑101 (MQ=38)
                                                                                                                                   atcatcGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCAc  >  1:346067/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                          
GACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCAC  >  NC_002947/1691822‑1692100

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: