Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,248,812 T→C 5.2% N134D (AAT→GAT)  U712_11310 ← Riboflavin biosynthesis protein RibD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,248,8120TC 5.2% 238.1 / 3.4 77N134D (AAT→GAT) U712_11310Riboflavin biosynthesis protein RibD
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (2/2);  total (39/38)
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.15e-01

AGCATCCTGTCTTGCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAACACTCTTTTGATACCAGAGTTG  >  CP006881/2248675‑2248948
                                                                                                                                         |                                                                                                                                        
aGCATCCTGTCTTGCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCAGTCCGCATAAAGTGCAGAAATTTttcatt                                                                                                                                          >  1:306423/1‑138 (MQ=255)
  cATCCTGTCTTGCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAg                                                                                                                                       >  8:71485/1‑139 (MQ=255)
     ccTGTCTTGCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCct                                                                                                                                    >  6:160977/1‑139 (MQ=255)
       tGTCTTGCAGCCTCTGACGTGATCCATTTGATGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGAACGGAAGGCCTGTCCTCATAAACTGCAGAAATTTTTCAATCAg                                                                                                                                       >  2:262385/1‑134 (MQ=255)
             gCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGc                                                                                                                              <  4:143080/137‑1 (MQ=255)
             gCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGc                                                                                                                              <  4:138296/137‑1 (MQ=255)
             gCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCt                                                                                                                            >  6:255175/1‑139 (MQ=255)
                gCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGt                                                                                                                         <  2:217222/139‑1 (MQ=255)
                 ccTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCGTTTTGCGTGTCGTACGGAAGGCCTGTCCTCATAAAGGGCAGAAATTTTTCCTTCCGCCTCTCCGCCGGGTc                                                                                                                        >  2:253807/1‑139 (MQ=255)
                   tctGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTg                                                                                                                      >  4:200905/1‑139 (MQ=255)
                   tctGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTg                                                                                                                      >  5:26567/1‑139 (MQ=255)
                       aCGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGcc                                                                                                                    <  4:91260/137‑1 (MQ=255)
                        cGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAgg                                                                                                                 <  1:226507/139‑1 (MQ=255)
                           gATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATg                                                                                                              >  6:41969/1‑139 (MQ=255)
                            aTCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGc                                                                                                             >  8:82154/1‑139 (MQ=255)
                            aTCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGc                                                                                                             >  8:274031/1‑139 (MQ=255)
                             tCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGcc                                                                                                            <  7:82154/139‑1 (MQ=255)
                              ccATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           <  8:215835/139‑1 (MQ=255)
                               cATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           <  7:274031/138‑1 (MQ=255)
                               cATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           <  1:253807/138‑1 (MQ=255)
                               cATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           <  3:200905/138‑1 (MQ=255)
                               cATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTAAGGAAGGCCAGTCCACATAAAGTGCAGAAATATTACATTAAGCCTCTCCGCCTGGTCAGCCAGGATGCCtt                                                                                                          >  2:226507/1‑139 (MQ=255)
                                 tttGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGcc                                                                                                                             >  1:265331/1‑118 (MQ=255)
                                 tttGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGcc                                                                                                                             <  2:265331/118‑1 (MQ=255)
                                    gCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCt                                                                                                      <  6:168105/138‑1 (MQ=255)
                                    gCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCt                                                                                                      >  8:69012/1‑138 (MQ=255)
                                        tCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCt                                                                                                 <  2:306423/139‑1 (MQ=255)
                                               tGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGcc                                                                                          >  3:153858/1‑139 (MQ=255)
                                                    gTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGctt                                                                                     >  5:244520/1‑139 (MQ=255)
                                                     tAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGcttc                                                                                    <  8:34168/139‑1 (MQ=255)
                                                           tCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTtca                                                                               <  3:310832/138‑1 (MQ=255)
                                                            cTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGCCGTACGGAAGGCCGGTCCTCATAAAGGGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGCTGCCTTCCCGTACCGCAATGCCAGCTTCTTtcagc                                                                             >  2:269170/1‑137 (MQ=255)
                                                               gCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTtcat                                                                              <  6:104655/135‑1 (MQ=255)
                                                                ccGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGc                                                                         >  6:285115/1‑139 (MQ=255)
                                                                   tCAAGGCTGGCAGCCGCTTTTAGCGTGCCCTACGGAAGGCATGTCCTCAAAAAATGCAGAAATTTTTAATTCAGCCTCTCCGCCTGGTCTGCCAGGGTGCCGTCCCATGCCGCAATGCCCGCTTCGTACATCATGCTgc                                                                      >  4:246204/1‑138 (MQ=255)
                                                                     aaGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCTCCTGGTCTGCCAGGATGc                                                                                                             <  8:18487/98‑1 (MQ=255)
                                                                     aaGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGc                                                                                                             >  7:18487/1‑98 (MQ=255)
                                                                      aGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATcc                                                                   <  1:149303/139‑1 (MQ=255)
                                                                            gCAGCCGCTTTTATCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTAATTAAGCCGCTCCGCCGGGGGTGCCAGGCAGGCGTCCCTTACCTCAATGCCCGCGTCGTTGATAAAGGAGATCCctcgtg                                                             >  2:149303/1‑136 (MQ=255)
                                                                            gCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTc                                                             <  7:69012/139‑1 (MQ=255)
                                                                                      ttAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCcaca                                                     <  3:114562/137‑1 (MQ=255)
                                                                                      ttAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACaa                                                    <  5:289390/138‑1 (MQ=255)
                                                                                       tAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAg                                                   <  4:257535/138‑1 (MQ=255)
                                                                                       tAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAg                                                   <  6:26567/138‑1 (MQ=255)
                                                                                       tAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAg                                                   <  5:41969/138‑1 (MQ=255)
                                                                                       tAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCCGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAg                                                   <  5:305686/138‑1 (MQ=255)
                                                                                       tAGCGTGACGTACGGAAGGCCTGTACTCATAAAGTGCAGAATTTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAg                                                   <  1:262385/138‑1 (MQ=255)
                                                                                             gccGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATAATGCTGATCCCGCTTCCCGCAACACGCGGAttt                                            >  8:72541/3‑138 (MQ=255)
                                                                                              aCGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTCCCTCAATGCCAGCTTCTTTCCTCATGCTCATCCCTCTGCCAGCACCAACCCGATTAg                                           >  2:135552/1‑139 (MQ=255)
                                                                                              aCGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATCCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAg                                                                                        <  4:179564/94‑1 (MQ=255)
                                                                                              aCGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATCCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAg                                                                                        >  3:175494/1‑94 (MQ=255)
                                                                                              aCGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATCCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAg                                                                                        <  4:175494/94‑1 (MQ=255)
                                                                                              aCGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATCCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAg                                                                                        >  3:179564/1‑94 (MQ=255)
                                                                                                 tACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGAt                                        <  1:135552/139‑1 (MQ=255)
                                                                                                        ggCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCGTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCAt                                 >  8:129782/1‑139 (MQ=255)
                                                                                                           cTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGc                              <  4:153858/139‑1 (MQ=255)
                                                                                                            tGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATc                                                                    >  1:237314/1‑100 (MQ=255)
                                                                                                            tGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATc                                                                    <  2:237314/100‑1 (MQ=255)
                                                                                                            tGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGc                              >  5:155214/1‑138 (MQ=255)
                                                                                                              tCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGAt                                                                     >  4:292495/1‑97 (MQ=255)
                                                                                                              tCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGAt                                                                     <  3:292495/97‑1 (MQ=255)
                                                                                                              tCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCAc                           >  1:24349/1‑139 (MQ=255)
                                                                                                              tCCTCATAAAGTGCAGAAATTTTTAATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCCGATTCTTTCATCATGCAGATCCCTCTTCCAGCCACAAGAGGATTAGGATCTCTCATCTCCAc                           >  7:294807/1‑139 (MQ=255)
                                                                                                                cTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTATCATCGCCACg                          >  2:73509/1‑138 (MQ=255)
                                                                                                                cTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTGCTTTCTTCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACg                          >  8:70337/1‑138 (MQ=255)
                                                                                                                   aTAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           <  1:172572/54‑1 (MQ=255)
                                                                                                                   aTAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           >  2:172572/1‑54 (MQ=255)
                                                                                                                   aTAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGaa                        <  6:259806/137‑1 (MQ=255)
                                                                                                                   aTAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAaca                      >  5:259806/1‑139 (MQ=255)
                                                                                                                        gTGCAGAAATATTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGACTTCCCTTACATCAATGACAGCTTCTTTCATCATTCTGATCCCTCTACCAGCCACAAGCGGATTAGGATCTCTCAAAGCCACGAACActcta                 >  4:203595/1‑138 (MQ=255)
                                                                                                                         tGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAACACTCttt                <  8:26266/139‑1 (MQ=255)
                                                                                                                           cAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTtcat                                                                              <  3:263604/75‑1 (MQ=255)
                                                                                                                           cAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTtcat                                                                              >  4:263604/1‑75 (MQ=255)
                                                                                                                           cAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAacac                     >  8:227509/1‑132 (MQ=255)
                                                                                                                           cAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAacac                     <  7:227509/132‑1 (MQ=255)
                                                                                                                           cAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAGCACTCtttt               >  4:87045/1‑138 (MQ=255)
                                                                                                                                   tttCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           <  1:52609/38‑1 (MQ=255)
                                                                                                                                   tttCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCt                                                                                                           >  2:52609/1‑38 (MQ=255)
                                                                                                                                       attcaGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAACACTCTTTTGATACCAGAGTTg  <  5:160977/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
AGCATCCTGTCTTGCAGCCTCTGACGTGATCCATTTGCTGTCACCCGTGCTGGTAGCTATCTTGCCGTCAAGGCTGGCAGCCGCTTTTAGCGTGACGTACGGAAGGCCTGTCCTCATAAAGTGCAGAAATTTTTCATTCAGCCTCTCCGCCTGGTCTGCCAGGATGCCTTCCCTTACCTCAATGCCAGCTTCTTTCATCATGCTGATCCCTCTTCCAGCCACAAGCGGATTAGGATCTCTCATCGCCACGAACACTCTTTTGATACCAGAGTTG  >  CP006881/2248675‑2248948

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: