Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,519,942 C→A 6.7% W230C (TGG→TGT narG ← nitrate reductase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,519,9420CA 6.7% 255.3 / 5.0 75W230C (TGG→TGTnarGnitrate reductase alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base A (3/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

ATGACTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTG  >  CP000730/2519802‑2520078
                                                                                                                                            |                                                                                                                                        
aTGGCTTTTGTACCTTGAGAGCGGTCTTCAGTCATAAAATGTGCAGCCGGAGTACGTGTTGAAGGTACATTAGAGCCCCACCTAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTccc                                                                                                                                          <  2:527205/141‑1 (MQ=255)
aTGACTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTccc                                                                                                                                          <  1:23985/141‑1 (MQ=255)
  gACTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCa                                                                                                                                        <  2:661202/141‑1 (MQ=255)
   aCTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCaa                                                                                                                                       <  2:439744/141‑1 (MQ=255)
   aCTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCaa                                                                                                                                       <  2:326455/141‑1 (MQ=255)
    cTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCaaa                                                                                                                                      <  2:265328/141‑1 (MQ=255)
    cTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCaaa                                                                                                                                      >  2:539270/1‑141 (MQ=255)
    cTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACAGTAGAGCCCCACATAATAATGTAGGATGCGTTATACCAGTCACTTGATGCAGGCACATCTGTTTGCTCTCCCCaaa                                                                                                                                      >  1:469969/1‑141 (MQ=255)
     ttttGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAAt                                                                                                                                     >  2:479938/1‑141 (MQ=255)
       ttGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAAttt                                                                                                                                   >  1:104766/1‑141 (MQ=255)
            ccTTTATAGCGGCCTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCACAAATTTGGGga                                                                                                                              <  2:289768/141‑1 (MQ=255)
            ccTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGga                                                                                                                              <  1:750450/141‑1 (MQ=255)
               ttATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGagag                                                                                                                           <  2:257923/141‑1 (MQ=255)
                     cGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggt                                                                                                                     <  2:107343/141‑1 (MQ=255)
                      ggACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAggtg                                                                                                                    >  1:79620/1‑141 (MQ=255)
                          ttCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaa                                                                                                                >  2:585572/1‑141 (MQ=255)
                           tCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTaaa                                                                                                               >  1:381002/1‑141 (MQ=255)
                              gTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCt                                                                                                            >  1:395803/1‑141 (MQ=255)
                               tCATAAAATGTGCATCCGGAGTACGTGTTAAGGTTACATTATAGGCCCAGAGGATAATGTATTATTCTGTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTg                                                                                                           <  1:662278/141‑1 (MQ=255)
                                cATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGc                                                                                                          <  1:323234/141‑1 (MQ=255)
                                  aaaaaTGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCAt                                                                                                        <  2:433489/140‑1 (MQ=255)
                                      aTGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCa                                                                                                    >  2:384073/1‑141 (MQ=255)
                                       tgtgCATCCGGAGTACGGGTTGAAGGTACTTTAGAGCCCCACATTATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAg                                                                                                   <  2:231532/141‑1 (MQ=255)
                                        gtgCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGt                                                                                                  >  2:301581/1‑141 (MQ=255)
                                         tgCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTc                                                                                                 <  1:719967/141‑1 (MQ=255)
                                           cATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                                                               >  1:596133/1‑141 (MQ=255)
                                           cATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                                                               <  2:325559/141‑1 (MQ=255)
                                              ccGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCaaa                                                                                                                                      >  2:523007/1‑99 (MQ=255)
                                              ccGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCaaa                                                                                                                                      <  1:523007/99‑1 (MQ=255)
                                                      cGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGc                                                                                    >  2:62403/1‑141 (MQ=255)
                                                            aaaGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTca                                                                              >  1:209140/1‑141 (MQ=255)
                                                             aaGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGCAGTTCCACATAGAACTGATAAAAACTATGCAGTTcgc                                                                             >  1:424605/1‑139 (MQ=255)
                                                                ggaCATTAGAACCGCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTCAATCTGCCTACCAGTCATAAAAACTaaaaacctgacatc                                                                          >  2:28984/3‑129 (MQ=255)
                                                                 tACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCaccacc                                                                         <  1:467947/141‑1 (MQ=255)
                                                                   cATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCaa                                                                       <  2:478075/141‑1 (MQ=255)
                                                                   cATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCaa                                                                       >  2:640663/1‑141 (MQ=255)
                                                                    aTTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAg                                                                      <  1:539270/141‑1 (MQ=255)
                                                                      tAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCa                                                                    <  1:61658/141‑1 (MQ=255)
                                                                         agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGga                                                                                                                              >  2:56931/1‑80 (MQ=255)
                                                                         agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGga                                                                                                                              <  1:56931/80‑1 (MQ=255)
                                                                         agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTCATTCAAGCACAACTGTTTGCTCTCCCCACATCTGTGGAGAGGCCGGTGGTAAATCTGCATACAAGTCATAAAGACTAAGCCTTTCACCACCAAGCAAAt                                                                 >  2:503097/1‑141 (MQ=255)
                                                                         agCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTCAGTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAAt                                                                 >  1:574191/1‑141 (MQ=255)
                                                                            cccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtga                                                              <  2:54739/141‑1 (MQ=255)
                                                                             ccACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtgat                                                             <  1:581940/141‑1 (MQ=255)
                                                                                   aataatGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCg                                                       <  2:602408/141‑1 (MQ=255)
                                                                                        tGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCAc                                                  <  2:487566/141‑1 (MQ=255)
                                                                                          tATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACct                                                >  2:135048/1‑141 (MQ=255)
                                                                                           atgatgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctg                                               <  2:730504/141‑1 (MQ=255)
                                                                                              atgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAATACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg                                            >  1:336103/1‑141 (MQ=255)
                                                                                              atgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACctgctg                                            >  1:451782/1‑141 (MQ=255)
                                                                                               tgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc                                           <  2:745443/141‑1 (MQ=255)
                                                                                               tgCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCtgctgc                                           >  1:431360/1‑141 (MQ=255)
                                                                                                gCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCa                                          >  1:720089/1‑141 (MQ=255)
                                                                                                   ttATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGGGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATaa                                       >  2:283308/1‑141 (MQ=255)
                                                                                                   ttATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGCGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATaa                                       <  1:589535/141‑1 (MQ=255)
                                                                                                   ttATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATaa                                       <  1:229628/141‑1 (MQ=255)
                                                                                                     ataCCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACt                                     >  1:456601/1‑141 (MQ=255)
                                                                                                           gTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCAt                               >  2:680803/1‑141 (MQ=255)
                                                                                                             cACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTg                             <  1:538820/141‑1 (MQ=255)
                                                                                                              aCTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGa                            <  2:765008/141‑1 (MQ=255)
                                                                                                                 tGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACAt                         <  1:135048/141‑1 (MQ=255)
                                                                                                                  gATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTATGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATc                        >  2:416763/1‑141 (MQ=255)
                                                                                                                     tCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCt                     <  2:211927/141‑1 (MQ=255)
                                                                                                                     tCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCt                     >  2:113439/1‑141 (MQ=255)
                                                                                                                        gatttaccacctgcctctCCACAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGa                  >  2:84764/15‑141 (MQ=38)
                                                                                                                            caTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAg              >  1:325968/1‑141 (MQ=255)
                                                                                                                             aTCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAgg             >  2:46841/1‑141 (MQ=255)
                                                                                                                               cTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGgtg           <  2:230220/141‑1 (MQ=255)
                                                                                                                               cTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGgtg           <  2:706482/141‑1 (MQ=255)
                                                                                                                                tGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGgtgt          >  2:352589/1‑141 (MQ=255)
                                                                                                                                  tttGCGCTCCCCAAAGTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTaa        <  2:224087/141‑1 (MQ=255)
                                                                                                                                   ttGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCAGGatttacacctat       <  1:109261/141‑12 (MQ=255)
                                                                                                                                      ctctCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATcc    >  1:241099/1‑141 (MQ=255)
                                                                                                                                       tctcCACAAATTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGCATAGGGGTAAATCCt   >  2:71848/1‑141 (MQ=255)
                                                                                                                                        ctcCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTg  <  2:622624/141‑1 (MQ=255)
                                                                                                                                        ctcCACAAGTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtgat                                                             >  1:307945/1‑82 (MQ=255)
                                                                                                                                        ctcCACAAGTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATtgat                                                             <  2:307945/82‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
ATGACTTTTGTACCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTG  >  CP000730/2519802‑2520078

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