Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,113,450 A→T 100% I161I (ATT→ATA yceI ← secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,113,4500AT98.0% 815.3 / 4.2 248I161I (ATT→ATAyceIsecreted protein
Reads supporting (aligned to +/- strand):  ref base A (1/4);  new base T (131/112);  total (132/116)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 248 total aligned reads displayed.

CAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTC  >  CP009273/1113305‑1113594
                                                                                                                                                 |                                                                                                                                                
  aCAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAttt                                                                                                                                              <  1:1474625/148‑1 (MQ=255)
   cAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTAtttt                                                                                                                                             >  1:727925/1‑148 (MQ=255)
     gAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGc                                                                                                                                           <  1:355943/148‑1 (MQ=255)
      aataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGcc                                                                                                                                          <  1:949931/148‑1 (MQ=255)
      aataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGcc                                                                                                                                          >  2:2692379/1‑148 (MQ=255)
       ataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCt                                                                                                                                         <  1:1132179/148‑1 (MQ=255)
       ataaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCt                                                                                                                                         <  1:2027491/148‑1 (MQ=255)
         aaataTCATTGCCGGATTCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTc                                                                                                                                       >  1:1415433/1‑148 (MQ=255)
         aaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTc                                                                                                                                       >  2:645387/1‑148 (MQ=255)
         aaataTCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTGCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTc                                                                                                                                       >  2:666879/1‑148 (MQ=255)
           atatCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCgg                                                                                                                                     >  2:3166861/1‑148 (MQ=255)
              tCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCt                                                                                                                                  >  1:272587/1‑148 (MQ=255)
              tCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCt                                                                                                                                  >  1:567504/1‑148 (MQ=255)
                 ttGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGa                                                                                                                               >  2:364639/1‑148 (MQ=255)
                  tGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGaa                                                                                                                              >  1:154857/1‑148 (MQ=255)
                    ccGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGc                                                                                                                            <  2:1702354/148‑1 (MQ=255)
                        aTGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGc                                                                                                                        >  1:1788222/1‑148 (MQ=255)
                        aTGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGc                                                                                                                        <  1:2053740/148‑1 (MQ=255)
                        aTGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGATTTATTTTGACTTCTGCCTCGAAGCCTGc                                                                                                                        >  2:3434973/1‑148 (MQ=255)
                         tGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCa                                                                                                                       >  1:797440/1‑148 (MQ=255)
                            gTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGt                                                                                                                    <  2:2570454/148‑1 (MQ=255)
                             tGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGTACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGtt                                                                                                                   <  1:3409648/148‑1 (MQ=255)
                                 aTCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTacc                                                                                                               <  2:743032/148‑1 (MQ=255)
                                  tCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTacca                                                                                                              <  1:1826576/148‑1 (MQ=255)
                                  tCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTacca                                                                                                              >  2:49363/1‑148 (MQ=255)
                                      gCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCAcccc                                                                                                          >  2:2821914/1‑148 (MQ=255)
                                        aaCATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCAt                                                                                                        <  2:1372676/148‑1 (MQ=255)
                                         aCATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATg                                                                                                       <  2:2927489/148‑1 (MQ=255)
                                         aCATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATg                                                                                                       >  2:1069915/1‑148 (MQ=255)
                                            tttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgt                                                                                                    >  2:2208132/1‑148 (MQ=255)
                                             ttCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtc                                                                                                   <  1:3368139/148‑1 (MQ=255)
                                              tCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtcg                                                                                                  <  1:1612333/148‑1 (MQ=255)
                                               cACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGgtcgt                                                                                                 <  1:1875254/148‑1 (MQ=255)
                                                    ttACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGccc                                                                                            <  1:421116/148‑1 (MQ=255)
                                                             ctgTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaatt                                                                                   >  1:2975325/1‑148 (MQ=255)
                                                             ctgTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaatt                                                                                   <  2:1741120/148‑1 (MQ=255)
                                                                tACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaat                                                                                <  2:845553/148‑1 (MQ=255)
                                                                tACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaat                                                                                <  2:408304/148‑1 (MQ=255)
                                                                tACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaat                                                                                <  1:1764741/148‑1 (MQ=255)
                                                                tACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaat                                                                                <  2:286818/148‑1 (MQ=255)
                                                                tACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCaattaat                                                                                >  1:2567720/1‑148 (MQ=255)
                                                                   gCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTc                                                                             <  2:2774873/148‑1 (MQ=255)
                                                                      ttCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCt                                                                          >  2:2034334/1‑148 (MQ=255)
                                                                       tCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCtt                                                                         >  1:2843761/1‑148 (MQ=255)
                                                                       tCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCtt                                                                         <  1:202254/148‑1 (MQ=255)
                                                                        ccACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTc                                                                        <  1:3112865/148‑1 (MQ=255)
                                                                        ccACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTc                                                                        >  2:3126244/1‑148 (MQ=255)
                                                                            tGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGc                                                                    >  2:2364116/1‑148 (MQ=255)
                                                                             gAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCg                                                                   <  1:479486/148‑1 (MQ=255)
                                                                               aataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTg                                                                 <  1:807113/148‑1 (MQ=255)
                                                                               aataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTg                                                                 <  1:2589259/148‑1 (MQ=255)
                                                                                ataatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGa                                                                <  1:1006022/148‑1 (MQ=255)
                                                                                 taatCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGAc                                                               >  2:2248780/1‑148 (MQ=255)
                                                                                   atCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACgg                                                             >  1:619570/1‑148 (MQ=255)
                                                                                     cAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGt                                                           >  2:554391/1‑148 (MQ=255)
                                                                                     cAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGt                                                           <  2:2388609/148‑1 (MQ=255)
                                                                                      aGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGtt                                                          >  2:1196012/1‑148 (MQ=255)
                                                                                       gATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGttt                                                         >  2:3043163/1‑148 (MQ=255)
                                                                                        aTCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTg                                                        >  2:3416821/1‑148 (MQ=255)
                                                                                        aTCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTg                                                        >  1:1250469/1‑148 (MQ=255)
                                                                                         tCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTgg                                                       >  2:610511/1‑148 (MQ=255)
                                                                                           cACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGtt                                                     >  2:2150261/1‑148 (MQ=255)
                                                                                            aCTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTa                                                    >  2:3538678/1‑148 (MQ=255)
                                                                                                    agaAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTc                                            >  2:2229608/1‑148 (MQ=255)
                                                                                                    agaAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTc                                            >  2:755002/1‑148 (MQ=255)
                                                                                                        gCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCg                                        >  2:166120/1‑148 (MQ=255)
                                                                                                        gCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCg                                        >  2:1970110/1‑148 (MQ=255)
                                                                                                           ggACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCa                                     >  2:1409811/1‑148 (MQ=255)
                                                                                                           ggACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCa                                     <  2:1733593/148‑1 (MQ=255)
                                                                                                           ggACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCa                                     <  2:1605386/148‑1 (MQ=255)
                                                                                                            gACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAg                                    >  1:3017879/1‑148 (MQ=255)
                                                                                                               cTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATc                                 >  2:3214151/1‑148 (MQ=255)
                                                                                                               cTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATc                                 >  2:1444229/1‑148 (MQ=255)
                                                                                                                tAAATCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCg                                <  1:1230784/148‑1 (MQ=255)
                                                                                                                 aaaTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGc                               >  2:1286125/1‑148 (MQ=255)
                                                                                                                 aaaTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGc                               >  2:3496780/1‑148 (MQ=255)
                                                                                                                  aaTCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGcc                              >  1:655678/1‑148 (MQ=255)
                                                                                                                    tCTGTCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCgg                            <  2:2567825/148‑1 (MQ=255)
                                                                                                                        tCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAat                        <  1:229358/148‑1 (MQ=255)
                                                                                                                        tCTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAat                        >  2:3203782/1‑148 (MQ=255)
                                                                                                                         cTTGATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAata                       >  2:2521861/1‑148 (MQ=255)
                                                                                                                            gATATTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATcc                    <  2:1135125/148‑1 (MQ=255)
                                                                                                                             atatTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCa                   >  2:1720444/1‑148 (MQ=255)
                                                                                                                               atTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGt                 <  1:2875653/148‑1 (MQ=255)
                                                                                                                               atTGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTACAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGt                 >  2:184578/1‑148 (MQ=255)
                                                                                                                                ttGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGtt                <  2:797719/148‑1 (MQ=255)
                                                                                                                                ttGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGtt                <  2:2469413/148‑1 (MQ=255)
                                                                                                                                 tGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTc               >  2:619721/1‑148 (MQ=255)
                                                                                                                                 tGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTc               >  2:2404798/1‑148 (MQ=255)
                                                                                                                                 tGAAGTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTc               <  2:1299247/148‑1 (MQ=255)
                                                                                                                                     gTCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCa           <  1:3142490/148‑1 (MQ=255)
                                                                                                                                      tCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCAc          <  1:1634180/148‑1 (MQ=255)
                                                                                                                                      tCCTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCAc          <  1:1634952/148‑1 (MQ=255)
                                                                                                                                       ccTTGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCAcc         >  1:233133/1‑148 (MQ=255)
                                                                                                                                          tGAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtc      >  1:1312344/1‑148 (MQ=255)
                                                                                                                                           gAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtct     <  1:2370030/148‑1 (MQ=255)
                                                                                                                                           gAGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtct     <  1:503867/148‑1 (MQ=255)
                                                                                                                                            aGCTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtctt    >  1:2391955/1‑148 (MQ=255)
                                                                                                                                              cTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtctttc  <  2:3026566/148‑1 (MQ=255)
                                                                                                                                              cTTTATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGtctttc  >  1:1940168/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAGCTGGACCTAAATCTGTCTTGATATTGAAGTCCTTGAGCTTAATTTTGCCTTCGGCCTCGAAGCCTGCACGTTTACCACCCCATGGGTCGTCGCCCTGACCAATTAATTTCGCTTCCAGCGTGACGGGTTTGGTTACGCCATTCAGCGTCAGATCGCCGGTAATATCCAGTTCGTCACCGTCTTTC  >  CP009273/1113305‑1113594

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: