Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,234,371 C→T 6.4% G236S (GGC→AGC)  U712_16665 ← Iron(3+)‑hydroxamate import system permease protein fhuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,234,3710CT 6.4% 199.9 / 2.2 62G236S (GGC→AGC) U712_16665Iron(3+)‑hydroxamate import system permease protein fhuG
Reads supporting (aligned to +/- strand):  ref base C (22/36);  new base T (2/2);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

GGTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGCCCCCCATATTGAGCCGG  >  CP006881/3234240‑3234508
                                                                                                                                   |                                                                                                                                         
ggTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCg                                                                                                                                    >  7:289211/1‑139 (MQ=255)
 gTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCgg                                                                                                                                   >  3:326021/1‑139 (MQ=255)
 gTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCgg                                                                                                                                   >  6:220647/1‑139 (MQ=255)
 gTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCgg                                                                                                                                   >  5:207848/1‑139 (MQ=255)
 gTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCAGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCgg                                                                                                                                   >  8:210706/1‑139 (MQ=255)
            aGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTgg                                                                                                                        >  6:69019/1‑139 (MQ=255)
             gCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAAtt                                                                                                                          <  8:58431/136‑1 (MQ=255)
             gCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAAtt                                                                                                                          <  4:111843/136‑1 (MQ=255)
               cacaTGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGtct                                                                                                                     <  7:210706/139‑1 (MQ=255)
                acaTGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCTCAGCCCGGTCGCCAATTGGtctc                                                                                                                    >  6:44835/1‑139 (MQ=255)
                acaTGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGtctc                                                                                                                    >  8:293450/1‑139 (MQ=255)
                  aTGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCTCAGCCCGGTCGCCAATTGGTCTCCg                                                                                                                  >  3:82713/1‑139 (MQ=255)
                    gcggcgCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGGCTCCGAg                                                                                                                <  5:69019/139‑1 (MQ=255)
                     cggcgCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAg                                                                                                                <  3:93498/138‑1 (MQ=255)
                      ggcgCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCt                                                                                                              <  6:111483/139‑1 (MQ=255)
                      ggcgCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGGGCCCAGCCCGGGCGCCAATTGGTCTCCGAGCt                                                                                                              <  1:256702/139‑1 (MQ=255)
                          cgATCAATCCGAGAAACGCAATGCCCCCCGCAGCTTCCACGCTTGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCAt                                                                                                          <  2:87661/139‑1 (MQ=255)
                          cgATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCAt                                                                                                          <  4:19335/139‑1 (MQ=255)
                          cgATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAAACCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCAt                                                                                                          <  1:251230/139‑1 (MQ=255)
                            aTCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATaa                                                                                                        <  6:140617/139‑1 (MQ=255)
                            aTCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATaa                                                                                                        >  5:140617/1‑139 (MQ=255)
                            aTCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATaa                                                                                                        >  5:108863/1‑139 (MQ=255)
                            aTCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGGCGCCAATTGGTCTCCGAGCTGCATaa                                                                                                        <  4:73751/139‑1 (MQ=255)
                                aTCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGtt                                                                                                    >  5:304364/1‑139 (MQ=255)
                                   cGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCaa                                                                                                 <  1:22268/139‑1 (MQ=255)
                                   cGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCa                                                                                                  <  6:68840/138‑1 (MQ=255)
                                        aaCGCAATGCCCCCCGCAGCTGCCACGCATGCAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGc                                                                                            <  7:293450/139‑1 (MQ=255)
                                        aaCGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGTTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGc                                                                                            <  8:289211/139‑1 (MQ=255)
                                        aaCGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGc                                                                                            <  6:108863/139‑1 (MQ=255)
                                        aaCGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGc                                                                                            <  5:220647/139‑1 (MQ=255)
                                         aCGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCg                                                                                           <  3:205559/139‑1 (MQ=255)
                                          cGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGt                                                                                          <  6:142148/139‑1 (MQ=255)
                                            cAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGc                                                                                        <  6:99763/139‑1 (MQ=255)
                                             aaggCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGcc                                                                                       <  6:207848/136‑1 (MQ=255)
                                             aaTGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGcc                                                                                       <  4:326021/139‑1 (MQ=255)
                                                 ccccccGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGATTCGGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTtat                                                                                   <  5:262571/139‑1 (MQ=255)
                                                       gCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTAAAATAGCTTGCCTTATTTAAAg                                                                             >  2:240047/1‑139 (MQ=255)
                                                                   cATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGc                                                                  >  5:118652/1‑138 (MQ=255)
                                                                    aTGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCaa                                                                <  5:125381/139‑1 (MQ=255)
                                                                        aGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTaaa                                                                         <  8:68040/126‑1 (MQ=255)
                                                                        aGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTaaa                                                                         >  7:68040/1‑126 (MQ=255)
                                                                            gCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAAt                                                                                                       >  6:11748/1‑92 (MQ=255)
                                                                            gCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAAt                                                                                                       <  5:11748/92‑1 (MQ=255)
                                                                                   tGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAgg                                                 <  5:75822/139‑1 (MQ=255)
                                                                                     aCCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTGTGCGCGCGGCATAATGTTCAAAAAGCGTGCCGTATATAAAGTAAACACAAGCAATATAACAATCaaaag                                                 >  3:295077/1‑135 (MQ=255)
                                                                                     aCCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAgg                                                 >  5:339270/1‑137 (MQ=255)
                                                                                     aCCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGGCGCCGCGCGGCAGAATGTGCCAATAGCGGGCCTGATATAAAGTAAACAGAAGCAAGCTCCCAATCCAAgggc                                               >  1:234062/1‑138 (MQ=255)
                                                                                          tGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATc                                          >  8:50238/1‑139 (MQ=255)
                                                                                             cAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCTCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGcc                                       <  4:82713/139‑1 (MQ=255)
                                                                                              agaagaagaATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGGCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGccc                                      <  7:231478/139‑1 (MQ=255)
                                                                                                        tCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATAtt                            >  8:36296/1‑139 (MQ=255)
                                                                                                          cGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTc                           >  7:1517/1‑138 (MQ=255)
                                                                                                             cTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGt                       <  6:339270/139‑1 (MQ=255)
                                                                                                                tCCTTTTCATCAGCGGTGCCCCGCCCGGTCGGCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTATCCTGAAGCAAGATCACAATCCAAGGGTGAATCGCCCAGATCATATTCCAGt                       <  4:295077/136‑1 (MQ=255)
                                                                                                                tCCTTTTCAACAGCGGTGCTCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTc                     <  5:44835/138‑1 (MQ=255)
                                                                                                                tCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTTGTCTCCGAGCTGCATAATGTTAAAATAGCGTGCCTTATATAAAGTAAACAGAACCAAGAACACAATACAAGGGAGATTGGCCCAGATCATATTCCAGttag                    >  5:325214/1‑137 (MQ=255)
                                                                                                                  cTTTTCCACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGc                   <  2:234062/138‑1 (MQ=255)
                                                                                                                  cTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGc                   <  6:304364/138‑1 (MQ=255)
                                                                                                                  cTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGc                   <  1:240047/138‑1 (MQ=255)
                                                                                                                    tttCAACAGCGGTGCCCAGCCCGGGCGCCCATTGGTCGCCGCGGGGGGTTATGTTGAAATATCGGGGCTGTTTAAAAGTTATAAAAAGCAAGAGCACAAAACAAGGGTGGATGGCCCACCTCAGTTTTCAGGTGGcccc                >  8:307830/1‑139 (MQ=255)
                                                                                                                    cagcAACAGCGGTGCCCAGCCCGGTCGCCAATTGGACTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGcccc                <  3:146343/136‑1 (MQ=255)
                                                                                                                     ttCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGcccc                <  8:1517/138‑1 (MQ=255)
                                                                                                                                  cccAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGCCCCCCATATTGAGCCgg  <  4:145025/139‑1 (MQ=255)
                                                                                                                                  cccAGCCCGGTCGCCAATTGGACTCCGCGCTGCATAATTGTCCAATAGCCTGCCGTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGCAGAATCGCCCAGAGCCTAATCCAGATCGCCCCCCATAATGATcctg  >  4:146343/1‑137 (MQ=255)
                                                                                                                                   |                                                                                                                                         
GGTCAGCCTTCTAGCCACATGCGGCGCGATCAATCCGAGAAACGCAATGCCCCCCGCAGCTGCCACGCATGAAGCGGCAAGAGTGACCGCTGCCAGAAGAAGAATCCGTCTTTCCTTTTCAACAGCGGTGCCCAGCCCGGTCGCCAATTGGTCTCCGAGCTGCATAATGTTCAAATAGCGTGCCTTATATAAAGTAAACAGAAGCAAGATCACAATCCAAGGGAGAATCGCCCAGATCATATTCCAGTTCGCCCCCCATATTGAGCCGG  >  CP006881/3234240‑3234508

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