Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,295,783 A→G 6.2% W180R (TGG→CGG)  ccmC ← heme export ABC transporter permease; CcmE‑interacting protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,295,7830AG 6.2% 191.5 / 2.6 64W180R (TGG→CGG) ccmCheme export ABC transporter permease; CcmE‑interacting protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAA  >  NC_000913/2295637‑2295931
                                                                                                                                                  |                                                                                                                                                    
cATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTccac                                                                                                                                                     >  5:480704/1‑148 (MQ=255)
cATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTccac                                                                                                                                                     >  2:51282/1‑148 (MQ=255)
 aTCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTccacca                                                                                                                                                   >  2:132546/1‑149 (MQ=255)
    aGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTc                                                                                                                                                >  3:204767/1‑149 (MQ=255)
       aaaaTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGCACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCAc                                                                                                                                             >  5:474783/1‑149 (MQ=255)
         aaTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGCCAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACg                                                                                                                                            >  5:153171/1‑148 (MQ=255)
         aaTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGCTCCCTGATGCAGGGTGTTCCACCACTCCACg                                                                                                                                            >  3:372957/1‑148 (MQ=255)
         aaTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCGCCACTCCACg                                                                                                                                            >  7:499948/1‑148 (MQ=255)
         aaTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCGCCACTCCACg                                                                                                                                            >  3:422133/1‑148 (MQ=255)
          aTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACg                                                                                                                                            >  7:378369/1‑147 (MQ=255)
          aTCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGa                                                                                                                                          <  1:210927/149‑1 (MQ=255)
             aaaTTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCCCGGAGt                                                                                                                                        <  7:102249/148‑1 (MQ=255)
               aTTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGTATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAAt                                                                                                                                     <  4:45184/149‑1 (MQ=255)
               aTTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAAt                                                                                                                                     <  3:93035/149‑1 (MQ=255)
                   cgcatccgcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGaa                                                                                                                                  <  5:292807/148‑1 (MQ=255)
                      atccgcagcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGTTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGaataat                                                                                                                              <  5:326848/149‑1 (MQ=255)
                         cgcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCGCCACTCCACGGAGTAATGAATAATCg                                                                                                                            <  8:499948/148‑1 (MQ=255)
                         cgcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCGCCACTCCACGGAGTAATGAATAATCg                                                                                                                            <  4:422133/148‑1 (MQ=255)
                         cgcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCgg                                                                                                                           <  3:295933/149‑1 (MQ=255)
                         cgcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCgg                                                                                                                           <  5:408347/149‑1 (MQ=255)
                          gcatcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCgg                                                                                                                           >  7:99102/1‑148 (MQ=255)
                           catcAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCa                                                                                                                         <  4:204767/149‑1 (MQ=255)
                                  gtcagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcac                                                                                                                  <  7:203554/149‑1 (MQ=255)
                                    cagcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcacca                                                                                                                >  2:476638/1‑149 (MQ=255)
                                      gcgtGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACg                                                                                                              >  8:152964/1‑149 (MQ=255)
                                          ggCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCaa                                                                                                          >  6:92756/1‑149 (MQ=255)
                                            cAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATc                                                                                                        <  5:92756/149‑1 (MQ=255)
                                              gACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAg                                                                                                      <  1:132546/149‑1 (MQ=255)
                                               aCAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGc                                                                                                     <  1:33466/149‑1 (MQ=255)
                                                      cAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAAt                                                                                                                                     <  3:63260/110‑1 (MQ=255)
                                                      cAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAAt                                                                                                                                     >  4:63260/1‑110 (MQ=255)
                                                        ggAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATa                                                                                            <  6:153171/149‑1 (MQ=255)
                                                             ccAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCCTCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGc                                                                                       >  7:83397/1‑149 (MQ=255)
                                                               aaaaaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCg                                                                                     >  2:502495/1‑149 (MQ=255)
                                                                 aaaTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCa                                                                                   <  5:415979/149‑1 (MQ=255)
                                                                   aTCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACg                                                                                 <  7:298241/149‑1 (MQ=255)
                                                                           gcgcagcgGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGc                                                                           >  3:27992/1‑147 (MQ=255)
                                                                           gcgcagcgGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCa                                                                         <  8:99102/149‑1 (MQ=255)
                                                                                         gCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGgtcgtc                                                            >  4:333749/1‑148 (MQ=255)
                                                                                           aTCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGtcgtcg                                                           >  6:194500/1‑147 (MQ=255)
                                                                                            tCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAg                                                        <  8:378369/149‑1 (MQ=255)
                                                                                               ccGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCg                                                     <  8:83397/149‑1 (MQ=255)
                                                                                                cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcaccac                                                                                                               >  6:433416/1‑90 (MQ=255)
                                                                                                cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTcaccac                                                                                                               <  5:433416/90‑1 (MQ=255)
                                                                                                cGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGt                                                    <  1:51282/149‑1 (MQ=255)
                                                                                                 ggATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTg                                                   >  2:31404/1‑149 (MQ=255)
                                                                                                   aTCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGc                                                  >  8:131087/1‑148 (MQ=255)
                                                                                                       aTACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAg                                             <  6:474783/149‑1 (MQ=255)
                                                                                                        tACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAgg                                            >  1:329491/1‑149 (MQ=255)
                                                                                                          cTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAggg                                           <  6:480704/148‑1 (MQ=255)
                                                                                                             tgctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCAAGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCaa                                        <  1:388277/148‑1 (MQ=255)
                                                                                                               ctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAg                                                                                                                                         <  4:78715/49‑1 (MQ=255)
                                                                                                               ctgcATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAg                                                                                                                                         >  3:78715/1‑49 (MQ=255)
                                                                                                                  cATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCacac                                  >  1:400126/1‑149 (MQ=255)
                                                                                                                       gcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCcaca                             <  5:194500/149‑1 (MQ=255)
                                                                                                                       gcgTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCcaca                             >  4:63333/1‑149 (MQ=255)
                                                                                                                          ttGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATAc                          <  3:63333/149‑1 (MQ=255)
                                                                                                                                ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCCCGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAc                    >  1:166588/1‑149 (MQ=255)
                                                                                                                                ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACaaaaa                     <  3:333749/148‑1 (MQ=255)
                                                                                                                                ccTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACAAACaaaaa                     >  5:309651/1‑148 (MQ=255)
                                                                                                                                ccTGATGCAGGGTGTACCACCACTCCACGGAGTAATGAATAAACCGCAGATACACCACGCCAATAAGCACCAGGATACCTGCCCCAAGTCCCGCCACACGGCGGTCGTAGAAGGCGTGCCACAGCGAAATCACACCCACATACaaaaa                     >  4:293516/1‑148 (MQ=255)
                                                                                                                                 cTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCcaca                             <  3:343898/139‑1 (MQ=255)
                                                                                                                                 cTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCcaca                             >  4:343898/1‑139 (MQ=255)
                                                                                                                                   gATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAAcag                  >  5:350324/1‑148 (MQ=255)
                                                                                                                                       cAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAgcagca             >  6:370195/1‑149 (MQ=255)
                                                                                                                                              ttCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACATGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt       <  1:502495/148‑1 (MQ=255)
                                                                                                                                               tCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGtt       >  4:222223/1‑147 (MQ=255)
                                                                                                                                                 caccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAg    <  7:131087/148‑1 (MQ=255)
                                                                                                                                                 caccacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCCTGCCACAGGGCAATCACACCCACAAACAAAAACAGCAGCACCAGTTCAg    >  5:159649/1‑148 (MQ=255)
                                                                                                                                                  accacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAg    >  4:166052/1‑147 (MQ=255)
                                                                                                                                                  accacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGAACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAg    >  6:443398/1‑147 (MQ=255)
                                                                                                                                                  accacTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCAGCACGGCCCGCCAGACGGCGGTCGACGAAGGAGAGCCAAAGGGCAATCACACCCACATACAAAAACAGCAGCACCACTTCAGaa  >  7:83695/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CATCAGCAAAATCAAATTACGCATCCGCATCAGCGTCAGCGTGGCAGACAGGAGCAGGAAGCCAAAAATCGACCAGCGCAGCGGCGAACGCATCGCCGGATCGATACTTTGCTGCATCCGCGTTGATCCCTGATGCAGGGTGTTCCACCACTCCACGGAGTAATGAATAATCGGCAGATTCACCACGCCAATCAGCACCAGGATACCTGCCGCACGGCCCGCCAGACGGCGGTCGTCGAAGGCGTGCCACAGGGCAATCACACCCACATACAAAAACAGCAGCACCAGTTCAGAA  >  NC_000913/2295637‑2295931

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