Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,989,295 C→A 91.0% pseudogene (85/261 nt) ygeN ← pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,989,2950CA91.0% 547.2 / 41.6 178pseudogene (85/261 nt)ygeNpseudogene
Reads supporting (aligned to +/- strand):  ref base C (9/7);  new base A (64/98);  total (73/105)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 182 total aligned reads displayed.

AGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTCAACA  >  CP009273/2989148‑2989441
                                                                                                                                                   |                                                                                                                                                   
   gTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCa                                                                                                                                                <  2:2603030/150‑1 (MQ=255)
    tataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCaa                                                                                                                                               <  1:883461/150‑1 (MQ=255)
     ataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAAt                                                                                                                                              >  2:5059773/1‑150 (MQ=255)
     ataAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAAt                                                                                                                                              <  1:1681657/150‑1 (MQ=255)
          aTTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATc                                                                                                                                         <  1:4548756/150‑1 (MQ=255)
          aTTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATc                                                                                                                                         <  2:4031968/150‑1 (MQ=255)
            tATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAg                                                                                                                                       <  2:3115459/150‑1 (MQ=255)
                 aTCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAAtt                                                                                                                                  >  2:2181741/1‑150 (MQ=255)
                     catcaATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTc                                                                                                                              >  1:3470747/1‑150 (MQ=255)
                      atcaATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCt                                                                                                                             >  1:3815494/1‑150 (MQ=255)
                       tcaATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCtt                                                                                                                            <  1:904822/150‑1 (MQ=255)
                        caATACTTCCTTTTGTATATCCATCATAATAACCTCGAATAAATTTTTCTTCTCTCTTTTCCTCTATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCttt                                                                                                                           <  2:3495166/150‑1 (MQ=255)
                        caATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCttt                                                                                                                           >  1:477093/1‑150 (MQ=255)
                        caATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCttt                                                                                                                           <  2:192769/150‑1 (MQ=255)
                            aatGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAAt                                                                                                                       <  2:520566/150‑1 (MQ=255)
                             atGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTTTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATc                                                                                                                      <  1:4069880/150‑1 (MQ=255)
                             atGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATc                                                                                                                      >  1:3496724/1‑150 (MQ=255)
                              tGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATca                                                                                                                     <  1:2184674/150‑1 (MQ=255)
                                 cTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCacaa                                                                                                                  >  1:3899052/1‑150 (MQ=255)
                                     tGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACAcc                                                                                                              <  1:4683523/150‑1 (MQ=255)
                                        atatCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTTCGGTTGTTTCAAGAATGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAtt                                                                                                           <  2:362683/150‑1 (MQ=255)
                                        ataGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAtt                                                                                                           >  1:2472424/1‑150 (MQ=255)
                                          aGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTATTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTATTCACAACACCATTTg                                                                                                         >  2:3575219/1‑150 (MQ=255)
                                              atcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAgg                                                                                                     >  2:4226301/1‑150 (MQ=255)
                                              atcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAgg                                                                                                     >  1:196011/1‑150 (MQ=255)
                                                  aaataaCCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAt                                                                                                 >  1:3322952/2‑150 (MQ=255)
                                                   aataaCCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAtt                                                                                                >  1:2307357/1‑150 (MQ=255)
                                                     taaCCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAtttt                                                                                              <  1:2875526/150‑1 (MQ=255)
                                                     taaCCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGAtttt                                                                                              <  2:1025140/150‑1 (MQ=255)
                                                      aaCCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTa                                                                                             <  2:1581699/150‑1 (MQ=255)
                                                        ccTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATc                                                                                           >  2:4233313/1‑150 (MQ=255)
                                                           tGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAg                                                                                        >  1:1450256/1‑150 (MQ=255)
                                                             aaTAAATTTTTCTTCTCTCTTGTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCt                                                                                      <  1:3441662/150‑1 (MQ=255)
                                                                aaaTTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCAt                                                                                   <  2:442302/150‑1 (MQ=255)
                                                                  aTTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCAtct                                                                                 <  2:3496724/150‑1 (MQ=255)
                                                                  aTTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCAtct                                                                                 >  2:3128274/1‑150 (MQ=255)
                                                                   tttttCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCAtctc                                                                                <  2:2153747/150‑1 (MQ=255)
                                                                    ttttCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATctct                                                                               >  2:2381052/1‑150 (MQ=255)
                                                                    ttttCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATctct                                                                               <  2:1043992/150‑1 (MQ=255)
                                                                     tttCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCATCTCTa                                                                              >  1:2463944/1‑150 (MQ=255)
                                                                       tcttctCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAtt                                                                           >  2:4704349/1‑150 (MQ=255)
                                                                        cttctCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAttt                                                                          <  1:4199981/150‑1 (MQ=255)
                                                                        cttctCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAttt                                                                          >  2:3347643/1‑150 (MQ=255)
                                                                          tctctcTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGATCATCTCTAttttt                                                                        >  1:1328215/1‑150 (MQ=255)
                                                                               ctTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCg                                                                   >  1:857184/1‑150 (MQ=255)
                                                                               ctTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCg                                                                   <  1:4740864/150‑1 (MQ=255)
                                                                               ctTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCg                                                                   >  2:3695280/1‑150 (MQ=255)
                                                                                ttttCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGc                                                                  <  1:5109398/150‑1 (MQ=255)
                                                                                ttttCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGc                                                                  >  2:898405/1‑150 (MQ=255)
                                                                                  ttCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCAt                                                                <  1:2377607/150‑1 (MQ=255)
                                                                                  ttCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCAt                                                                >  1:3779133/1‑150 (MQ=255)
                                                                                      tCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATaaaa                                                            >  1:1734506/1‑150 (MQ=255)
                                                                                       cAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAATCTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAAt                                                           >  1:4334628/1‑150 (MQ=255)
                                                                                       cAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAAt                                                           <  2:1147525/150‑1 (MQ=255)
                                                                                       cAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAAt                                                           <  2:1733028/150‑1 (MQ=255)
                                                                                        aaTACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATg                                                          <  1:2088790/150‑1 (MQ=255)
                                                                                           aCTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTaa                                                       <  1:4382236/150‑1 (MQ=255)
                                                                                            ctctGAATTTTTGCGGTTGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAAt                                                      <  2:2671215/150‑1 (MQ=255)
                                                                                            ctctGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAAt                                                      <  1:5059773/150‑1 (MQ=255)
                                                                                            ctctGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAAt                                                      <  1:4245628/150‑1 (MQ=255)
                                                                                                gAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATtctc                                                  <  2:1389373/150‑1 (MQ=255)
                                                                                                   tttttGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTctctct                                               >  2:405114/1‑150 (MQ=255)
                                                                                                   tttttGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTctctct                                               >  2:405457/1‑150 (MQ=255)
                                                                                                    ttttGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTa                                              <  1:1694204/150‑1 (MQ=255)
                                                                                                    ttttGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTa                                              <  2:2539342/150‑1 (MQ=255)
                                                                                                    ttttGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTa                                              <  1:3686120/150‑1 (MQ=255)
                                                                                                    ttttGCGGTGGTTTCAAGAAGGATTTTTTATTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTa                                              <  2:4463440/150‑1 (MQ=255)
                                                                                                      ttGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTaaa                                            <  1:1359755/150‑1 (MQ=255)
                                                                                                       tGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTaaaa                                           >  2:2766648/1‑150 (MQ=255)
                                                                                                         cGGTGGTTTCAAGAAGGATTTTTTCTTTAAGATAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAgt                                         <  2:2012005/150‑1 (MQ=255)
                                                                                                          ggtggtTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAgtg                                        <  1:1448270/150‑1 (MQ=255)
                                                                                                          ggtggtTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTAGCATAAAAAGTAATTCTCTCTAAAAgtg                                        >  1:4438871/1‑150 (MQ=255)
                                                                                                            tggtTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGtt                                      <  1:3736429/150‑1 (MQ=255)
                                                                                                              gtTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGa                                    <  1:4723997/150‑1 (MQ=255)
                                                                                                               tttCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGAt                                   <  1:3592777/150‑1 (MQ=255)
                                                                                                                 tCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGc                                 <  1:1530679/150‑1 (MQ=255)
                                                                                                                    agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGAt                                   <  2:1791012/145‑1 (MQ=255)
                                                                                                                    agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGAt                                   >  1:1791012/1‑145 (MQ=255)
                                                                                                                    agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAg                              >  1:565282/1‑150 (MQ=255)
                                                                                                                    agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAg                              >  1:3172994/1‑150 (MQ=255)
                                                                                                                        ggATTTTTTTTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGt                          >  1:3962819/1‑150 (MQ=255)
                                                                                                                         gATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTa                         >  2:3923132/1‑150 (MQ=255)
                                                                                                                          aTTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTaa                        >  1:5065568/1‑150 (MQ=255)
                                                                                                                            tttttCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACg                      <  1:2043699/150‑1 (MQ=255)
                                                                                                                                  tttAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGAtt                <  2:4574005/150‑1 (MQ=255)
                                                                                                                                   ttAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTg               >  1:4106382/1‑150 (MQ=255)
                                                                                                                                     aaGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTa             >  1:1679269/1‑150 (MQ=255)
                                                                                                                                      agatAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTa             <  2:4924763/150‑1 (MQ=255)
                                                                                                                                      agagAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTa             <  2:4924401/150‑1 (MQ=255)
                                                                                                                                       gagaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtg           >  1:2560509/1‑150 (MQ=255)
                                                                                                                                       gagaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtg           <  2:3238097/150‑1 (MQ=255)
                                                                                                                                       gagaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtg           >  1:4095331/1‑150 (MQ=255)
                                                                                                                                        agaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtgt          >  1:980778/1‑150 (MQ=255)
                                                                                                                                        agaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtgt          <  2:633561/150‑1 (MQ=255)
                                                                                                                                         gaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTTCAGGAGATTTTATCAAGCTCATATCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtgtg         <  2:4306993/150‑1 (MQ=255)
                                                                                                                                         gaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtgtg         <  2:4307134/150‑1 (MQ=255)
                                                                                                                                          aaaCTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTgg        <  2:1449953/150‑1 (MQ=255)
                                                                                                                                            aCTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTc      >  2:3781674/1‑150 (MQ=255)
                                                                                                                                            aCTATTTATTTCAATCCATCAGAAATGTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTc      <  1:4661983/150‑1 (MQ=255)
                                                                                                                                                tttatttcAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTcaaca  <  2:1258262/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTCAACA  >  CP009273/2989148‑2989441

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: