Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 250,909 A→T 100% L432Q (CTG→CAG)  pepD ← aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273250,9090AT99.3% 901.4 / ‑2.8 282L432Q (CTG→CAG) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (161/117);  total (162/120)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
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 282 total aligned reads displayed.

TCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTT  >  CP009273/250762‑251053
                                                                                                                                                   |                                                                                                                                                
tctttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTgg                                                                                                                                                >  2:1076476/1‑150 (MQ=255)
  tttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGcc                                                                                                                                              >  2:5227847/1‑150 (MQ=255)
     cagcagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCcgc                                                                                                                                           >  2:1438407/1‑150 (MQ=255)
        cagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGAGACGGTGTTTTGAACAGACCACATTCCTGGCCCGCGTg                                                                                                                                        >  2:3032311/1‑150 (MQ=255)
            tcagtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATa                                                                                                                                    <  2:4410896/150‑1 (MQ=255)
              agtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAAt                                                                                                                                  <  1:3789085/150‑1 (MQ=255)
                tcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCt                                                                                                                                >  2:4912897/1‑150 (MQ=255)
                tcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCt                                                                                                                                >  2:2893838/1‑150 (MQ=255)
                  agcagTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCttg                                                                                                                              >  2:1711674/1‑148 (MQ=255)
                    cagTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGAt                                                                                                                            <  2:3378452/150‑1 (MQ=255)
                    cagTGTCCAGGAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGAt                                                                                                                            <  2:327262/150‑1 (MQ=255)
                      gtgtCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGt                                                                                                                          <  1:2021289/150‑1 (MQ=255)
                       tgtCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGtt                                                                                                                         <  2:4977367/150‑1 (MQ=255)
                        gtCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTc                                                                                                                        <  2:1373119/150‑1 (MQ=255)
                         tCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCg                                                                                                                       >  2:2756855/1‑150 (MQ=255)
                         tCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCg                                                                                                                       >  2:4485871/1‑150 (MQ=255)
                         tCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCg                                                                                                                       <  1:2942604/150‑1 (MQ=255)
                         tCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCg                                                                                                                       <  2:116810/150‑1 (MQ=255)
                           cAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGc                                                                                                                     >  2:2737289/1‑150 (MQ=255)
                            aGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCg                                                                                                                    >  1:770342/1‑150 (MQ=255)
                            aGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCg                                                                                                                    >  2:3471785/1‑150 (MQ=255)
                            aGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCg                                                                                                                    <  1:3682068/150‑1 (MQ=255)
                               aaTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCt                                                                                                                 <  2:1768494/150‑1 (MQ=255)
                                aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGTCGTCtt                                                                                                                >  1:4922964/1‑150 (MQ=255)
                                aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCtt                                                                                                                >  2:383630/1‑150 (MQ=255)
                                    ccTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCttgttg                                                                                                            >  1:1647381/1‑150 (MQ=255)
                                     cTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGa                                                                                                           >  2:2188104/1‑150 (MQ=255)
                                      tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGaa                                                                                                          >  2:537900/1‑150 (MQ=255)
                                       aCGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAAc                                                                                                         >  1:1092940/1‑150 (MQ=255)
                                       aCGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAAc                                                                                                         <  1:3570515/150‑1 (MQ=255)
                                        cGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACa                                                                                                        >  2:464674/1‑150 (MQ=255)
                                         gCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAg                                                                                                       <  2:4694379/150‑1 (MQ=255)
                                          cTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAgg                                                                                                      <  2:2086809/150‑1 (MQ=255)
                                            ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGAACGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGTATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTTGATAATCTGGATGTTCGGCTTCTTGTTGAACAGgcg                                                                                                    <  2:182866/150‑1 (MQ=255)
                                                 tgtgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGAt                                                                                               <  2:225542/150‑1 (MQ=255)
                                                   tgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAg                                                                                             <  2:4129489/150‑1 (MQ=255)
                                                    gAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAGGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAgg                                                                                            <  2:1960776/150‑1 (MQ=255)
                                                       cTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGttt                                                                                         <  1:4416058/150‑1 (MQ=255)
                                                            tCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTCGAACAGGCGCTGATAGGTTTCACGt                                                                                    >  1:2666203/1‑150 (MQ=255)
                                                                 cGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAg                                                                               >  2:4273001/1‑150 (MQ=255)
                                                                  ggAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGa                                                                              >  1:4628838/1‑150 (MQ=255)
                                                                  ggAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGa                                                                              >  2:417487/1‑150 (MQ=255)
                                                                   gagaGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGAt                                                                             <  1:4275024/150‑1 (MQ=255)
                                                                    agagTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATg                                                                            >  1:2445168/1‑150 (MQ=255)
                                                                          tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAc                                                                      >  2:5150417/1‑150 (MQ=255)
                                                                           ggACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAcc                                                                     >  1:2775405/1‑150 (MQ=255)
                                                                           ggACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAcc                                                                     >  2:164148/1‑150 (MQ=255)
                                                                                ggTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGgaga                                                                >  1:586540/1‑150 (MQ=255)
                                                                                ggTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGgaga                                                                >  2:870523/1‑150 (MQ=255)
                                                                                  tGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAAt                                                              >  2:5161669/1‑150 (MQ=255)
                                                                                    aTAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTa                                                            <  2:5302063/150‑1 (MQ=255)
                                                                                       gTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCg                                                         <  1:848859/150‑1 (MQ=255)
                                                                                              cGATAGACACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCt                                                  <  2:3831061/150‑1 (MQ=255)
                                                                                              cGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCt                                                  <  1:1773208/150‑1 (MQ=255)
                                                                                               gATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTg                                                 >  1:1443954/1‑150 (MQ=255)
                                                                                                aTAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgc                                                >  1:1347942/1‑150 (MQ=255)
                                                                                                aTAGAAACCATGTCCATTTCCGGATACGGTAATTAGAACAGACCACATACCAGGCCCGCGTGGATAATCTGGATGTACTGCGACTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAACGGAGAATTAGCGTCCGGCTgc                                                >  2:2725681/1‑150 (MQ=255)
                                                                                                 tAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcc                                               >  1:2377554/1‑150 (MQ=255)
                                                                                                  aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcca                                              >  2:624749/1‑150 (MQ=255)
                                                                                                  aGAAACCAGGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcca                                              <  2:5044553/150‑1 (MQ=255)
                                                                                                   gAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccag                                             >  1:5231289/1‑150 (MQ=255)
                                                                                                   gAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccag                                             >  2:3921426/1‑150 (MQ=255)
                                                                                                    aaaCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagc                                            >  2:907958/1‑150 (MQ=255)
                                                                                                     aaCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcc                                           >  1:5210724/1‑150 (MQ=255)
                                                                                                      aCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcca                                          >  1:2544460/1‑150 (MQ=255)
                                                                                                      aCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcca                                          >  2:750914/1‑150 (MQ=255)
                                                                                                       ccATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGccagccag                                         >  2:3126277/1‑150 (MQ=255)
                                                                                                         aTGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGa                                       >  2:1558289/1‑150 (MQ=255)
                                                                                                          tGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGat                                      <  1:3394837/150‑1 (MQ=255)
                                                                                                          tGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGat                                      <  1:4362661/150‑1 (MQ=255)
                                                                                                          tGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGat                                      <  1:1273943/150‑1 (MQ=255)
                                                                                                            tCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGatat                                    <  1:4150851/150‑1 (MQ=255)
                                                                                                             ccATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTACACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATg                                   <  1:4989416/150‑1 (MQ=255)
                                                                                                              cATTTCCGGCTACGGTTTTTTGAACAGCCCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGATGATAGGTTTCACGTACCAGATGCATCACCGGAGATTTAGCGTCCGGCTGCCAGCCAGGATATgc                                  <  1:2148735/150‑1 (MQ=255)
                                                                                                               aTTTCCGGATACGTTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 <  1:5069959/150‑1 (MQ=255)
                                                                                                               aTTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 >  2:3077315/1‑150 (MQ=255)
                                                                                                               aTTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 >  1:46740/1‑150 (MQ=255)
                                                                                                               aTTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 >  2:1645615/1‑150 (MQ=255)
                                                                                                               aTTTCCGGATACGGTTTTTGGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcg                                 <  1:4808338/150‑1 (MQ=255)
                                                                                                                  tCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCt                              <  1:2893838/150‑1 (MQ=255)
                                                                                                                     ggATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTc                           <  1:715721/150‑1 (MQ=255)
                                                                                                                      gATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCg                          <  1:2235486/150‑1 (MQ=255)
                                                                                                                       aTACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGc                         >  2:3648269/1‑150 (MQ=255)
                                                                                                                        tACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCt                        <  1:1124535/150‑1 (MQ=255)
                                                                                                                               ttttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGttt                 <  2:1730093/150‑1 (MQ=255)
                                                                                                                                 ttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTc               >  1:2026083/1‑150 (MQ=255)
                                                                                                                                   gAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTcgc             >  1:4035128/1‑150 (MQ=255)
                                                                                                                                      cAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcc          <  1:494422/150‑1 (MQ=255)
                                                                                                                                      cAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcc          >  2:1097706/1‑150 (MQ=255)
                                                                                                                                        gACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccag        >  1:517060/1‑150 (MQ=255)
                                                                                                                                           cacaTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccagcca     >  2:2466096/1‑150 (MQ=255)
                                                                                                                                            acaTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGccagccag    >  2:4016285/1‑150 (MQ=255)
                                                                                                                                             caTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGt   <  1:1998184/150‑1 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  <  1:383630/150‑1 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  1:3048509/1‑150 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  2:2118448/1‑150 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  1:225108/1‑150 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  <  1:699037/150‑1 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  1:932223/1‑150 (MQ=255)
                                                                                                                                              aTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt  >  2:2901402/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTT  >  CP009273/250762‑251053

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: ‑