Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,641,937 T→C 5.6% S50G (AGC→GGC)  U712_18775 ← Nitrate reductase gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,641,9370TC 5.6% 213.4 / 3.0 71S50G (AGC→GGC) U712_18775Nitrate reductase gamma chain
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (2/2);  total (35/36)
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.21e-01

CAAGCCAGCGCCAATCGCCATTTTGTGATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATACCCCAGAGGAT  >  CP006881/3641803‑3642072
                                                                                                                                      |                                                                                                                                       
cAAGCCAGCGCCAATCGCCATTTTGTGATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAATAAGTTtgctgc                                                                                                                                       >  6:307844/1‑137 (MQ=255)
             aTCGCCATTTTGTGATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGttttt                                                                                                                        >  2:82823/1‑139 (MQ=255)
             aTCGCCATTTTGTGATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCACCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGttttt                                                                                                                        >  4:273244/1‑139 (MQ=255)
                       tgtgATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTaa                                                                                                              >  4:184104/1‑139 (MQ=255)
                         tgATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTaa                                                                                                              >  3:379047/1‑137 (MQ=255)
                              cATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCg                                                                                                       >  2:343706/1‑139 (MQ=255)
                                    ctctGTAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCTCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGa                                                                                                  <  6:247013/138‑1 (MQ=255)
                                      ctGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGAttt                                                                                               >  3:389184/1‑139 (MQ=255)
                                       tGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTc                                                                                              >  6:201964/1‑139 (MQ=255)
                                       tGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTc                                                                                              >  5:330944/1‑139 (MQ=255)
                                       tGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTc                                                                                              >  3:192258/1‑139 (MQ=255)
                                       tGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGATCGCTTGAtt                                                                                                >  3:61186/1‑137 (MQ=255)
                                         aaaTGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCCGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGc                                                                                            >  6:207447/1‑139 (MQ=255)
                                         aaaTGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCCGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGc                                                                                            >  4:287453/1‑139 (MQ=255)
                                           aTGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCg                                                                                          <  6:227937/139‑1 (MQ=255)
                                              ccAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTc                                                                                        <  4:85673/138‑1 (MQ=255)
                                                 aGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGc                                                                                    >  5:5765/1‑139 (MQ=255)
                                                       gCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACa                                                                                                            >  6:210343/1‑109 (MQ=255)
                                                       gCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACa                                                                                                            <  5:210343/109‑1 (MQ=255)
                                                       gCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCTGTCCAGCCAAAtt                                                                              >  6:187711/1‑139 (MQ=255)
                                                       gCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTGTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAAtt                                                                              >  1:98714/1‑139 (MQ=255)
                                                             acGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCAt                                                                        <  1:79680/139‑1 (MQ=255)
                                                                  ttCTGGGATCAGAATCCCCATTACATGCCAGCAGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGa                                                                   >  1:12564/1‑139 (MQ=255)
                                                                  ttCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGa                                                                   >  6:127328/1‑139 (MQ=255)
                                                                    cTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCAGAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGAt                                                                  <  5:307844/138‑1 (MQ=255)
                                                                    cTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCTCCGTCCAGCCAAATTGGTCATGCTGAt                                                                  <  3:51080/138‑1 (MQ=255)
                                                                    cTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGAt                                                                  <  4:389184/138‑1 (MQ=255)
                                                                    cTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCAAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGAt                                                                  <  6:330944/138‑1 (MQ=255)
                                                                       ggATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCCGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCg                                                              <  5:207447/139‑1 (MQ=255)
                                                                       ggATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCCGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCg                                                              <  3:287453/139‑1 (MQ=255)
                                                                         aTCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCgg                                                             >  7:327215/1‑138 (MQ=255)
                                                                           cAGAATCCCCATGACATTCCCGACGACAACGCACAGCAGTCCCAACTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTATTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGAGGTa                                                           >  3:271495/1‑138 (MQ=255)
                                                                           cAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTa                                                           >  1:391365/1‑138 (MQ=255)
                                                                           cAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTa                                                           >  1:84718/1‑138 (MQ=255)
                                                                                  cccATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGgc                                                   <  6:5765/139‑1 (MQ=255)
                                                                                   ccATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGAt                                                                  <  7:85531/123‑1 (MQ=255)
                                                                                   ccATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGAt                                                                  >  8:85531/1‑123 (MQ=255)
                                                                                    cATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGGCCAGCCATATTGGTCATGCTGATATCGGTAAATATGgccg                                                 <  4:61186/139‑1 (MQ=255)
                                                                                      tGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGgccgcc                                               >  6:215549/1‑139 (MQ=255)
                                                                                        aCATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGccgccg                                              >  7:257369/1‑138 (MQ=255)
                                                                                          aTGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGAt                                            >  7:320158/1‑138 (MQ=255)
                                                                                            gCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCACATTGGTCAGGCTGATAGCGGTAAATATGGCCGCCGATa                                           >  7:254486/1‑137 (MQ=255)
                                                                                            gCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATa                                           >  6:112192/1‑137 (MQ=255)
                                                                                            gCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATa                                           >  7:85468/1‑137 (MQ=255)
                                                                                             cccGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAg                                        <  2:84718/139‑1 (MQ=255)
                                                                                                  cgACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGt                                   <  1:82823/139‑1 (MQ=255)
                                                                                                  cgACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGt                                   <  3:273244/139‑1 (MQ=255)
                                                                                                  cgACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGt                                   <  8:257369/139‑1 (MQ=255)
                                                                                                  cgACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAAAATGTCCGCCGATAAAGATTGt                                   <  2:70646/139‑1 (MQ=255)
                                                                                                    aCAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCa                                 <  2:391365/139‑1 (MQ=255)
                                                                                                        cGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATa                              >  2:391774/1‑138 (MQ=255)
                                                                                                         gCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            <  3:350335/139‑1 (MQ=255)
                                                                                                         gCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            <  4:74957/139‑1 (MQ=255)
                                                                                                         gCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            <  8:320158/139‑1 (MQ=255)
                                                                                                            cagcagTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTGTTCTAACCGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATg                         >  6:289264/1‑139 (MQ=255)
                                                                                                              gcagTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCATCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       <  3:147590/139‑1 (MQ=255)
                                                                                                              gcagTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       <  5:46392/139‑1 (MQ=255)
                                                                                                              gcagTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       <  6:64867/139‑1 (MQ=255)
                                                                                                              gcagTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       <  3:184104/139‑1 (MQ=255)
                                                                                                              gcagTCCCCAGTGAAAAAGTGTGCTGACAGCCTCAAGTGTTTTCTTGTCTAACAGCTCGCTTGATTTCGCCGTCCTGCCAAATTGGTCATGCTGTTAGCGGTATATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       <  5:185372/139‑1 (MQ=255)
                                                                                                                  tCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATgg                    >  2:8922/1‑138 (MQ=255)
                                                                                                                    cccAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       <  6:50187/133‑1 (MQ=255)
                                                                                                                    cccAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTa                       >  5:50187/1‑133 (MQ=255)
                                                                                                                      cAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCataa               <  5:215549/139‑1 (MQ=255)
                                                                                                                      cAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCGGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCataa               <  4:243057/139‑1 (MQ=255)
                                                                                                                       aGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGt                                   >  5:230735/1‑118 (MQ=255)
                                                                                                                       aGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGt                                   <  6:230735/118‑1 (MQ=255)
                                                                                                                       aGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            >  4:360091/1‑125 (MQ=255)
                                                                                                                       aGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            <  3:360091/125‑1 (MQ=255)
                                                                                                                           aaaaaGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATAccc          <  1:343706/139‑1 (MQ=255)
                                                                                                                            aaaaGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCTGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATAccc          <  5:187711/138‑1 (MQ=255)
                                                                                                                                gtgtGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATACCCCagag     <  4:115270/139‑1 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGgc                                                   >  5:249489/1‑90 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGgc                                                   <  6:249489/90‑1 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCataa               >  1:325742/1‑126 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCataa               <  4:163388/126‑1 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCataa               >  3:163388/1‑126 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCataa               <  2:325742/126‑1 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATACCCCa        >  6:765/1‑133 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGAATGTCAATACAATGTATGGCATAATACCCCa        <  5:765/133‑1 (MQ=255)
                                                                                                                                   tgctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGAATTCGCAGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGTCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATACCCCAGAGGAt  >  4:242815/1‑139 (MQ=255)
                                                                                                                                     ctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            <  5:105095/111‑1 (MQ=255)
                                                                                                                                     ctgcCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTAATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACa                            >  6:105095/1‑111 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CAAGCCAGCGCCAATCGCCATTTTGTGATACATATGCTCTGAAATGCCAAGGGAAGCATACACGCCTTCTGGGATCAGAATCCCCATGACATGCCCGCCGACAACGCACAGCAGTCCCCAGTGAAAAAGTGTGCTGCCAGCCGCAAGTTTTTTCTTTTCTAACAGCTCGCTTGATTTCGCCGTCCAGCCAAATTGGTCATGCTGATAGCGGTAAATATGGCCGCCGATAAAGATTGTCAATACAATGTATGGCATAATACCCCAGAGGAT  >  CP006881/3641803‑3642072

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: ‑