Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 378,097 T→A 5.7% intergenic (+37/+28) PP_0314 → / ← gbcA hypothetical protein/glycine‑betaine dioxygenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947378,0970TA 5.7% 210.2 / 3.2 70intergenic (+37/+28)PP_0314/gbcAhypothetical protein/glycine‑betaine dioxygenase subunit
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base A (2/2);  total (35/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

GCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGC  >  NC_002947/377950‑378246
                                                                                                                                                   |                                                                                                                                                      
gCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGATATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCtgctg‑c                                                                                                                                                     <  1:425052/150‑1 (MQ=255)
  cacaGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCg                                                                                                                                                   >  2:294653/1‑150 (MQ=255)
       ttttGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCa                                                                                                                                              >  1:597065/1‑150 (MQ=255)
          tGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGG‑CAGCAGACCGAAGCATgg                                                                                                                                           <  2:710952/150‑1 (MQ=255)
               gCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGt                                                                                                                                      <  1:919873/150‑1 (MQ=255)
               gCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGt                                                                                                                                      <  1:537018/150‑1 (MQ=255)
               gCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGt                                                                                                                                      >  2:1209743/1‑150 (MQ=255)
                    tGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATc                                                                                                                                 >  2:734538/1‑150 (MQ=255)
                       gCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGt                                                                                                                              <  2:975292/150‑1 (MQ=255)
                          ttACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGATC‑CCGAAGCATGGCGGGTGCATCAGTtac                                                                                                                           <  1:998408/150‑1 (MQ=255)
                               gATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTg                                                                                                                      >  1:232319/1‑150 (MQ=255)
                                    gCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGG‑CAGCAGACCGAAGCATGGCgcttt                                                                                                                                     <  1:1277848/130‑5 (MQ=38)
                                    gCAACGTCGTCTGACGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGG‑CAGCAGACCGAAGCATGGCgcttt                                                                                                                                     >  2:1277848/1‑126 (MQ=37)
                                          tcgtcTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCActt                                                                                                           >  2:1197733/1‑150 (MQ=255)
                                            gtcTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATTAGTTAATACTGCGCCTGCActtct                                                                                                         >  1:984481/1‑150 (MQ=255)
                                            gtcTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCActtct                                                                                                         >  2:262226/1‑150 (MQ=255)
                                               tGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCtttt                                                                                                      >  2:155294/1‑149 (MQ=255)
                                                  cGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAgg                                                                                                   <  1:262226/150‑1 (MQ=255)
                                                    aGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTa                                                                                                 <  2:984481/150‑1 (MQ=255)
                                                        gcgcTATCCCCTACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGgc                                                                                             <  2:597065/150‑1 (MQ=255)
                                                        gcgcTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                                    >  2:1278476/1‑143 (MQ=255)
                                                        gcgcTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                                    <  1:1278476/143‑1 (MQ=255)
                                                          gcTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGgcgc                                                                                           <  1:294653/150‑1 (MQ=255)
                                                              tCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                                   <  1:667224/122‑1 (MQ=255)
                                                              tCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                                   >  2:667224/1‑122 (MQ=255)
                                                                ccATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                                   >  2:413328/1‑120 (MQ=255)
                                                                ccATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                                   <  1:413328/120‑1 (MQ=255)
                                                                   tACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGc                                                                                  <  1:1209743/150‑1 (MQ=255)
                                                                    aCAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGcc                                                                                 >  2:1061139/1‑150 (MQ=255)
                                                                      aGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCggg                                                                                                                                       >  2:662420/1‑94 (MQ=255)
                                                                      aGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCggg                                                                                                                                       <  1:662420/94‑1 (MQ=255)
                                                                         cGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATgg                                                                                                                                           >  2:933652/1‑87 (MQ=255)
                                                                         cGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATgg                                                                                                                                           <  1:933652/87‑1 (MQ=255)
                                                                            tGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGgtt                                                                         <  1:1061139/150‑1 (MQ=255)
                                                                                ttATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                            >  1:791112/1‑143 (MQ=255)
                                                                                ttATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTCCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                            <  2:791112/143‑1 (MQ=255)
                                                                                 tATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCa                                                                    <  2:394525/150‑1 (MQ=255)
                                                                                 tATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCGGCCCCCAGGTTGTTCa                                                                    <  2:232319/150‑1 (MQ=255)
                                                                                  aTGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGTGCCGGCTCTGCCCCCAGGTTGTTCAg                                                                   >  2:301589/1‑150 (MQ=255)
                                                                                    gagaTACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCa                                                                 >  1:923378/1‑150 (MQ=255)
                                                                                       aTACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTACCCCCAGGTTGTTCAGCAcgc                                                              >  1:1060715/1‑150 (MQ=255)
                                                                                           gtacgGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                            <  1:98448/132‑1 (MQ=255)
                                                                                           gtacgGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                            >  2:98448/1‑132 (MQ=255)
                                                                                           gtacgGGAGTACAGCCTTGAGACGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGccc                                                          >  2:1130926/1‑150 (MQ=255)
                                                                                                 gAGTCCAGCCTTGAGATGCAGTGGTGTGGCAGGGCGATGCTACGGTCTGCTG‑CCGAAGCATGGAGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGCCTCTGCCCCCAGGTTGTTCAGCGCGCGCCCGCGGTa                                                    >  2:485904/1‑150 (MQ=255)
                                                                                                 gAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTACCCCCAGGTTGTTCAGCACGCGCCCGCTGTa                                                    <  2:1060715/150‑1 (MQ=255)
                                                                                                    tACAGCCTTGAGACGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCa                                                 <  1:1130926/150‑1 (MQ=255)
                                                                                                       aGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCttt                                                                                                       <  1:1217252/93‑1 (MQ=255)
                                                                                                       aGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCttt                                                                                                       >  2:1217252/1‑93 (MQ=255)
                                                                                                           ttGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCTCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATg                                          >  2:1090900/1‑150 (MQ=255)
                                                                                                              agaTGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAg                                       <  1:734538/150‑1 (MQ=255)
                                                                                                                  gCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTca                                   <  1:468314/150‑1 (MQ=255)
                                                                                                                    aGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGTCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGCTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTcacc                                 <  1:485904/150‑1 (MQ=255)
                                                                                                                        gtgtGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGc                                                                                                                                    >  1:1159092/1‑47 (MQ=255)
                                                                                                                        gtgtGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGc                                                                                                                                    <  2:1159092/47‑1 (MQ=255)
                                                                                                                        gtgtGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATgg                                                                                              >  1:578386/1‑85 (MQ=255)
                                                                                                                        gtgtGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATgg                                                                                              <  2:578386/85‑1 (MQ=255)
                                                                                                                         tgtgAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTcaccacacc                            <  1:912966/150‑1 (MQ=255)
                                                                                                                           tgcacccgcCATGCTTCGG‑CAGCAGACCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCAAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCaa                          >  2:172416/7‑150 (MQ=38)
                                                                                                                             aaaGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAAc                        >  1:274002/1‑150 (MQ=255)
                                                                                                                             aaaGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAAc                        >  1:532613/1‑150 (MQ=255)
                                                                                                                               aGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGTGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTc                      <  1:301589/150‑1 (MQ=255)
                                                                                                                                 cgcCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGt                    <  1:859983/150‑1 (MQ=255)
                                                                                                                                  gcCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                                    <  2:567768/69‑1 (MQ=255)
                                                                                                                                  gcCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                                    >  1:567768/1‑69 (MQ=255)
                                                                                                                                   ccATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCTCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAGGTTCACCACACCAAACTCGTaa                  <  1:155294/150‑1 (MQ=255)
                                                                                                                                     aTGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGt                >  1:507710/1‑150 (MQ=255)
                                                                                                                                         ttCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTc            >  2:477989/1‑150 (MQ=255)
                                                                                                                                          tCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCAcgc                                                              <  2:994817/99‑1 (MQ=255)
                                                                                                                                          tCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCAcgc                                                              >  1:994817/1‑99 (MQ=255)
                                                                                                                                              tCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATa       >  2:430681/1‑150 (MQ=255)
                                                                                                                                                   tg‑cCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGc  <  2:532613/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTG‑CCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATACGGGC  >  NC_002947/377950‑378246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑