Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,312,055 C→A 100% G20G (GGC→GGA)  argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,312,0550CA97.6% 793.1 / 8.0 249G20G (GGC→GGA) argGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (137/106);  total (141/108)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 249 total aligned reads displayed.

TAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGCAC  >  CP009273/3311908‑3312200
                                                                                                                                                   |                                                                                                                                                 
tAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAgg                                                                                                                                                 >  1:3294141/1‑150 (MQ=255)
 aaaTGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGt                                                                                                                                                <  1:4637261/150‑1 (MQ=255)
   aTGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCt                                                                                                                                              >  1:4992443/1‑150 (MQ=255)
     gAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTgg                                                                                                                                            <  2:2919751/150‑1 (MQ=255)
       aaaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGac                                                                                                                                          >  1:4112202/1‑150 (MQ=255)
        aaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGaca                                                                                                                                         <  1:2029887/150‑1 (MQ=255)
           tCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCa                                                                                                                                      >  1:2547496/1‑150 (MQ=255)
            cATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAg                                                                                                                                     <  2:2502374/150‑1 (MQ=255)
              tttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTg                                                                                                                                   >  1:1122348/1‑150 (MQ=255)
              tttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTg                                                                                                                                   >  1:3177489/1‑150 (MQ=255)
                    ttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCAc                                                                                                                             >  1:1937157/1‑150 (MQ=255)
                      gCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCTGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActg                                                                                                                           >  2:1736894/1‑150 (MQ=255)
                       cATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActgc                                                                                                                          <  1:3092141/150‑1 (MQ=255)
                        aTAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActgct                                                                                                                         >  1:1715018/1‑150 (MQ=255)
                          aaaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCtgt                                                                                                                       >  2:218632/1‑150 (MQ=255)
                           aaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCtgtg                                                                                                                      <  2:3643947/150‑1 (MQ=255)
                           aaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCtgtg                                                                                                                      <  1:3241264/150‑1 (MQ=255)
                            aaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTgg                                                                                                                     >  2:2752854/1‑150 (MQ=255)
                              aTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGAt                                                                                                                   >  2:1044057/1‑150 (MQ=255)
                               ttCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATg                                                                                                                  <  1:2552281/150‑1 (MQ=255)
                                tCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGc                                                                                                                 <  2:1142400/150‑1 (MQ=255)
                                 cAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCg                                                                                                                >  1:745844/1‑150 (MQ=255)
                                    tGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTGTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACa                                                                                                             >  1:3379792/1‑150 (MQ=255)
                                      agagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaa                                                                                                           <  1:3051296/150‑1 (MQ=255)
                                      agagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaa                                                                                                           >  2:130610/1‑150 (MQ=255)
                                        agCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAg                                                                                                         >  1:2810248/1‑150 (MQ=255)
                                          cGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTTTTATTTAATGACGACGATTCTCAAGCATCTCCAGGAAGTTCACCTTAGTGGTATCGCTTATTCTGGAGGTCTGTACACCAGTGCCGCACTGCTGTGGATGCGACAAAAggg                                                                                                       >  1:68321/1‑150 (MQ=255)
                                          cGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAggg                                                                                                       >  2:1589815/1‑150 (MQ=255)
                                           ggAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGa                                                                                                      >  1:3752500/1‑150 (MQ=255)
                                           ggAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGa                                                                                                      >  2:3297842/1‑150 (MQ=255)
                                           ggAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGa                                                                                                      >  2:3243379/1‑150 (MQ=255)
                                             aaaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGc                                                                                                    >  2:346140/1‑150 (MQ=255)
                                             aaaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGc                                                                                                    <  2:5135438/150‑1 (MQ=255)
                                             aaaTCCAGGCTCATCATCAGTTAATTAAGAAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGAATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTTCCGCACTGCTTTGGATGCGACAAAAGGGAGc                                                                                                    >  2:52586/1‑150 (MQ=255)
                                              aaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCg                                                                                                   >  1:3622440/1‑150 (MQ=255)
                                              aaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCg                                                                                                   <  1:966709/150‑1 (MQ=255)
                                               aTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCgg                                                                                                  <  2:2470376/150‑1 (MQ=255)
                                                   aGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTcc                                                                                              <  1:1587075/150‑1 (MQ=255)
                                                     gCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCtt                                                                                            >  1:2108937/1‑150 (MQ=255)
                                                       tcatcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTAt                                                                                          <  2:4070396/150‑1 (MQ=255)
                                                        catcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATg                                                                                         <  1:3937394/150‑1 (MQ=255)
                                                           catcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCat                                                                                      >  1:663368/1‑150 (MQ=255)
                                                              caGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCAtata                                                                                   >  2:236603/1‑150 (MQ=255)
                                                                    attaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCaa                                                                             >  2:1128875/1‑150 (MQ=255)
                                                                       aaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAAcc                                                                          >  1:775916/1‑150 (MQ=255)
                                                                          cAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTgg                                                                       <  2:1611833/150‑1 (MQ=255)
                                                                          cAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTgg                                                                       >  1:3456039/1‑150 (MQ=255)
                                                                          cAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTgg                                                                       <  1:2053484/150‑1 (MQ=255)
                                                                           aGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTggg                                                                      <  2:1963597/150‑1 (MQ=255)
                                                                           aGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTggg                                                                      <  1:5042721/150‑1 (MQ=255)
                                                                            gggTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTCTTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAATGGGAGCGGTTCCTTCTGCATATACTGCAAACCTGGgc                                                                     >  2:4784391/1‑150 (MQ=255)
                                                                               tgtTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccag                                                                  >  2:4094566/1‑150 (MQ=255)
                                                                               tgtTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccag                                                                  >  2:306718/1‑150 (MQ=255)
                                                                                  tATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagcca                                                               <  2:609475/150‑1 (MQ=255)
                                                                                   aTTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGccagccag                                                              <  1:233978/150‑1 (MQ=255)
                                                                                     tttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGAc                                                            >  1:2199180/1‑150 (MQ=255)
                                                                                     tttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGAc                                                            <  1:518429/150‑1 (MQ=255)
                                                                                       tATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGa                                                          <  1:1024099/150‑1 (MQ=255)
                                                                                        aTGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGaa                                                         >  1:5034575/1‑150 (MQ=255)
                                                                                           acgacgATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAagag                                                      >  1:3967362/1‑150 (MQ=255)
                                                                                            cgacgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                     <  2:2589765/150‑1 (MQ=255)
                                                                                            cgacgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                     <  2:863741/150‑1 (MQ=255)
                                                                                              acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAt                                                   >  2:1082373/1‑150 (MQ=255)
                                                                                              acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAt                                                   <  2:2370033/150‑1 (MQ=255)
                                                                                              acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAt                                                   >  1:2974355/1‑150 (MQ=255)
                                                                                               cgaTTCTCAAGCATCTCCCTGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAtt                                                  <  2:882948/150‑1 (MQ=255)
                                                                                               cgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAtt                                                  >  1:3775184/1‑150 (MQ=255)
                                                                                                gaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTa                                                 >  1:1288342/1‑150 (MQ=255)
                                                                                                 aTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTat                                                >  1:4710692/1‑150 (MQ=255)
                                                                                                  ttCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTatg                                               >  2:1540462/1‑150 (MQ=255)
                                                                                                   tctcAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTatga                                              >  1:344810/1‑150 (MQ=255)
                                                                                                       aaGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCg                                          >  1:1609121/1‑150 (MQ=255)
                                                                                                       aaGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCg                                          <  1:2844522/150‑1 (MQ=255)
                                                                                                          cATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATc                                       >  2:3475978/1‑150 (MQ=255)
                                                                                                              tcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtc                                   >  1:1799747/1‑150 (MQ=255)
                                                                                                              tcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtc                                   <  1:4549896/150‑1 (MQ=255)
                                                                                                               cccGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcg                                  >  1:4870625/1‑150 (MQ=255)
                                                                                                                 cGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtc                                <  1:171963/150‑1 (MQ=255)
                                                                                                                 cGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtc                                <  2:4292697/150‑1 (MQ=255)
                                                                                                                   gTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTcgtcgt                              >  2:2441215/1‑150 (MQ=255)
                                                                                                                    tAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTg                             <  2:2338/150‑1 (MQ=255)
                                                                                                                         cAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATg                        <  2:2167674/150‑1 (MQ=255)
                                                                                                                           aCGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGa                      >  2:4442331/1‑150 (MQ=255)
                                                                                                                            cGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGaa                     <  1:3860756/150‑1 (MQ=255)
                                                                                                                            cGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGaa                     <  1:2649876/150‑1 (MQ=255)
                                                                                                                             gTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAAt                    >  1:87929/1‑150 (MQ=255)
                                                                                                                              tATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGTGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                   >  2:67025/1‑150 (MQ=255)
                                                                                                                              tATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                   >  1:135972/1‑150 (MQ=255)
                                                                                                                               aTTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATAc                  <  2:1410646/150‑1 (MQ=255)
                                                                                                                                  ggTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGgc               >  1:2269864/1‑150 (MQ=255)
                                                                                                                                   gTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATTCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGgcg              <  2:2702790/150‑1 (MQ=255)
                                                                                                                                   gTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAACACGgcg              <  2:3643681/150‑1 (MQ=255)
                                                                                                                                    tATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGgcgc             <  1:4779153/150‑1 (MQ=255)
                                                                                                                                     aTCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGcgcg            >  1:4204332/1‑150 (MQ=255)
                                                                                                                                       cGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGga          <  2:3967362/150‑1 (MQ=255)
                                                                                                                                        gCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGgag         >  1:1059631/1‑150 (MQ=255)
                                                                                                                                         cTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGgaga        <  1:3163415/150‑1 (MQ=255)
                                                                                                                                         cTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGgaga        <  2:1715018/150‑1 (MQ=255)
                                                                                                                                           ttttCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAAc      <  2:2547496/150‑1 (MQ=255)
                                                                                                                                               cTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGCAc  >  1:2771834/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TAAATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGCAC  >  CP009273/3311908‑3312200

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