Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,005 A→G 7.8% intergenic (+40/+27) PP_5374 → / ← PP_5375 choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,127,0050AG 7.8% 187.7 / 3.9 64intergenic (+40/+27)PP_5374/PP_5375choline/carnitine/betaine transporter family protein/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base G (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACA  >  NC_002947/6126856‑6127152
                                                                                                                                                      |                                                                                                                                                   
gTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCa                                                                                                                                                     <  1:102033/150‑1 (MQ=255)
     cgcCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGt                                                                                                                                                <  1:666723/150‑1 (MQ=255)
        cTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCgg                                                                                                                                             >  2:485737/1‑150 (MQ=255)
         ttCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTc                                                                                                                                               >  2:473439/1‑147 (MQ=255)
         ttCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTc                                                                                                                                               <  1:473439/147‑1 (MQ=255)
         ttCCTATACCTGGGATTCACCCGAAATTCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGt                                                                                                                                            >  1:996575/1‑150 (MQ=255)
              ataCCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTCATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTACAGCCACCGTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtgt                                                                                                                                       <  2:64977/150‑1 (MQ=255)
                  cTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGttt                                                                                                                                   <  1:338212/150‑1 (MQ=255)
                    gggATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAAACACAAGACCGACAGCTGATGGCTGTCGGTCTTGTGtcttc                                                                                                                                 <  2:712313/150‑2 (MQ=38)
                          cACCCGCAAGGCATCGCGCCGGTGGGTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAgg                                                                                                                           <  1:128531/150‑1 (MQ=255)
                          cACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtgt                                                                                                                                       >  1:772006/1‑138 (MQ=255)
                          cACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTtgt                                                                                                                                       <  2:772006/138‑1 (MQ=255)
                           aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCTtgt                                                                                                                                       <  2:1131551/137‑1 (MQ=38)
                           aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCTGATGGCTGTCGGTCTtgt                                                                                                                                       >  1:1131551/1‑137 (MQ=38)
                           aCCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGc                                                                                                                          >  1:797353/1‑150 (MQ=255)
                            cccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAgg                                                                                                                           <  2:572199/148‑1 (MQ=255)
                            cccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAgg                                                                                                                           >  1:572199/1‑148 (MQ=255)
                             ccGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTa                                                                                                                        >  1:94379/1‑150 (MQ=255)
                                 aaGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTc                                                                                                                    >  1:876124/1‑150 (MQ=255)
                                 aaGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACAGACAGCCATCAGCTGTGGGTGTTGTGTTTTGGGCGGGCCCTTTc                                                                                                                    >  2:283200/1‑150 (MQ=255)
                                   ggCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGACTGACAGCGAAGGAAGCCAGATAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGAAGGCTATTTCTg                                                                                                                  >  2:616660/1‑150 (MQ=255)
                                     cATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCa                                                                                                                <  1:517316/150‑1 (MQ=255)
                                      aTCACGCTGGGGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCGAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCAc                                                                                                               >  2:322572/1‑150 (MQ=255)
                                         aCGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGttt                                                                                                                                   <  2:286569/127‑1 (MQ=255)
                                         aCGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGttt                                                                                                                                   >  1:286569/1‑127 (MQ=255)
                                         aCGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCGGGCTATTTCTGCACCGa                                                                                                            >  2:900730/1‑150 (MQ=255)
                                                tGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAgg                                                                                                     <  1:154191/150‑1 (MQ=255)
                                                    tGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCa                                                                                                 >  1:157631/1‑150 (MQ=255)
                                                         ccTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                            <  1:485737/150‑1 (MQ=255)
                                                           ttGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCa                                                                                          <  1:488351/150‑1 (MQ=255)
                                                            tGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCAtatcaagc                                                                                         <  2:797353/150‑8 (MQ=255)
                                                                cTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTccc                                                                                     >  1:755850/1‑150 (MQ=255)
                                                                  gCGNGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCgg                                                                                   >  2:698033/1‑150 (MQ=255)
                                                                        ggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGc                                                                             <  2:901709/150‑1 (MQ=255)
                                                                        ggCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGc                                                                             <  1:605568/150‑1 (MQ=255)
                                                                           ttGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCa                                                                                          >  2:637344/1‑134 (MQ=255)
                                                                           ttGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCa                                                                                          <  1:637344/134‑1 (MQ=255)
                                                                           ttGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGc                                                                          >  2:470401/1‑150 (MQ=255)
                                                                            tggtatGGGGCTCACAGCGAAGAAAGCCAGCTAACCTAAAGCCACACTGGAACGAACACGCCTCGCAGGCTGCATCCGTAGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCa                                                                         <  1:78703/147‑1 (MQ=255)
                                                                                   gggCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTc                                                                  >  1:817342/1‑150 (MQ=255)
                                                                                        gACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAg                                                             <  2:755850/150‑1 (MQ=255)
                                                                                           aGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGAc                                                          >  2:720741/1‑150 (MQ=255)
                                                                                                agaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGc                                                     >  1:17283/1‑150 (MQ=255)
                                                                                                 gaaAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCt                                                    <  2:17283/150‑1 (MQ=255)
                                                                                                      ccAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTg                                               <  2:94379/150‑1 (MQ=255)
                                                                                                       cAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATGTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTGTCACGAGGACCACGCTCGGTgg                                              >  2:420004/1‑150 (MQ=255)
                                                                                                         gCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                            <  2:817342/150‑1 (MQ=255)
                                                                                                         gCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGgac                                            <  1:470401/150‑1 (MQ=255)
                                                                                                            aaCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAg                                                                                                                            >  2:38345/1‑67 (MQ=255)
                                                                                                            aaCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAg                                                                                                                            <  1:38345/67‑1 (MQ=255)
                                                                                                             aCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                            >  2:1003781/1‑98 (MQ=255)
                                                                                                             aCCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                            <  1:1003781/98‑1 (MQ=255)
                                                                                                              ccTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTc                                       >  2:405237/1‑150 (MQ=255)
                                                                                                                    gCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCcagcag                                 <  2:141908/150‑1 (MQ=255)
                                                                                                                     ccACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                 <  2:787351/117‑1 (MQ=255)
                                                                                                                     ccACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCa                                                                 >  1:787351/1‑117 (MQ=255)
                                                                                                                          cTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGaa                           >  1:105931/1‑150 (MQ=255)
                                                                                                                                gACTCAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTCTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                     <  2:653747/150‑1 (MQ=255)
                                                                                                                                gACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGc                     <  2:876124/150‑1 (MQ=255)
                                                                                                                                  cacaAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                            >  2:13489/1‑77 (MQ=21)
                                                                                                                                  cacaAGACCGACAGCTGATGGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAAc                                                                                            <  1:13489/77‑1 (MQ=21)
                                                                                                                                   acaAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTc                  >  1:492952/1‑150 (MQ=255)
                                                                                                                                     aaGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCgg                >  1:441750/1‑150 (MQ=255)
                                                                                                                                            aCAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAg         >  2:1013785/1‑150 (MQ=255)
                                                                                                                                                    tCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACcaca  >  1:98740/1‑150 (MQ=255)
                                                                                                                                                      aGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTccc                                                                                     <  2:814560/65‑1 (MQ=255)
                                                                                                                                                      aGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTccc                                                                                     >  1:814560/1‑65 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GTGGGCGCCTTCCTATACCTGGGATTCACCCGCAAGGCATCACGCCGGTGGCTGGGCCCTTGGCCTGCGAGGGGCTTGATATGGGGCTGACAGCGAAGAAAGCCAGCTAACCTAAAGCCACCCTGGAAGACACAAGACCGACAGCCATCAGCTGTCGGTCTTGTGTTTTGGGCAGGCTATTTCTGCACCGAAGGCAGGGCCATCAACCAGTCCCGGAACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACA  >  NC_002947/6126856‑6127152

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