Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,030,762 C→T 12.5% P9L (CCG→CTG)  mgtA → ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,030,7620CT12.5% 153.7 / 12.1 64P9L (CCG→CTG) mgtAATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (5/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

ATGGATTGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCGAA  >  NC_002947/3030615‑3030910
                                                                                                                                                    |                                                                                                                                                    
aTGGATTGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCcgccgc                                                                                                                                                    >  1:1290479/1‑150 (MQ=255)
 tGGATTGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCTACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCa                                                                                                                                                   <  1:373099/150‑1 (MQ=255)
             gTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGca                                                                                                                                       >  1:894581/1‑150 (MQ=255)
             gTAGAAACTTCACAACCGACACAACTTCCGACTCTCGGCAGCCGACAAACAATTGATCCCTCCAGGTAGTTGGATTCACCGACCCGACATTCATTCGTTCACCGCCTTAGTGAAACGGCTGTGTGCTCCTGCCGCTGCGGCGGCacgaca                                                                                                                                        <  1:1390930/150‑6 (MQ=255)
               aGAAGGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGcaca                                                                                                                                     <  1:1310394/150‑1 (MQ=255)
               aGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGcaca                                                                                                                                     <  2:691171/150‑1 (MQ=255)
                gAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAg                                                                                                                                    >  1:751048/1‑150 (MQ=255)
                 aaaGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGctca                                                                                                                                     >  1:915692/1‑148 (MQ=255)
                 aaaGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGcaca                                                                                                                                     <  2:915692/148‑1 (MQ=255)
                 aaaGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGc                                                                                                                                   <  1:1043366/150‑1 (MQ=255)
                 aaaGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGc                                                                                                                                   >  1:532062/1‑150 (MQ=255)
                     ttCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACa                                                                                                                               <  2:162050/150‑1 (MQ=255)
                     ttCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACa                                                                                                                               >  1:209854/1‑150 (MQ=255)
                        gCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                            <  2:334149/150‑1 (MQ=255)
                          tGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACt                                                                                                                          <  2:1290479/150‑1 (MQ=255)
                             cGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTc                                                                                                                       >  1:1179377/1‑150 (MQ=255)
                                 acaACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                            >  1:963849/1‑141 (MQ=255)
                                 acaACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                            <  2:963849/141‑1 (MQ=255)
                                   aaCTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAAc                                                                                                                                <  1:191167/135‑1 (MQ=255)
                                   aaCTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAAc                                                                                                                                >  2:191167/1‑135 (MQ=255)
                                   aaCTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGCCTGTCAAATAc                                                                                                                 >  1:17889/1‑150 (MQ=255)
                                    aCTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGATCAAGGAAGTGAGAGCAGCCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATAcc                                                                                                                <  2:748839/150‑1 (MQ=255)
                                      ttCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGAGGCAGGAGCACAGCAACATGACTGTCAAATACCGt                                                                                                              >  1:1076956/1‑150 (MQ=255)
                                       tCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTa                                                                                                             >  2:1265113/1‑150 (MQ=255)
                                          gACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAaca                                                                                                          >  1:294250/1‑150 (MQ=255)
                                            ctctCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACAcc                                                                                                        <  2:651841/150‑1 (MQ=255)
                                             tctcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                            <  1:432760/129‑1 (MQ=255)
                                             tctcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGa                                                                                                                            >  2:432760/1‑129 (MQ=255)
                                               tcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCa                                                                                                                                                   <  1:714745/104‑1 (MQ=255)
                                               tcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCa                                                                                                                                                   >  2:714745/1‑104 (MQ=255)
                                               tcGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCa                                                                                                     >  1:1415693/1‑150 (MQ=255)
                                                  gCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCg                                                                                                  <  2:1415693/150‑1 (MQ=255)
                                                  gCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAACGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCg                                                                                                  >  1:1040474/1‑150 (MQ=255)
                                                        gacagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                            <  2:341732/150‑1 (MQ=255)
                                                         acagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAgg                                                                                           <  2:391008/150‑1 (MQ=255)
                                                         acagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAgg                                                                                           <  2:532062/150‑1 (MQ=255)
                                                          cagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGc                                                                                          <  2:294250/150‑1 (MQ=255)
                                                             acaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGa                                                                                       >  2:1343232/1‑150 (MQ=255)
                                                               aaCTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAc                                                                                     <  2:1076956/150‑1 (MQ=255)
                                                                   gTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCga                                                                                 <  2:209854/150‑1 (MQ=255)
                                                                      ggAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGAc                                                                              >  1:172437/1‑150 (MQ=255)
                                                                            aaCCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCAt                                                                                                    <  2:141813/122‑1 (MQ=255)
                                                                            aaCCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCAt                                                                                                    >  1:141813/1‑122 (MQ=255)
                                                                                       ttCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCTGCGGCGGCac                                                                                                                                            >  1:79212/1‑70 (MQ=255)
                                                                                       ttCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCTGCGGCGGCac                                                                                                                                            <  2:79212/71‑2 (MQ=255)
                                                                                           gCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGcgccgc                                                         >  1:558154/1‑150 (MQ=255)
                                                                                                       gCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAgcgcg                                             >  2:697328/1‑150 (MQ=255)
                                                                                                          ttCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCaa                                          >  1:1079672/1‑150 (MQ=255)
                                                                                                           tCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAAt                                         <  1:1343232/150‑1 (MQ=255)
                                                                                                             gTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                            <  2:938736/97‑1 (MQ=255)
                                                                                                             gTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                            >  1:938736/1‑97 (MQ=255)
                                                                                                                 aCCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAAAACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTg                                   >  1:1182495/1‑150 (MQ=255)
                                                                                                                  ccGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTgg                                  >  2:883971/1‑150 (MQ=255)
                                                                                                                   cGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGc                                 <  2:1315698/150‑1 (MQ=255)
                                                                                                                                  tCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCaa                  >  2:1394779/1‑150 (MQ=255)
                                                                                                                                    tgtgtg‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAAcc                >  1:1412341/1‑150 (MQ=255)
                                                                                                                                        tg‑tCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGa            >  1:966329/1‑150 (MQ=255)
                                                                                                                                        tg‑tCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGa            <  2:1079672/150‑1 (MQ=255)
                                                                                                                                         gaccttGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGa            >  2:1003200/6‑150 (MQ=255)
                                                                                                                                          ctcctGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGAt           >  2:381561/5‑150 (MQ=255)
                                                                                                                                           tcctGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTgg                                  >  1:1223396/4‑126 (MQ=255)
                                                                                                                                           tcctGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTgg                                  <  2:1223396/123‑1 (MQ=255)
                                                                                                                                           tcctGCCGCTGCGGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATg          >  2:535807/4‑150 (MQ=255)
                                                                                                                                            cGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGc         <  2:894581/150‑1 (MQ=255)
                                                                                                                                             gTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCGGGCGGTAAGCCTAGCCAGCCGGGATggc        >  1:556608/1‑148 (MQ=255)
                                                                                                                                                   ccgcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCGaa  <  2:966329/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGGATTGCGCGAGTAGAAAGTTCGCTGCCGACACAACTTCCGACTCTCGGCAGCGGACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTG‑TCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGATGCCAGCGAA  >  NC_002947/3030615‑3030910

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