Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 378,092 Δ1 bp 5.5% intergenic (+32/+33) 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,0920T. 5.5% 219.9 / 5.9 73intergenic (+32/+33)PP_0314/gbcAhypothetical protein/glycine‑betaine dioxygenase subunit
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base . (2/2);  total (35/38)
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.49e-01

CTCATTCGCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATA  >  NC_002947/377943‑378241
                                                                                                                                                     |                                                                                                                                                      
cTCATTCGCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGt                                                                                                                                                        >  2:537018/1‑150 (MQ=255)
 tCATTCGCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGACACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAGGCGCCATGCTTCGGTc                                                                                                                                                       <  1:759062/150‑1 (MQ=255)
       gCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGATATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCtgctgc                                                                                                                                                <  1:425052/150‑1 (MQ=255)
         cacaGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCg                                                                                                                                              >  2:294653/1‑150 (MQ=255)
              ttttGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCa                                                                                                                                         >  1:597065/1‑150 (MQ=255)
                 tGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGCAGCAGACCGAAGCATgg                                                                                                                                      <  2:710952/150‑1 (MQ=255)
                      gCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGt                                                                                                                                 >  2:1209743/1‑150 (MQ=255)
                      gCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGt                                                                                                                                 <  1:537018/150‑1 (MQ=255)
                      gCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGt                                                                                                                                 <  1:919873/150‑1 (MQ=255)
                           tGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATc                                                                                                                            >  2:734538/1‑150 (MQ=255)
                              gCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGt                                                                                                                         <  2:975292/150‑1 (MQ=255)
                                 ttACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGATCCCGAAGCATGGCGGGTGCATCAGTtac                                                                                                                      <  1:998408/150‑1 (MQ=255)
                                      gATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTg                                                                                                                 >  1:232319/1‑150 (MQ=255)
                                           gCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGCAGCAGACCGAAGCATGGCgcttt                                                                                                                                <  1:1277848/130‑5 (MQ=38)
                                           gCAACGTCGTCTGACGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGCAGCAGACCGAAGCATGGCgcttt                                                                                                                                >  2:1277848/1‑126 (MQ=37)
                                                 tcgtcTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCActt                                                                                                      >  2:1197733/1‑150 (MQ=255)
                                                   gtcTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATTAGTTAATACTGCGCCTGCActtct                                                                                                    >  1:984481/1‑150 (MQ=255)
                                                   gtcTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCActtct                                                                                                    >  2:262226/1‑150 (MQ=255)
                                                      tGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCtttt                                                                                                 >  2:155294/1‑149 (MQ=255)
                                                         cGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAgg                                                                                              <  1:262226/150‑1 (MQ=255)
                                                           aGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTa                                                                                            <  2:984481/150‑1 (MQ=255)
                                                               gcgcTATCCCCTACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGgc                                                                                        <  2:597065/150‑1 (MQ=255)
                                                               gcgcTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                               >  2:1278476/1‑143 (MQ=255)
                                                               gcgcTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                               <  1:1278476/143‑1 (MQ=255)
                                                                 gcTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGgcgc                                                                                      <  1:294653/150‑1 (MQ=255)
                                                                     tCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                              >  2:667224/1‑122 (MQ=255)
                                                                     tCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                              <  1:667224/122‑1 (MQ=255)
                                                                       ccATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                              >  2:413328/1‑120 (MQ=255)
                                                                       ccATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTgcgc                                                                                                              <  1:413328/120‑1 (MQ=255)
                                                                          tACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGc                                                                             <  1:1209743/150‑1 (MQ=255)
                                                                           aCAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGcc                                                                            >  2:1061139/1‑150 (MQ=255)
                                                                             aGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCggg                                                                                                                                  >  2:662420/1‑94 (MQ=255)
                                                                             aGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCggg                                                                                                                                  <  1:662420/94‑1 (MQ=255)
                                                                                cGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATgg                                                                                                                                      >  2:933652/1‑87 (MQ=255)
                                                                                cGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATgg                                                                                                                                      <  1:933652/87‑1 (MQ=255)
                                                                                   tGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGgtt                                                                    <  1:1061139/150‑1 (MQ=255)
                                                                                       ttATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                       >  1:791112/1‑143 (MQ=255)
                                                                                       ttATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTCCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                       <  2:791112/143‑1 (MQ=255)
                                                                                        tATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCa                                                               <  2:394525/150‑1 (MQ=255)
                                                                                        tATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCGGCCCCCAGGTTGTTCa                                                               <  2:232319/150‑1 (MQ=255)
                                                                                         aTGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGTGCCGGCTCTGCCCCCAGGTTGTTCAg                                                              >  2:301589/1‑150 (MQ=255)
                                                                                           gagaTACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCa                                                            >  1:923378/1‑150 (MQ=255)
                                                                                              aTACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTACCCCCAGGTTGTTCAGCAcgc                                                         >  1:1060715/1‑150 (MQ=255)
                                                                                                  gtacgGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTctgc                                                                                                                                                    <  1:1098149/56‑1 (MQ=255)
                                                                                                  gtacgGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTctgc                                                                                                                                                    >  2:1098149/1‑56 (MQ=255)
                                                                                                  gtacgGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                       >  2:98448/1‑132 (MQ=255)
                                                                                                  gtacgGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAg                                                                       <  1:98448/132‑1 (MQ=255)
                                                                                                  gtacgGGAGTACAGCCTTGAGACGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGccc                                                     >  2:1130926/1‑150 (MQ=255)
                                                                                                        gAGTCCAGCCTTGAGATGCAGTGGTGTGGCAGGGCGATGCTACGGTCTGCTGCCGAAGCATGGAGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGCCTCTGCCCCCAGGTTGTTCAGCGCGCGCCCGCGGTa                                               >  2:485904/1‑150 (MQ=255)
                                                                                                        gAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTACCCCCAGGTTGTTCAGCACGCGCCCGCTGTa                                               <  2:1060715/150‑1 (MQ=255)
                                                                                                           tACAGCCTTGAGACGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCa                                            <  1:1130926/150‑1 (MQ=255)
                                                                                                              aGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCttt                                                                                                  >  2:1217252/1‑93 (MQ=255)
                                                                                                              aGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCttt                                                                                                  <  1:1217252/93‑1 (MQ=255)
                                                                                                                  ttGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCTCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATg                                     >  2:1090900/1‑150 (MQ=255)
                                                                                                                     agaTGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAg                                  <  1:734538/150‑1 (MQ=255)
                                                                                                                         gCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTca                              <  1:468314/150‑1 (MQ=255)
                                                                                                                           aGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGTCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGCTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTcacc                            <  1:485904/150‑1 (MQ=255)
                                                                                                                               gtgtGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGc                                                                                                                               <  2:1159092/47‑1 (MQ=255)
                                                                                                                               gtgtGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGc                                                                                                                               >  1:1159092/1‑47 (MQ=255)
                                                                                                                               gtgtGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATgg                                                                                         >  1:578386/1‑85 (MQ=255)
                                                                                                                               gtgtGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATgg                                                                                         <  2:578386/85‑1 (MQ=255)
                                                                                                                                tgtgAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTcaccacacc                       <  1:912966/150‑1 (MQ=255)
                                                                                                                                  tgcacccgcCATGCTTCGGCAGCAGACCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCAAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCaa                     >  2:172416/7‑150 (MQ=38)
                                                                                                                                    aaaGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAAc                   >  1:274002/1‑150 (MQ=255)
                                                                                                                                    aaaGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAAc                   >  1:532613/1‑150 (MQ=255)
                                                                                                                                      aGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGTGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTc                 <  1:301589/150‑1 (MQ=255)
                                                                                                                                        cgcCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGt               <  1:859983/150‑1 (MQ=255)
                                                                                                                                         gcCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                               <  2:567768/69‑1 (MQ=255)
                                                                                                                                         gcCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAg                                                                                               >  1:567768/1‑69 (MQ=255)
                                                                                                                                          ccATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCTCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAGGTTCACCACACCAAACTCGTaa             <  1:155294/150‑1 (MQ=255)
                                                                                                                                            aTGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGt           >  1:507710/1‑150 (MQ=255)
                                                                                                                                                ttCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTc       >  2:477989/1‑150 (MQ=255)
                                                                                                                                                 tCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCAcgc                                                         <  2:994817/99‑1 (MQ=255)
                                                                                                                                                 tCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCAcgc                                                         >  1:994817/1‑99 (MQ=255)
                                                                                                                                                     tCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATa  >  2:430681/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CTCATTCGCCACAGTTTTGTCAGCTTATGAGCTTTACGGATTGGCAACGTCGTCTGCCGAGCGGCGCTATCCCATACAGCCGGTGGCTTATGAGATACGTACGGGAGTACAGCCTTGAGATGCAGTGGTGTGAAAGCGCCATGCTTCGGTCTGCTGCCGAAGCATGGCGGGTGCATCAGTTACTACTGCGCCTGCACTTCTTTCAGGTATGGCGCCGGCTCTGCCCCCAGGTTGTTCAGCACGCGCCCGCTGTACCAGTCAATGAAGTTCACCACACCAAACTCGTAAGTCTTCGAATA  >  NC_002947/377943‑378241

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: