Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,889,129 C→A 100% A4S (GCA→TCA)  ygcQ ← putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,889,1290CA96.8% 966.5 / 18.2 310A4S (GCA→TCA) ygcQputative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/8);  new base A (135/165);  total (137/173)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 312 total aligned reads displayed.

GCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCGGC  >  CP009273/2888982‑2889275
                                                                                                                                                   |                                                                                                                                                  
  cACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGata                                                                                                                                                <  2:3327534/150‑1 (MQ=255)
    cTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATAtt                                                                                                                                              <  2:1305954/150‑1 (MQ=255)
       ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTCCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTCTTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCAt                                                                                                                                           >  2:1467175/1‑150 (MQ=255)
       ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCAt                                                                                                                                           <  2:5209129/150‑1 (MQ=255)
        tCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATg                                                                                                                                          <  1:387075/150‑1 (MQ=255)
         ccaccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                         <  1:5077947/150‑1 (MQ=255)
           accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGCTATTCATGGCt                                                                                                                                       <  1:2732371/150‑1 (MQ=255)
            ccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTg                                                                                                                                      <  1:3731217/150‑1 (MQ=255)
            ccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTg                                                                                                                                      <  2:284967/150‑1 (MQ=255)
               aCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCAt                                                                                                                                   <  2:4049196/150‑1 (MQ=255)
               aCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCAt                                                                                                                                   >  1:1167143/1‑150 (MQ=255)
                  cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                                >  1:911307/1‑150 (MQ=255)
                  cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                                >  2:4606896/1‑150 (MQ=255)
                    tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCttt                                                                                                                              <  2:209004/150‑1 (MQ=255)
                     ccAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTg                                                                                                                             <  1:5215372/150‑1 (MQ=255)
                       aGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGc                                                                                                                           <  2:3932720/150‑1 (MQ=255)
                       aGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGc                                                                                                                           >  2:3805358/1‑150 (MQ=255)
                        gCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCg                                                                                                                          <  2:1634204/150‑1 (MQ=255)
                        gCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCg                                                                                                                          <  2:2878422/150‑1 (MQ=255)
                         cACCTGTTCCGCCTCTACAGGCTGAGGTTCAATCTGCCTGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGTCATTCTCCTGATTTATTGTCACACTTGAGATATTCATGGCTTCATCCTTTGGCGt                                                                                                                         <  2:2350045/150‑1 (MQ=255)
                             tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAAt                                                                                                                     <  1:2885935/150‑1 (MQ=255)
                             tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAAt                                                                                                                     >  1:3022072/1‑150 (MQ=255)
                             tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAAt                                                                                                                     <  1:5163909/150‑1 (MQ=255)
                              gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATa                                                                                                                    >  1:3463288/1‑150 (MQ=255)
                              gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATa                                                                                                                    >  2:3437559/1‑150 (MQ=255)
                                tCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGt                                                                                                                  <  2:4208860/150‑1 (MQ=255)
                                   gCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCt                                                                                                               <  1:2621479/150‑1 (MQ=255)
                                     cACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGc                                                                                                             <  2:1575445/150‑1 (MQ=255)
                                     cACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGc                                                                                                             <  2:2154348/150‑1 (MQ=255)
                                      aCTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGcc                                                                                                            <  2:913790/150‑1 (MQ=255)
                                       cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCa                                                                                                           >  1:2346024/1‑150 (MQ=255)
                                       cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCa                                                                                                           <  1:4462517/150‑1 (MQ=255)
                                        tACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCAt                                                                                                          <  2:3129364/150‑1 (MQ=255)
                                        tACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCAt                                                                                                          <  2:598877/150‑1 (MQ=255)
                                         aCAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATa                                                                                                         <  1:1467185/150‑1 (MQ=255)
                                          cAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                        <  1:4628292/150‑1 (MQ=255)
                                            ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCt                                                                                                      <  1:3921443/150‑1 (MQ=255)
                                            ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCt                                                                                                      <  2:1708407/150‑1 (MQ=255)
                                             gCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                                     <  2:4355153/150‑1 (MQ=255)
                                              cTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCttt                                                                                                    >  2:3310458/1‑150 (MQ=255)
                                              cTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCttt                                                                                                    <  1:2144064/150‑1 (MQ=255)
                                              cTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCttt                                                                                                    <  2:184761/150‑1 (MQ=255)
                                                gAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTg                                                                                                  >  1:3126381/1‑150 (MQ=255)
                                                      cattCTTCCAGTGCGCCAGGTTACAACCTCTTAAATCCTGCGCCGCCAGCCAGCTTGATATTTATTCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTTCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                            <  1:4176824/147‑1 (MQ=255)
                                                      cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                            >  2:3618826/1‑150 (MQ=255)
                                                      cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                            >  2:307045/1‑150 (MQ=255)
                                                      cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                            >  1:1349431/1‑150 (MQ=255)
                                                       aaTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCtttt                                                                                           <  1:4013013/150‑1 (MQ=255)
                                                        aTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCTGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttttt                                                                                          <  1:1813686/150‑1 (MQ=255)
                                                          cTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCt                                                                                        <  2:1667506/150‑1 (MQ=255)
                                                           tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTg                                                                                       <  2:4236893/150‑1 (MQ=255)
                                                               aGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTAc                                                                                   >  1:759359/1‑150 (MQ=255)
                                                               aGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTAc                                                                                   >  1:1099563/1‑150 (MQ=255)
                                                                  gcgcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGt                                                                                <  2:3840974/150‑1 (MQ=255)
                                                                  gcgcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGt                                                                                <  2:470700/150‑1 (MQ=255)
                                                                    gcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGttt                                                                              <  2:3517839/150‑1 (MQ=255)
                                                                    gcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGttt                                                                              >  2:4449718/1‑150 (MQ=255)
                                                                             cAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaa                                                                     >  1:4685503/1‑150 (MQ=255)
                                                                             cAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaa                                                                     >  2:3029478/1‑150 (MQ=255)
                                                                                ccGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatca                                                                  <  1:117107/150‑1 (MQ=255)
                                                                                   cTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCg                                                               >  2:5319266/1‑150 (MQ=255)
                                                                                         tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTc                                                         <  2:3582201/150‑1 (MQ=255)
                                                                                         tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTc                                                         >  1:2290771/1‑150 (MQ=255)
                                                                                                  gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGtt                                                >  1:2989468/1‑150 (MQ=255)
                                                                                                      gccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctc                                            >  1:670447/1‑150 (MQ=255)
                                                                                                       ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctct                                           <  2:4730166/150‑1 (MQ=255)
                                                                                                             tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCtt                                     <  2:1241101/150‑1 (MQ=255)
                                                                                                             tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCtt                                     >  2:1427831/1‑150 (MQ=255)
                                                                                                             tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCtt                                     >  2:2482861/1‑150 (MQ=255)
                                                                                                             tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCtt                                     <  2:2276394/150‑1 (MQ=255)
                                                                                                              gCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTg                                    <  2:3979554/150‑1 (MQ=255)
                                                                                                              gCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTg                                    <  2:4531825/150‑1 (MQ=255)
                                                                                                                        cATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAg                          <  1:3287818/150‑1 (MQ=255)
                                                                                                                         aTTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCTGa                         >  2:4094633/1‑150 (MQ=255)
                                                                                                                         aTTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                         >  2:497651/1‑150 (MQ=255)
                                                                                                                         aTTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                         >  2:2990056/1‑150 (MQ=255)
                                                                                                                         aTTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                         >  2:5006463/1‑150 (MQ=255)
                                                                                                                           tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGaca                       >  2:3241143/1‑150 (MQ=255)
                                                                                                                           tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGaca                       >  1:4856128/1‑150 (MQ=255)
                                                                                                                            ctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGacac                      <  1:3805358/150‑1 (MQ=255)
                                                                                                                             tcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTACATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                     <  2:2547533/150‑1 (MQ=255)
                                                                                                                             tcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                     <  1:1226148/150‑1 (MQ=255)
                                                                                                                             tcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                     <  2:2761380/150‑1 (MQ=255)
                                                                                                                             tcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                     >  1:251641/1‑150 (MQ=255)
                                                                                                                             tcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACt                     <  1:3618826/150‑1 (MQ=255)
                                                                                                                              ccTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTTAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTg                    >  1:1074785/1‑150 (MQ=255)
                                                                                                                              ccTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTg                    >  1:1941720/1‑150 (MQ=255)
                                                                                                                              ccTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTg                    >  1:3425900/1‑150 (MQ=255)
                                                                                                                                tGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa                  >  2:1210154/1‑150 (MQ=255)
                                                                                                                                 gATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCAt                 >  1:1631361/1‑150 (MQ=255)
                                                                                                                                     tATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCt             >  1:3776257/1‑150 (MQ=255)
                                                                                                                                      aTTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTg            <  2:5129814/150‑1 (MQ=255)
                                                                                                                                       ttGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTACTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTgg           >  1:3562342/1‑150 (MQ=255)
                                                                                                                                        tGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGc          <  2:4665830/150‑1 (MQ=255)
                                                                                                                                          tCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAt        <  2:729693/150‑1 (MQ=255)
                                                                                                                                           cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAtt       >  1:1587566/1‑150 (MQ=255)
                                                                                                                                              aaTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCg    <  1:256317/150‑1 (MQ=255)
                                                                                                                                              aaTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCg    >  1:4368130/1‑150 (MQ=255)
                                                                                                                                               aTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCgg   >  2:461992/1‑150 (MQ=255)
                                                                                                                                                ttGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTACCAGTTGCAGACACTGCATCGCTGGCATTTCGgc  >  2:2607310/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCGGC  >  CP009273/2888982‑2889275

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: