Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,661,524 C→T 7.6% G342D (GGC→GAC)  yhjA ← putative cytochrome C peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,661,5240CT 7.6% 899.2 / 43.6 263G342D (GGC→GAC) yhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base C (121/121);  new base T (12/8);  total (134/129)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 265 total aligned reads displayed.

CTTTAGTCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAA  >  CP009273/3661375‑3661672
                                                                                                                                                     |                                                                                                                                                    
  ttAGTCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCt                                                                                                                                                    <  2:1407044/150‑1 (MQ=255)
  ttAGTCACAGTCATACGCCCAATAGCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCt                                                                                                                                                    <  1:3325312/150‑1 (MQ=255)
   tAGTCACATTCATACGACCAATAGCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCtt                                                                                                                                                   <  2:3082258/150‑1 (MQ=255)
     gTCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCtttt                                                                                                                                                 <  2:2340047/150‑1 (MQ=255)
      tCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCttttt                                                                                                                                                <  1:4303444/150‑1 (MQ=255)
      tCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGACttttt                                                                                                                                                >  2:4203109/1‑150 (MQ=255)
      tCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACAACCATGACAAGTTGCACATTTATTATCTTTAAATCATTGATAGCCtgttt                                                                                                                                                >  2:3890609/1‑150 (MQ=255)
       cacaTTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCtttttt                                                                                                                                               <  1:4012494/150‑1 (MQ=255)
        acaTTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTc                                                                                                                                              <  2:225781/150‑1 (MQ=255)
           ttCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGt                                                                                                                                           <  1:2766673/150‑1 (MQ=255)
               tACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTgcg                                                                                                                                       >  1:2709244/1‑150 (MQ=255)
               tACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATATATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTgcg                                                                                                                                       >  2:1876560/1‑150 (MQ=255)
                 cGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGcc                                                                                                                                     <  1:665806/150‑1 (MQ=255)
                  gACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCg                                                                                                                                    <  2:3839740/150‑1 (MQ=255)
                       atatCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAgc                                                                                                                               >  2:3146585/1‑150 (MQ=255)
                       atatCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAgc                                                                                                                               >  1:1042445/1‑150 (MQ=255)
                        tatCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAgcg                                                                                                                              >  1:2791050/1‑150 (MQ=255)
                        tatCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAgcg                                                                                                                              <  1:1326424/150‑1 (MQ=255)
                              ccgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtttt                                                                                                                        <  2:1762529/150‑1 (MQ=255)
                              ccgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACAACCATGACAAGTTGCACATTTATTATCTTTAAATAATAGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCAtttt                                                                                                                        >  1:2531636/1‑150 (MQ=255)
                               cgccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTc                                                                                                                       >  2:4351129/1‑150 (MQ=255)
                                gccgTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCa                                                                                                                      >  2:804328/1‑150 (MQ=255)
                                 ccgTTATTTCCCCATATTTAAAGTCTTTTTTCATCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCAt                                                                                                                     <  2:1368673/150‑1 (MQ=255)
                                     aaTTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTcc                                                                                                                 <  1:1960691/150‑1 (MQ=255)
                                         tCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCgcgc                                                                                                             >  1:151224/1‑150 (MQ=255)
                                         tCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCgcgc                                                                                                             >  2:1963783/1‑150 (MQ=255)
                                         tCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCgcgc                                                                                                             >  1:152474/1‑150 (MQ=255)
                                          ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                            >  1:970720/1‑150 (MQ=255)
                                          ccccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCa                                                                                                            <  1:1067824/150‑1 (MQ=255)
                                            ccAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAAc                                                                                                          <  1:1067940/150‑1 (MQ=255)
                                             cAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTAACTTTAAATAATTGATAGCCTATGTTCTGTTGCGCCGACAGCGCATTATCATCTCCGCGCAAcc                                                                                                         >  2:1547581/1‑150 (MQ=255)
                                                 gTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCAttt                                                                                                     <  1:4295974/150‑1 (MQ=255)
                                                  ttAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTa                                                                                                    >  2:2538045/1‑150 (MQ=255)
                                                  ttAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTa                                                                                                    >  2:3790277/1‑150 (MQ=255)
                                                    aaaGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATc                                                                                                  <  1:1162427/150‑1 (MQ=255)
                                                      aGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCaa                                                                                                >  1:1035576/1‑150 (MQ=255)
                                                      aGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCaa                                                                                                >  2:2186421/1‑150 (MQ=255)
                                                              ttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGaa                                                                                        >  2:3234758/1‑150 (MQ=255)
                                                              ttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGaa                                                                                        >  1:1957531/1‑150 (MQ=255)
                                                              ttCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGaa                                                                                        >  1:2612905/1‑150 (MQ=255)
                                                               tCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGact                                                                                       >  1:485614/1‑148 (MQ=255)
                                                                 aGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATTATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATcc                                                                                     <  2:2369576/150‑1 (MQ=255)
                                                                 aGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCCTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATcc                                                                                     <  2:192625/150‑1 (MQ=255)
                                                                 aGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATcc                                                                                     >  2:4206613/1‑150 (MQ=255)
                                                                  gCCCCCACCGTTCAAACGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACAGTTAATATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCg                                                                                    <  2:575498/150‑1 (MQ=255)
                                                                   ccccAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCgg                                                                                   >  2:2157684/1‑150 (MQ=255)
                                                                      cAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGt                                                                                >  1:1628449/1‑150 (MQ=255)
                                                                       accGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGta                                                                               <  1:3341128/148‑1 (MQ=255)
                                                                       aaCGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGta                                                                               >  1:3754056/1‑150 (MQ=255)
                                                                       aaCGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGta                                                                               >  2:2927747/1‑150 (MQ=255)
                                                                          ggTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaatt                                                                            >  2:2627710/1‑150 (MQ=255)
                                                                          ggTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaatt                                                                            >  2:3450703/1‑150 (MQ=255)
                                                                          ggTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaatt                                                                            <  1:2370945/150‑1 (MQ=255)
                                                                            ttCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGtaattaa                                                                          <  1:3967566/150‑1 (MQ=255)
                                                                               aaaGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGt                                                                       >  1:3204060/1‑150 (MQ=255)
                                                                               aaaGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGt                                                                       >  1:831501/1‑150 (MQ=255)
                                                                                 aGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCTCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGttt                                                                     <  2:3635913/150‑1 (MQ=255)
                                                                                  ggAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGtttt                                                                    <  1:2321081/150‑1 (MQ=255)
                                                                                    aaCGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtct                                                                  >  1:3914094/1‑150 (MQ=255)
                                                                                     aCGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtctc                                                                 <  1:1825721/150‑1 (MQ=255)
                                                                                     aCGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTtctc                                                                 <  2:1447419/150‑1 (MQ=255)
                                                                                          cgccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAAtt                                                            >  2:1183685/1‑150 (MQ=255)
                                                                                            ccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCa                                                          >  2:964495/1‑150 (MQ=255)
                                                                                               agaATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCa                                                       <  2:2323031/150‑1 (MQ=255)
                                                                                                gaATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCaa                                                      <  1:2888103/150‑1 (MQ=255)
                                                                                                gaATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCaa                                                      >  2:2704736/1‑150 (MQ=255)
                                                                                                 aataatACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAAt                                                     <  2:3422525/150‑1 (MQ=255)
                                                                                                  ataataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATg                                                    <  1:3835576/150‑1 (MQ=255)
                                                                                                     ataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCa                                                 >  1:534726/1‑150 (MQ=255)
                                                                                                     ataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCa                                                 >  1:534572/1‑150 (MQ=255)
                                                                                                        ccaccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCa                                              >  1:1328067/1‑150 (MQ=255)
                                                                                                          accaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGt                                            >  2:3669488/1‑150 (MQ=255)
                                                                                                            caTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTaa                                          >  2:4142121/1‑150 (MQ=255)
                                                                                                              tGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAata                                        >  2:704770/1‑150 (MQ=255)
                                                                                                                aCAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATAtt                                      >  1:3985102/1‑150 (MQ=255)
                                                                                                                aCAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATAtt                                      >  1:21678/1‑150 (MQ=255)
                                                                                                                  aaGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATAtttt                                    <  1:862542/150‑1 (MQ=255)
                                                                                                                  aaGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATAtttt                                    <  1:352602/150‑1 (MQ=255)
                                                                                                                    gTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCg                                  <  1:3047100/150‑1 (MQ=255)
                                                                                                                      tGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGGAATATTTTCGcc                                >  2:2234983/1‑150 (MQ=255)
                                                                                                                       gCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCa                               <  1:3244213/150‑1 (MQ=255)
                                                                                                                        cacaTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCAc                              >  2:3564762/1‑150 (MQ=255)
                                                                                                                        cacaTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATGTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCAc                              <  1:888844/150‑1 (MQ=255)
                                                                                                                           aTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGa                           <  1:1691082/150‑1 (MQ=255)
                                                                                                                            tttATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGaa                          <  2:3905673/150‑1 (MQ=255)
                                                                                                                              tattatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAAc                        >  1:1289304/1‑150 (MQ=255)
                                                                                                                               attatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGAATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAAcc                       <  1:1876560/150‑1 (MQ=255)
                                                                                                                                 tatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCtt                     >  2:3547681/1‑150 (MQ=255)
                                                                                                                                  atCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTg                    >  2:1377602/1‑150 (MQ=255)
                                                                                                                                    cTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCg                  <  2:1610703/150‑1 (MQ=255)
                                                                                                                                       tAAATAATTGATAGTCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGAt               >  2:2069238/1‑150 (MQ=255)
                                                                                                                                        aaaTAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATa              >  1:1206116/1‑150 (MQ=255)
                                                                                                                                          ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGa            >  2:2727683/1‑150 (MQ=255)
                                                                                                                                            aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt          <  1:1739821/150‑1 (MQ=255)
                                                                                                                                            aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt          <  1:906215/150‑1 (MQ=255)
                                                                                                                                              ttGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTc        <  1:3521507/150‑1 (MQ=255)
                                                                                                                                               tGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCg       >  1:1181358/1‑150 (MQ=255)
                                                                                                                                                gATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGa      <  2:865553/150‑1 (MQ=255)
                                                                                                                                                   aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGa   <  2:4018316/150‑1 (MQ=255)
                                                                                                                                                   aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGa   >  1:3101033/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTTAGTCACATTCATACGACCAATATCCGCCGCCGTAATTTCCCCAAAGTTAAAGTCTTTTTTCAGCCCCAACGGTTCAAAGGAACGTCCGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAA  >  CP009273/3661375‑3661672

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: