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,9090AT98.8% 818.4 / ‑0.7 250L432Q (CTG→CAG) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (118/128);  total (121/129)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
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 253 total aligned reads displayed.

TTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCC  >  CP009273/250760‑251058
                                                                                                                                                     |                                                                                                                                                     
tttctttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCa                                                                                                                                                       <  1:4903658/150‑1 (MQ=255)
  tctttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTgg                                                                                                                                                     <  1:358324/150‑1 (MQ=255)
   ctttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGc                                                                                                                                                    >  2:5099883/1‑150 (MQ=255)
     ttcAGCAGTTCAGTCAGCAGTGTCCAGTAATGTCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGAGAGTTGGCCCTATAGAAACCATGTCCATTTCCGGATACGGTGTTGTGAACAGTCCCCATTCCTGGccc                                                                                                                                                  >  2:4392378/1‑150 (MQ=255)
       cagcagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCcgc                                                                                                                                                >  1:3300002/1‑150 (MQ=255)
         gcagTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGt                                                                                                                                              >  2:2535240/1‑150 (MQ=255)
           agTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTgg                                                                                                                                            >  1:3204238/1‑150 (MQ=255)
           agTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTgg                                                                                                                                            <  1:2529845/150‑1 (MQ=255)
              tcagtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATa                                                                                                                                         >  1:1241863/1‑150 (MQ=255)
               cagtcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATaa                                                                                                                                        <  1:1778826/150‑1 (MQ=255)
                  tcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCt                                                                                                                                     >  2:3445529/1‑150 (MQ=255)
                  tcagCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCt                                                                                                                                     >  1:1420519/1‑150 (MQ=255)
                       agTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATg                                                                                                                                >  1:626623/1‑150 (MQ=255)
                            ccAGTAATGACCTACGCTTTCGATGTGAACGTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCgg                                                                                                                           <  1:4362910/150‑1 (MQ=255)
                               gTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGt                                                                                                                        >  2:3203105/1‑150 (MQ=255)
                               gTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGt                                                                                                                        >  2:3201850/1‑150 (MQ=255)
                                tAATGTACTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCTGTGATAGTTGGCCCGCTAGAAACCATGTCCATTTCCTGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGAGTTGATAATCTGGATGTTCGGCGTc                                                                                                                       <  2:4965318/150‑1 (MQ=255)
                                tAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTc                                                                                                                       <  2:2888270/150‑1 (MQ=255)
                                tAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTc                                                                                                                       >  1:1803140/1‑150 (MQ=255)
                                 aaTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCt                                                                                                                      >  2:3003080/1‑150 (MQ=255)
                                     aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCttgtt                                                                                                                  <  2:4320295/150‑1 (MQ=255)
                                      ccTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGTATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCttgttg                                                                                                                 >  1:4521155/1‑150 (MQ=255)
                                       cTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGa                                                                                                                >  2:614942/1‑150 (MQ=255)
                                          ctcTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACa                                                                                                             <  1:633305/148‑1 (MQ=255)
                                          cGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACa                                                                                                             <  1:1475309/150‑1 (MQ=255)
                                           gcgttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAg                                                                                                            <  2:2750258/147‑1 (MQ=255)
                                             tttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGgc                                                                                                          >  2:4899798/1‑150 (MQ=255)
                                                cGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCt                                                                                                       <  1:1238910/150‑1 (MQ=255)
                                                cGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCt                                                                                                       >  1:1683025/1‑150 (MQ=255)
                                                  aTGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGa                                                                                                     <  1:5105305/150‑1 (MQ=255)
                                                   tgtgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGAt                                                                                                    <  1:1080768/150‑1 (MQ=255)
                                                   tgtgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGAt                                                                                                    <  1:4829043/150‑1 (MQ=255)
                                                   tgtgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGAt                                                                                                    <  2:1403781/150‑1 (MQ=255)
                                                     tgAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAg                                                                                                  <  1:1885526/150‑1 (MQ=255)
                                                      gAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAgg                                                                                                 <  2:1023574/150‑1 (MQ=255)
                                                        aCTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGtt                                                                                               <  1:3053357/150‑1 (MQ=255)
                                                        aCTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGtt                                                                                               <  1:685248/150‑1 (MQ=255)
                                                           tGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCa                                                                                            <  2:4161827/150‑1 (MQ=255)
                                                             cTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACg                                                                                          >  1:871822/1‑150 (MQ=255)
                                                             cTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAAGCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGGTAGGTTTCACg                                                                                          >  2:4718676/1‑150 (MQ=255)
                                                              tCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCAGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCActt                                                                                         >  1:2640413/1‑148 (MQ=255)
                                                                    ggAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGa                                                                                   <  2:2681202/150‑1 (MQ=255)
                                                                      agagTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATg                                                                                 <  1:1534523/150‑1 (MQ=255)
                                                                      agagTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATg                                                                                 >  1:716160/1‑150 (MQ=255)
                                                                       gagTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGc                                                                                <  1:4860568/150‑1 (MQ=255)
                                                                        agTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCa                                                                               >  2:1212136/1‑150 (MQ=255)
                                                                        agTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAGGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCa                                                                               <  2:3790159/150‑1 (MQ=255)
                                                                         gtgtGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCAt                                                                              >  2:2072086/1‑150 (MQ=255)
                                                                           gtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCa                                                                            >  2:4554142/1‑150 (MQ=255)
                                                                            tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAc                                                                           <  2:1719119/150‑1 (MQ=255)
                                                                            tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCAc                                                                           <  2:4694245/150‑1 (MQ=255)
                                                                              gACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCg                                                                         <  1:2409205/150‑1 (MQ=255)
                                                                               aCCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCgg                                                                        <  1:1701751/150‑1 (MQ=255)
                                                                               aCCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCgg                                                                        <  1:3719583/150‑1 (MQ=255)
                                                                                 cGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGgag                                                                      <  1:1095226/150‑1 (MQ=255)
                                                                                  ggTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGgaga                                                                     >  1:2506001/1‑150 (MQ=255)
                                                                                  ggTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGgaga                                                                     >  2:181051/1‑150 (MQ=255)
                                                                                   gTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGaa                                                                    >  2:3767659/1‑150 (MQ=255)
                                                                                      aTAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTa                                                                 <  1:3338690/150‑1 (MQ=255)
                                                                                      aTAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTa                                                                 <  2:4848464/150‑1 (MQ=255)
                                                                                       tAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAg                                                                >  1:779379/1‑150 (MQ=255)
                                                                                       tAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAg                                                                >  2:4009027/1‑150 (MQ=255)
                                                                                          ttGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGt                                                             <  2:1480560/150‑1 (MQ=255)
                                                                                           tGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTc                                                            <  1:2698280/150‑1 (MQ=255)
                                                                                                cGATAGAACCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACTTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCt                                                       <  1:521768/150‑1 (MQ=255)
                                                                                                cGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCt                                                       <  1:2781503/150‑1 (MQ=255)
                                                                                                cGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCt                                                       <  1:1309482/150‑1 (MQ=255)
                                                                                                cGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCt                                                       >  1:4413183/1‑150 (MQ=255)
                                                                                                   tAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcc                                                    <  1:2840701/150‑1 (MQ=255)
                                                                                                   tAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcc                                                    >  1:3279369/1‑150 (MQ=255)
                                                                                                    aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcca                                                   <  1:344924/150‑1 (MQ=255)
                                                                                                    aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcca                                                   <  2:427485/150‑1 (MQ=255)
                                                                                                    aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgcca                                                   <  2:3126998/150‑1 (MQ=255)
                                                                                                       aaCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcc                                                >  2:2402777/1‑150 (MQ=255)
                                                                                                        aCCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTgccagcca                                               >  2:1377848/1‑150 (MQ=255)
                                                                                                             gTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGata                                          >  1:3862878/1‑150 (MQ=255)
                                                                                                              tCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGatat                                         <  2:3111472/150‑1 (MQ=255)
                                                                                                                  tttCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcgc                                     >  2:1213107/1‑150 (MQ=255)
                                                                                                                  tttCCGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATgcgc                                     >  1:4312100/1‑150 (MQ=255)
                                                                                                                     ccGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCtt                                  <  1:3202073/150‑1 (MQ=255)
                                                                                                                     ccGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCtt                                  >  1:4062890/1‑150 (MQ=255)
                                                                                                                      cGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTTGATGTTCGGCGTCTTGTTGAACAGGCGATGATATGTTTCACGTACAAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCttt                                 <  1:4281689/150‑1 (MQ=255)
                                                                                                                      cGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCttt                                 >  2:1844793/1‑150 (MQ=255)
                                                                                                                      cGGATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCttt                                 <  1:4009027/150‑1 (MQ=255)
                                                                                                                       ggATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTc                                >  2:304154/1‑150 (MQ=255)
                                                                                                                        gATACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCg                               >  1:1255442/1‑150 (MQ=255)
                                                                                                                         aTACGGTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGc                              >  1:2031942/1‑150 (MQ=255)
                                                                                                                             ggTTTTTTGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCg                          >  2:2274187/1‑150 (MQ=255)
                                                                                                                                  tttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGtttt                     >  1:4220642/1‑150 (MQ=255)
                                                                                                                                   ttGAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTc                    <  2:3300002/150‑1 (MQ=255)
                                                                                                                                     gAACAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTcgc                  >  2:3878484/1‑150 (MQ=255)
                                                                                                                                       aCAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTGTCGCTTCGGTTTTCgcgc                >  1:1440836/1‑150 (MQ=255)
                                                                                                                                        cAGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcc               >  2:3348067/1‑150 (MQ=255)
                                                                                                                                         aGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcca              <  1:2964392/150‑1 (MQ=255)
                                                                                                                                         aGACCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgcca              <  2:4413183/150‑1 (MQ=255)
                                                                                                                                           aCCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccagc            <  1:4788333/150‑1 (MQ=255)
                                                                                                                                           aCCACATTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCgccagc            >  1:2872191/1‑150 (MQ=255)
                                                                                                                                              acaTTCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGccagccag         >  2:3613786/1‑150 (MQ=255)
                                                                                                                                                aTTCCGGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTGGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGtt       <  1:3111844/150‑1 (MQ=255)
                                                                                                                                                 ttCCTGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGttt      >  2:4107309/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTACCAGATGCATCACCGGAGAATTAGCGTCCGGCTGCCAGCCAGGATATGCGCCTTTCGCTTCGGTTTTCGCGCCAGCCAGTTTACCC  >  CP009273/250760‑251058

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