Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,519,951 T→G 12.8% P227P (CCA→CCC narG ← nitrate reductase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,519,9510TG12.8% 166.7 / 24.3 71P227P (CCA→CCCnarGnitrate reductase alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base G (7/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTGG  >  CP000730/2519814‑2520091
                                                                                                                                         |                                                                                                                                            
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                                                                                                            <  2:325559/141‑1 (MQ=255)
                               cATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCAt                                                                                                            >  1:596133/1‑141 (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:640663/1‑141 (MQ=255)
                                                       cATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCaa                                                                                    <  2:478075/141‑1 (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:706482/141‑1 (MQ=255)
                                                                                                                   cTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGgtg                        <  2:230220/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)
                                                                                                                                   aaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTGATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCg        >  1:240974/1‑141 (MQ=255)
                                                                                                                                   aaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCg        >  1:641213/1‑141 (MQ=255)
                                                                                                                                   aaTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCg        >  2:252379/1‑141 (MQ=255)
                                                                                                                                    aTTTGGGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGa       >  1:710146/1‑141 (MQ=255)
                                                                                                                                      ttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATc     >  2:778944/1‑141 (MQ=255)
                                                                                                                                      ttGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATc     >  1:131858/1‑141 (MQ=255)
                                                                                                                                        gtgGATAGGCAGGTGGTATATCTGCATACCAGTCATAAAAACTGAGCATTTCACCACCAAGCAAATTGATGAATCGAGAACCTGCTGCATAACTGATCATTGACATCGCTGGAATAGGGGTAAATCCTGCGTTTCTTTCTg   >  2:574824/1‑141 (MQ=255)
                                                                                                                                        gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTg   >  2:88175/1‑141 (MQ=255)
                                                                                                                                        gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTg   >  1:673352/1‑141 (MQ=255)
                                                                                                                                        gtgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTg   >  2:651598/1‑141 (MQ=255)
                                                                                                                                         tgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTgg  >  2:462012/1‑141 (MQ=255)
                                                                                                                                         tgGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTgg  >  2:699259/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                            
CCTTTATAGCGGACTTCAGTCATAAAATGTGCATCCGGAGTACGTGTTAAAGGTACATTAGAGCCCCACATAATAATGTATGATGCGTTATACCAGTCACTTGATTCAGGCACATCTGTTTGCTCTCCCCAAATTTGTGGAGAGGCAGGTGGTAAATCTGCATACCAGTCATAAAAACTAAGCATTTCACCACCAAGCAAATTGATGAATCGAGCACCTGCTGCATAACTAATCATTGACATCGCTGGAATAGGTGTAAATCCTGCGATTCGATCTGG  >  CP000730/2519814‑2520091

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: