Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 284,604:1 +G 8.1% intergenic (+50/+13) betT‑II → / ← atsK choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947284,6041.G 8.1% 240.6 / 14.9 62intergenic (+50/+13)betT‑II/atsKcholine/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base . (31/26);  new base G (2/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTGCCAACGCCAGCGAATG  >  NC_002947/284455‑284752
                                                                                                                                                      |                                                                                                                                                    
ccAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCa                                                                                                                                                       <  2:1066157/150‑1 (MQ=255)
 cAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACACGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAg                                                                                                                                                     >  2:820317/1‑150 (MQ=255)
       tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGc                                                                                                                                               >  1:1023993/1‑150 (MQ=255)
         tGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCTGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCAGGACAggaga                                                                                                                                             <  1:686117/150‑4 (MQ=37)
             cTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGc                                                                                                                                         >  1:16471/1‑150 (MQ=255)
                    gCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAg                                                                                                                                  <  2:1268530/150‑1 (MQ=255)
                     cTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAgg                                                                                                                                 <  1:1319830/150‑1 (MQ=255)
                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACGTGATCAGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                                                <  2:1046965/150‑1 (MQ=255)
                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                                                >  2:1122376/1‑150 (MQ=255)
                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                                                >  2:1164334/1‑150 (MQ=255)
                      tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACACGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                                                <  1:820317/150‑1 (MQ=255)
                       aCAGCAAGGAGCAGATCATCAACTACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTTACCAAGCTTTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTg                                                                                                                               <  1:1112001/150‑1 (MQ=255)
                                    gATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGAt                                                                                                                  >  1:735848/1‑150 (MQ=255)
                                     atcatcAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATg                                                                                                                 <  1:1164334/150‑1 (MQ=255)
                                       catcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGCCCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGCAGAGGCCAGACAGGCAACCGCTAGCAAGGTGCATCCCCAAGGATgct                                                                                                               >  2:190585/1‑148 (MQ=255)
                                       catcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAAcc                                                                                                                                           >  2:357174/1‑122 (MQ=255)
                                       catcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAAcc                                                                                                                                           <  1:357174/122‑1 (MQ=255)
                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAgg                                                                                                                                                >  2:1117802/1‑115 (MQ=255)
                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAgg                                                                                                                                                <  1:1117802/115‑1 (MQ=255)
                                         tcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCg                                                                                                             >  1:758546/1‑150 (MQ=255)
                                                          cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTACTGTCATGGGCTCTTAGCGGGATTGCCGGCGAGGAGGCCAGACAGGCAACGGCTAGAAAGTTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgc                                                                                            >  1:531338/1‑150 (MQ=255)
                                                                 aGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcg                                                                                     <  1:190585/150‑1 (MQ=255)
                                                                   cGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                                                   <  2:758546/150‑1 (MQ=255)
                                                                                 ccTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGGGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAc                                                                     >  2:1140351/1‑150 (MQ=255)
                                                                                  cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgcc                                                                                           <  2:766038/127‑1 (MQ=255)
                                                                                  cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgcc                                                                                           >  1:766038/1‑127 (MQ=255)
                                                                                  cTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCAcc                                                                    >  2:118941/1‑150 (MQ=255)
                                                                                       cttGGTCAGCGGACCAGGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCaaa                                                               >  2:966036/3‑150 (MQ=255)
                                                                                       ccTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCaaa                                                               >  1:79122/1‑150 (MQ=255)
                                                                                           gTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAccc                                                                                              <  2:253978/115‑1 (MQ=255)
                                                                                           gTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAccc                                                                                              >  1:253978/1‑115 (MQ=255)
                                                                                             ccGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCGCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTg                                                         <  2:1023993/150‑1 (MQ=255)
                                                                                             ccGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTg                                                         <  1:1122376/150‑1 (MQ=255)
                                                                                                   aCCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcg                                                                                 <  1:1248411/120‑1 (MQ=255)
                                                                                                   aCCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcg                                                                                 >  2:1248411/1‑120 (MQ=255)
                                                                                                     cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAcgc                                                 >  1:445568/1‑150 (MQ=255)
                                                                                                          tctcCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTg                                            >  1:669760/1‑150 (MQ=255)
                                                                                                             cctttcctgGCCTCTTCGCGGTCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGTTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                                         <  2:531338/150‑1 (MQ=255)
                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTgg                                                        <  2:886587/134‑1 (MQ=255)
                                                                                                              ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTgg                                                        >  1:886587/1‑134 (MQ=255)
                                                                                                               tgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCa                                       <  2:79122/150‑1 (MQ=255)
                                                                                                                  cctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGaa                                    <  2:445568/150‑1 (MQ=255)
                                                                                                                   ctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAAc                                   >  1:967142/1‑150 (MQ=255)
                                                                                                                   ctgGCCCCGTCGCGGCCTTGCCCGCGAAGAGGCCAGACAGGCAACCCCTAGAAAGGTGCATCCCCAAGGATGGTCACCCGGTGCATCATCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAAc                                   <  2:148914/150‑1 (MQ=255)
                                                                                                                      gCCTCTTCGCGGGCAAGCCCACGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGa                                                                                                                   >  1:637478/1‑67 (MQ=25)
                                                                                                                      gCCTCTTCGCGGGCAAGCCCACGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGa                                                                                                                   <  2:637478/67‑1 (MQ=25)
                                                                                                                        ctctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc                              <  2:735848/150‑1 (MQ=255)
                                                                                                                        ctctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGc                                                                                                     >  2:384322/1‑79 (MQ=37)
                                                                                                                        ctctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGc                                                                                                     <  1:384322/79‑1 (MQ=37)
                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGtt                                             >  2:281062/1‑135 (MQ=39)
                                                                                                                         tctTCGCGGGCAAGCCCGCGAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGtt                                             <  1:281062/135‑1 (MQ=39)
                                                                                                                          ctTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGt                            >  2:227205/1‑150 (MQ=255)
                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGaa                                    <  2:1338401/141‑1 (MQ=255)
                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGaa                                    >  1:1338401/1‑141 (MQ=255)
                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc                           >  1:203555/1‑150 (MQ=255)
                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCt                                                          >  2:955292/1‑120 (MQ=38)
                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCt                                                          <  1:955292/120‑1 (MQ=38)
                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGa                                                                                                                   >  1:244382/1‑62 (MQ=25)
                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGa                                                                                                                   <  2:244382/62‑1 (MQ=25)
                                                                                                                                   cTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                                       >  1:81572/1‑66 (MQ=255)
                                                                                                                                   cTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                                       <  2:81572/66‑1 (MQ=255)
                                                                                                                                    ttGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTg                  >  1:587571/1‑150 (MQ=255)
                                                                                                                                                  gCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTGCCAACGCCAGCGaa    >  1:359213/1‑150 (MQ=255)
                                                                                                                                                    cAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTACCAGAACGCCACGTCATTCTCCTGCCAACGCCAGCGAATg  >  1:670807/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCCTGCCAACGCCAGCGAATG  >  NC_002947/284455‑284752

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