Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,994,519 G→A 40.3% intergenic (‑103/‑267) yigE ← / → corA DUF2233 family protein/magnesium/nickel/cobalt transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,994,5190GA40.3% 35.4 / 61.4 57intergenic (‑103/‑267)yigE/corADUF2233 family protein/magnesium/nickel/cobalt transporter
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base A (11/12);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCA  >  CP009273/3994280‑3994727
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                
tGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTGCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACg                                                                                                                                                                               >  1:204506/1‑275 (MQ=255)
         aGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGCTTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCCTTATAGGTCCTGGATGTGGGATTTTTTTATCCTTTTAGCGACCTTTACGAGTACCCAAAAGCGCGCAGTTCCACTCTTGTTCTGTGTTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTACCCTGCGCGCCCGATCActc                                                                                                                                                                        >  2:167597/1‑273 (MQ=255)
                          tGCGAAATCATCAGCGGCAAAAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCt                                                                                                                                                       <  1:167600/273‑1 (MQ=255)
                           gcgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTCTCGTGTTGACGGAAAAATTTTGCTGGCGTAACATgcg                                                                                                                                                                                     >  2:618583/1‑242 (MQ=255)
                            cgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAgg                                                                                                                                                                 >  2:233609/1‑261 (MQ=255)
                            cgcAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCa                                                                                                                                                                                 >  1:542096/1‑245 (MQ=255)
                              cAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGcc                                                                                                                                                        <  1:391775/268‑1 (MQ=255)
                               aaTCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGTAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGTGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTTTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCa                                                                                                                                                                                 <  1:618594/242‑1 (MQ=255)
                                atcatcAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACAt                                                                                                                                                             <  1:233613/261‑1 (MQ=255)
                                atcatcAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGAt                                                                                                                                                                             <  2:542085/245‑1 (MQ=255)
                                                                           agaAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCa                                                                                                                                                                                 >  2:182540/1‑198 (MQ=255)
                                                                               aaaTCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGAt                                                                                                                                                                             <  1:182544/198‑1 (MQ=255)
                                                                                                      aTCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAAt                                                                                                                                                                                                        >  1:300457/1‑148 (MQ=255)
                                                                                                      aTCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGACCAATATCCTGTCATCCAGGACCAACCCCTAAGGCTATCCCTTTTTTCTgtt                                                                           >  2:480126/1‑271 (MQ=255)
                                                                                                          ttCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTg                                                                                                                                                                                                    <  2:300452/148‑1 (MQ=255)
                                                                                                                       gAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTc                                                                                                                             >  1:126455/1‑206 (MQ=255)
                                                                                                                           aGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTa                                                                                                                         <  2:126454/206‑1 (MQ=255)
                                                                                                                                aTGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAAGATGCGCGCACGATAACTCTAAGAGGACAt                                                                                                                                                             <  1:482120/165‑1 (MQ=255)
                                                                                                                                aTGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTAGCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACAt                                                                                                                                                             <  1:24669/165‑1 (MQ=255)
                                                                                                                                        ttATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCt                                                                                                                              >  1:340981/1‑188 (MQ=255)
                                                                                                                                            aGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACt                                                                                                                          <  2:340976/188‑1 (MQ=255)
                                                                                                                                              gTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTCTTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACg                                                                                                                                                                               >  2:367720/1‑133 (MQ=255)
                                                                                                                                                  tGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCa                                                                                                                                                                           <  1:367725/133‑1 (MQ=255)
                                                                                                                                                                          ttAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCtt                                                                                    >  1:45400/1‑196 (MQ=255)
                                                                                                                                                                             gCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACNTGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCTCGCTGTTTAATGGTCTGAAAcaca  <  1:480134/275‑1 (MQ=255)
                                                                                                                                                                              cGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTg                                                                                <  2:45400/196‑1 (MQ=255)
                                                                                                                                                                                 ccTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACa                                                                                                                                                                                         >  2:610101/1‑88 (MQ=255)
                                                                                                                                                                                    tGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTc                                                                                                                                                           >  2:361762/1‑115 (MQ=255)
                                                                                                                                                                                     gACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATgcg                                                                                                                                                                                     <  1:610112/88‑1 (MQ=255)
                                                                                                                                                                                        gAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCt                                                                                                                                                       <  1:361767/115‑1 (MQ=255)
                                                                                                                                                                                           tACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCtt                                                                                                                                                      >  1:302926/1‑113 (MQ=255)
                                                                                                                                                                                           tACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGa                              >  2:631078/1‑233 (MQ=255)
                                                                                                                                                                                               aaaaaGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGac                                                                                                                                                  <  2:302921/113‑1 (MQ=255)
                                                                                                                                                                                               aaaaaGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCa                          <  1:631089/233‑1 (MQ=255)
                                                                                                                                                                                                aaaaGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACa                                                                                                                                                                                         >  1:98377/1‑73 (MQ=255)
                                                                                                                                                                                                    gcgcGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATgcg                                                                                                                                                                                     <  2:98377/73‑1 (MQ=255)
                                                                                                                                                                                                       agAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATg                                                                                                                                                                                       >  2:214830/2‑68 (MQ=255)
                                                                                                                                                                                                           gTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATgcgcg                                                                                                                                                                                   <  1:214834/68‑1 (MQ=255)
                                                                                                                                                                                                                      gTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGa                                                                                                                                                                              >  2:508367/1‑62 (MQ=255)
                                                                                                                                                                                                                         cTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAgg                                                                                                                                                                 >  2:98929/1‑72 (MQ=255)
                                                                                                                                                                                                                          tgtgGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCAc                                                                                                                                                                          <  1:508377/62‑1 (MQ=255)
                                                                                                                                                                                                                            tggtgTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTaa                                                                                                                                                                                           >  2:75423/1‑43 (MQ=255)
                                                                                                                                                                                                                            tggtgTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTgg                                                                                                                                >  2:120724/1‑102 (MQ=255)
                                                                                                                                                                                                                            tggtgTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTAt                                                                                                                        >  2:514657/1‑110 (MQ=255)
                                                                                                                                                                                                                            tggtgTTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCtt                                                                                                                                                      >  2:531121/1‑80 (MQ=255)
                                                                                                                                                                                                                             ggtgTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACAt                                                                                                                                                             <  1:98929/72‑1 (MQ=255)
                                                                                                                                                                                                                                gtTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATg                                                                                                                                                                                       <  1:75423/43‑1 (MQ=255)
                                                                                                                                                                                                                                gtTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCa                                                                                                                            <  1:120725/102‑1 (MQ=255)
                                                                                                                                                                                                                                gtTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTg                                                                                                                    <  1:514667/110‑1 (MQ=255)
                                                                                                                                                                                                                                gtTCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGac                                                                                                                                                  <  1:531132/80‑1 (MQ=255)
                                                                                                                                                                                                                                 ttCTGTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTa                                                                                                                                                                      >  1:128043/1‑59 (MQ=255)
                                                                                                                                                                                                                                     gTTGCGTGTTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAagag                                                                                                                                                                  <  2:128042/59‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCa                                                                                                                                                                                 >  1:332617/1‑43 (MQ=255)
                                                                                                                                                                                                                                          gtgtTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGAt                                                                                                                                                                             <  2:332612/43‑1 (MQ=255)
                                                                                                                                                                                                                                          gtgtTGACGGCAAAATTTTGCTCGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGACTTAGCGGTTGTCAGCGATCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTaa                <  1:355877/200‑1 (MQ=255)
                                                                                                                                                                                                                                           tgtTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGaca                                                                                                                                                 >  1:526697/1‑70 (MQ=255)
                                                                                                                                                                                                                                           tgtTAACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAg                                                                                                            >  1:444171/1‑107 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                
TGTACGGTCAGCGTCGGATCTGAGAGTGCGCAATCATCAGCGGCAACAGCAAACAAGGGAAGTAAGGTGAGGGCGAGAAAAATCCGTTTGAGATTCAAGGTGATCATTCCTTTACCAATGAGTAGCTGATGCGCCATTATAGGTCCTGGATGTGGGATTTTTTTATCCTGTTAGCGACCTTGACGAGTACCAAAAAGCGCGAAGTTCAACTATTGTTCTGTGGTGTTCTGTTGCGTGTTGACGGCAAAATTTTGCTGGCGTAACATGCGCGCACGATCACTCTAAGAGGACATTCGCCTTGGACACACCCAGTAGATACTGGCTCACTATCCTGTCATCCAGGATCAACTCCTAAGGCTATCCCTTTTTGCTGATAGCCTTAGCGGTTGTCAGCGACCTCAATTTTTCCCGTCGCGCTGAGTCAGGCTGTTTAATGGTCTGAAACCCA  >  CP009273/3994280‑3994727

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