Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,691,970 C→G 6.2% L23V (CTG→GTG)  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,9700CG 6.2% 201.0 / 3.7 65L23V (CTG→GTG) PP_1488methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

GACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCACCAGCCAGCAACTGGTGG  >  NC_002947/1691822‑1692117
                                                                                                                                                      |                                                                                                                                                   
gacAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCAtgctgc                                                                                                                                                     <  6:311186/149‑1 (MQ=255)
   aaTAACTTCAATACTGTAAGTACTCTTGTGTTTATTTCGTCTATTTTTCGCTCTCGTATTGTATTGCTGAGTTCTGAATGTGAAGAAGTGGACCTTGTGCCAACTGATCCTGGCAAGTTTCGCCGTTATCATCGGCATCATgcttctta                                                                                                                                                  <  8:325950/149‑3 (MQ=255)
         aTCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCg                                                                                                                                            >  7:275826/1‑149 (MQ=255)
         aTCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCg                                                                                                                                            >  4:243491/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                                                                                                                                          >  3:181667/1‑131 (MQ=255)
                             gggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTg                                                                                                                                          <  4:181667/131‑1 (MQ=255)
                              ggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTgg                                                                                                                       <  6:235489/149‑1 (MQ=255)
                              ggCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTg                                                                                                                        <  5:127310/148‑1 (MQ=255)
                               tcTATACTTCGTCGATTTTTGCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAATAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCAGCCTACTCGCGGTTGGt                                                                                                                      <  8:165940/148‑1 (MQ=255)
                               gCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCAGCAGCATGATGGCGATGATca                                                                                                                                            <  1:324696/127‑2 (MQ=255)
                               gCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCAGCAGCATGATGGCGATGATca                                                                                                                                            >  2:324696/1‑126 (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                                                                                                                     <  6:148853/107‑1 (MQ=255)
                                                                          ggCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTg                                                                                                                     >  5:148853/1‑107 (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                                  <  5:220059/149‑1 (MQ=255)
                                                                                                                   gCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGt                                  >  8:304824/1‑149 (MQ=255)
                                                                                                                   gCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGt                                  <  5:215062/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)
                                                                                                                                                    tgctgATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCACCAGCCAGCAACtggtg   <  4:174606/149‑1 (MQ=255)
                                                                                                                                                    tgctgATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCACCAGCCAGCAACtggtg   >  6:141915/1‑149 (MQ=255)
                                                                                                                                                      ctgATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCACCAGCCAGCAACtggtgg  <  2:1628/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAGTTTCGCCGTGATCATCGCCATCATGCTGCTGATGATCGTGGCCGCCTACTCGCGGTTGGTGACTATCGAAACCGCCGAAGAGGCGGTGGGGACCGACAGCATTCCCGGCGTCTACTACAGCTCGATGATCCGCAGCGCCTGGGTCGACAGCTACGTCACCAGCCAGCAACTGGTGG  >  NC_002947/1691822‑1692117

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: