Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,945,181 G→T 94.8% pseudogene (280/663 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,945,1810GT94.8% 537.2 / 21.4 191pseudogene (280/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base T (69/112);  total (72/119)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 191 total aligned reads displayed.

CAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTC  >  CP009273/3945035‑3945325
                                                                                                                                                  |                                                                                                                                                
    tgtgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTc                                                                                                                                             <  1:892642/148‑1 (MQ=255)
    tgtgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTc                                                                                                                                             <  2:68096/148‑1 (MQ=255)
      tgGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTCCCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTg                                                                                                                                           >  2:1956463/1‑148 (MQ=255)
       ggggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGt                                                                                                                                          <  2:1356767/148‑1 (MQ=255)
        gggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTa                                                                                                                                         <  1:359907/148‑1 (MQ=255)
        gggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTa                                                                                                                                         <  1:2188789/148‑1 (MQ=255)
         ggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTNGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTAt                                                                                                                                        >  1:1805597/1‑148 (MQ=255)
         ggCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTAt                                                                                                                                        <  2:2208363/148‑1 (MQ=255)
           cagcaCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGTAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTAtct                                                                                                                                      <  2:1793671/148‑1 (MQ=255)
            agcaCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTAtctc                                                                                                                                     <  2:630603/148‑1 (MQ=255)
             gcaCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTcc                                                                                                                                    <  2:2739977/148‑1 (MQ=255)
             gcaCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTcc                                                                                                                                    >  1:1670767/1‑148 (MQ=255)
              caCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCg                                                                                                                                   >  1:1110194/1‑148 (MQ=255)
              caCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCg                                                                                                                                   >  1:699494/1‑148 (MQ=255)
              caCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCg                                                                                                                                   <  1:1016969/148‑1 (MQ=255)
              caCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCg                                                                                                                                   >  1:2272496/1‑148 (MQ=255)
               aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCgg                                                                                                                                  >  1:3080103/1‑148 (MQ=255)
               aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCgg                                                                                                                                  <  1:2355089/148‑1 (MQ=255)
               aCCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCgg                                                                                                                                  <  1:564593/148‑1 (MQ=255)
                ccAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGt                                                                                                                                 <  2:2156560/148‑1 (MQ=255)
                ccAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGt                                                                                                                                 >  1:2339139/1‑148 (MQ=255)
                  aGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGa                                                                                                                               >  2:1359367/1‑148 (MQ=255)
                       tgGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCTAACGATACCTTTGTCTGTATCTCCGGTGACGGct                                                                                                                          <  2:2963914/148‑1 (MQ=255)
                        gggCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGctc                                                                                                                         >  2:705156/1‑148 (MQ=255)
                         ggCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGctct                                                                                                                        >  2:1005187/1‑148 (MQ=255)
                         ggCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGctct                                                                                                                        >  2:1005115/1‑148 (MQ=255)
                         ggCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGctct                                                                                                                        <  1:6393/148‑1 (MQ=255)
                         ggCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGctct                                                                                                                        <  1:1003439/148‑1 (MQ=255)
                         ggCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTGTGTCTGTATCTCCGGTGACGGctct                                                                                                                        >  2:2550581/1‑148 (MQ=255)
                          gCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCtt                                                                                                                       <  1:257299/148‑1 (MQ=255)
                                cagcaCATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCatga                                                                                                                 <  1:3438240/148‑1 (MQ=255)
                                  gcaCATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCatgatg                                                                                                               <  1:3339720/148‑1 (MQ=255)
                                   cacaTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCAtgatga                                                                                                              >  2:1753786/1‑148 (MQ=255)
                                   cacaTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCAtgatga                                                                                                              <  2:536833/148‑1 (MQ=255)
                                   cacaTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCAtgatga                                                                                                              >  2:796082/1‑148 (MQ=255)
                                    acaTCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGaa                                                                                                             <  1:412595/148‑1 (MQ=255)
                                              cacTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAagag                                                                                                   <  1:2718653/148‑1 (MQ=255)
                                                cTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCt                                                                                                 <  2:996551/148‑1 (MQ=255)
                                                cTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCt                                                                                                 >  1:2584702/1‑148 (MQ=255)
                                                cTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCt                                                                                                 >  1:1321359/1‑148 (MQ=255)
                                                  cGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTgg                                                                                               >  1:1603665/1‑148 (MQ=255)
                                                     ccGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCa                                                                                            <  1:2242355/148‑1 (MQ=255)
                                                      cGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCAc                                                                                           >  1:2561900/1‑148 (MQ=255)
                                                         aaaaTTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGt                                                                                        >  1:225540/1‑148 (MQ=255)
                                                         aaaaTTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGt                                                                                        >  1:740098/1‑148 (MQ=255)
                                                            aTTTCATCACCTCCAGCGGTTTAGGGACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTaaa                                                                                     <  2:2969417/148‑1 (MQ=255)
                                                             tttCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTaaaa                                                                                    <  1:1859783/148‑1 (MQ=255)
                                                               tcatcaCCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACg                                                                                  >  1:2516307/1‑148 (MQ=255)
                                                                 atcaCCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCa                                                                                >  2:3475188/1‑148 (MQ=255)
                                                                   caCCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAg                                                                              <  2:1050471/148‑1 (MQ=255)
                                                                      cTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAg                                                                           <  2:2034952/148‑1 (MQ=255)
                                                                          aGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTAc                                                                       <  1:422762/148‑1 (MQ=255)
                                                                              gTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGtt                                                                   >  2:460188/1‑148 (MQ=255)
                                                                                ttAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGa                                                                 >  2:1638068/1‑148 (MQ=255)
                                                                                  aGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGaaa                                                               <  1:1904160/148‑1 (MQ=255)
                                                                                  aGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGaaa                                                               <  1:1276719/148‑1 (MQ=255)
                                                                                   ggTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGaaaa                                                              <  2:1763193/148‑1 (MQ=255)
                                                                                   ggTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGaaaa                                                              >  1:3251419/1‑148 (MQ=255)
                                                                                      aCCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCg                                                           <  2:280789/148‑1 (MQ=255)
                                                                                       ccATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGt                                                          >  2:1548949/1‑148 (MQ=255)
                                                                                        cATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTc                                                         <  1:1005115/148‑1 (MQ=255)
                                                                                          tGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCtt                                                       <  2:3448217/148‑1 (MQ=255)
                                                                                           gggTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTa                                                      >  1:1240669/1‑148 (MQ=255)
                                                                                           gggTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTa                                                      >  1:2729366/1‑148 (MQ=255)
                                                                                               tttGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCg                                                  >  2:1506742/1‑148 (MQ=255)
                                                                                                  ggTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATa                                               <  2:1414768/148‑1 (MQ=255)
                                                                                                   gTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATaa                                              >  1:1599399/1‑148 (MQ=255)
                                                                                                     ttACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAAcc                                            >  1:3050836/1‑148 (MQ=255)
                                                                                                       aCCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCaa                                          <  2:2561916/148‑1 (MQ=255)
                                                                                                                ggTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAggg                                 <  1:1359367/148‑1 (MQ=255)
                                                                                                                ggTTGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAggg                                 <  1:2124083/148‑1 (MQ=255)
                                                                                                                  ttGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGAt                               <  2:22639/148‑1 (MQ=255)
                                                                                                                   tGGCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATg                              <  2:579513/148‑1 (MQ=255)
                                                                                                                    ggCGCACAAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATgg                             <  2:2144563/148‑1 (MQ=255)
                                                                                                                        cacaAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCg                         <  2:2635100/148‑1 (MQ=255)
                                                                                                                          caAGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGAc                       >  1:2218111/1‑148 (MQ=255)
                                                                                                                           aaGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACa                      <  2:3097666/148‑1 (MQ=255)
                                                                                                                            aGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACaa                     <  1:529399/148‑1 (MQ=255)
                                                                                                                            aGTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACaa                     <  1:3322900/148‑1 (MQ=255)
                                                                                                                             gTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAAt                    <  1:1965870/148‑1 (MQ=255)
                                                                                                                             gTCGCGCGACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAAt                    >  2:383981/1‑148 (MQ=255)
                                                                                                                                gcgcgACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGgc                 <  2:3262608/148‑1 (MQ=255)
                                                                                                                                 cgcgACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGgca                >  1:2583219/1‑148 (MQ=255)
                                                                                                                                 cgcgACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGgca                >  2:2542233/1‑148 (MQ=255)
                                                                                                                                  gcgACCGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGgcag               <  2:1904166/148‑1 (MQ=255)
                                                                                                                                  gcgACCGAACGAGACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGgcag               <  1:1515674/148‑1 (MQ=255)
                                                                                                                                      ccGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAAc           >  1:291054/1‑148 (MQ=255)
                                                                                                                                      ccGAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAAc           >  1:3157972/1‑148 (MQ=255)
                                                                                                                                       cGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACt          >  1:2530143/1‑148 (MQ=255)
                                                                                                                                        gAACGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTg         >  1:3101773/1‑148 (MQ=255)
                                                                                                                                         aaCGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGt        <  2:2341962/148‑1 (MQ=255)
                                                                                                                                         aaCGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGt        >  1:2128084/1‑148 (MQ=255)
                                                                                                                                         aaCGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGt        <  2:1321359/148‑1 (MQ=255)
                                                                                                                                           cGATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGttt      >  1:926455/1‑148 (MQ=255)
                                                                                                                                            gATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttt     <  1:1487652/148‑1 (MQ=255)
                                                                                                                                            gATACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttt     <  1:1487489/148‑1 (MQ=255)
                                                                                                                                              tACCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGtttttt   <  2:713434/148‑1 (MQ=255)
                                                                                                                                               aCCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTc  <  2:405095/148‑1 (MQ=255)
                                                                                                                                               aCCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTc  >  2:1047044/1‑148 (MQ=255)
                                                                                                                                               aCCTTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTc  >  1:3280858/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CAGATGTGGGGCAGCACCAGATGTGGGCTGCGCAGCACATCGCCCACACTCGCCCGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCAAGAGCTGGGCACCGTAAAACGCAAGCAGTTACCGTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTC  >  CP009273/3945035‑3945325

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