Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,661,610 A→T 100% F313L (TTT→TTA yhjA ← putative cytochrome C peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,661,6100AT97.0% 739.9 / 9.3 230F313L (TTT→TTAyhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base T (109/114);  total (113/117)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 233 total aligned reads displayed.

CGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAAC  >  CP009273/3661465‑3661754
                                                                                                                                                 |                                                                                                                                                
cgccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCtaa                                                                                                                                                <  2:2304470/148‑1 (MQ=255)
 gccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAAt                                                                                                                                               <  1:2403608/148‑1 (MQ=255)
  ccgAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAAtt                                                                                                                                              <  1:280057/148‑1 (MQ=255)
    gagaATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCa                                                                                                                                            >  1:1395815/1‑148 (MQ=255)
       aataatACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCa                                                                                                                                         <  1:2117217/148‑1 (MQ=255)
        ataataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCaa                                                                                                                                        >  1:2855354/1‑148 (MQ=255)
        ataataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCaa                                                                                                                                        <  1:2296050/148‑1 (MQ=255)
          aataCCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATg                                                                                                                                      >  1:1799034/1‑148 (MQ=255)
             accaccATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCa                                                                                                                                   <  2:3412448/148‑1 (MQ=255)
              ccaccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCAt                                                                                                                                  <  1:555603/148‑1 (MQ=255)
              ccaccaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCAt                                                                                                                                  <  1:1409110/148‑1 (MQ=255)
                accaTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCa                                                                                                                                <  2:2558778/148‑1 (MQ=255)
                  caTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGt                                                                                                                              >  1:1499174/1‑148 (MQ=255)
                   aTGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTa                                                                                                                             <  1:761611/148‑1 (MQ=255)
                    tGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTACTTCAGCAATGGCAGCAGGaa                                                                                                                            >  2:2850329/1‑148 (MQ=255)
                    tGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTaa                                                                                                                            <  2:3543254/148‑1 (MQ=255)
                    tGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTATCTAATTCAGCAATGGCATCAGTaa                                                                                                                            >  2:2772138/1‑148 (MQ=255)
                      aCAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAata                                                                                                                          >  2:2227064/1‑148 (MQ=255)
                        aaGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATAtt                                                                                                                        >  2:928389/1‑148 (MQ=255)
                        aaGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGAAATGGCATCAGTAATAtt                                                                                                                        >  2:2706491/1‑148 (MQ=255)
                         aGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCACTGGCATCAGTAATAttt                                                                                                                       >  2:1932803/1‑148 (MQ=255)
                         aGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATAttt                                                                                                                       >  2:678582/1‑148 (MQ=255)
                          gTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATAtttt                                                                                                                      >  2:2737496/1‑148 (MQ=255)
                            tGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCg                                                                                                                    <  2:482080/148‑1 (MQ=255)
                            tGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCg                                                                                                                    >  2:1168207/1‑148 (MQ=255)
                             gCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAAACATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGc                                                                                                                   <  1:1321416/148‑1 (MQ=255)
                                caTTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCAc                                                                                                                <  2:2640987/148‑1 (MQ=255)
                                   ttattaTCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGa                                                                                                             >  2:312001/1‑148 (MQ=255)
                                       tatCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAAcc                                                                                                         <  1:613046/148‑1 (MQ=255)
                                        atCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCt                                                                                                        <  2:3161205/148‑1 (MQ=255)
                                        atCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCt                                                                                                        >  1:1026564/1‑148 (MQ=255)
                                         tCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCtt                                                                                                       <  2:2513697/148‑1 (MQ=255)
                                         tCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCtt                                                                                                       >  1:3164810/1‑148 (MQ=255)
                                         tCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCtt                                                                                                       >  2:3402479/1‑148 (MQ=255)
                                         tCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCtt                                                                                                       >  2:3402621/1‑148 (MQ=255)
                                          cTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTg                                                                                                      <  1:1323321/148‑1 (MQ=255)
                                                ataatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATa                                                                                                >  1:1227488/1‑148 (MQ=255)
                                                  aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTAGCCACTGAAACCTTGCGGATAGa                                                                                              >  2:2963032/1‑148 (MQ=255)
                                                  aatTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATAAATGTTTTCTCTAATTCAGCAATGGCATAAGTAATATAAACGCCACTGAAACCTTGCGGATAGa                                                                                              >  2:2963212/1‑148 (MQ=255)
                                                    ttGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACt                                                                                            >  2:3014578/1‑148 (MQ=255)
                                                        tAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGa                                                                                        >  1:3464623/1‑148 (MQ=255)
                                                        tAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGa                                                                                        <  1:2676379/148‑1 (MQ=255)
                                                         aGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAg                                                                                       >  1:483882/1‑148 (MQ=255)
                                                          gCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGa                                                                                              <  2:1161299/140‑1 (MQ=255)
                                                            cTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGaa                                                                                    <  1:911043/148‑1 (MQ=255)
                                                             ttttttCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAAc                                                                                   >  1:2947568/1‑148 (MQ=255)
                                                               ttgtCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTg                                                                                 <  2:1961410/145‑1 (MQ=255)
                                                                 ttCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCg                                                                               <  2:1580155/148‑1 (MQ=255)
                                                                  tCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGt                                                                              <  1:2445861/148‑1 (MQ=255)
                                                                    tGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGttt                                                                            >  1:1134477/1‑148 (MQ=255)
                                                                     gTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGtttt                                                                           <  2:1778665/148‑1 (MQ=255)
                                                                     gTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGtttt                                                                           >  2:326678/1‑148 (MQ=255)
                                                                       tGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTaa                                                                         >  2:86232/1‑148 (MQ=255)
                                                                         cgcCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGc                                                                       >  1:3372522/1‑148 (MQ=255)
                                                                           ccGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTg                                                                     <  1:452047/148‑1 (MQ=255)
                                                                           ccGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTg                                                                     >  1:3520046/1‑148 (MQ=255)
                                                                             gTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCg                                                                   <  2:285232/148‑1 (MQ=255)
                                                                              tCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCgg                                                                  <  1:184596/148‑1 (MQ=255)
                                                                               cAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGa                                                                 <  2:1396844/148‑1 (MQ=255)
                                                                                aGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGAt                                                                <  1:2076558/148‑1 (MQ=255)
                                                                                 gcgcATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCTTAATAAAAGTTTTCTCTAATTCAGCAATGGCATCAGTACTATTTTCGCCACTGAACCCTTGCGGATAGACTTCGTGGACCTGCGTTTTAAGCTGCGGATc                                                               >  2:418228/1‑148 (MQ=255)
                                                                                  cgcATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCt                                                              <  1:678582/148‑1 (MQ=255)
                                                                                        ttCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTcc                                                        >  1:272114/1‑148 (MQ=255)
                                                                                           aTCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGc                                                     <  2:3232505/148‑1 (MQ=255)
                                                                                             ctcCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCtt                                                   <  1:2018985/148‑1 (MQ=255)
                                                                                             ctcCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCtt                                                   <  2:1682416/148‑1 (MQ=255)
                                                                                             ctcCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCtt                                                   <  1:1542313/148‑1 (MQ=255)
                                                                                                cgcgCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGc                                                <  1:837529/148‑1 (MQ=255)
                                                                                                  cgcAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCaa                                              >  1:1101000/1‑148 (MQ=255)
                                                                                                    cAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCaata                                            <  1:326678/148‑1 (MQ=255)
                                                                                                      aCCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCaataat                                          <  1:1716174/148‑1 (MQ=255)
                                                                                                      aCCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCaataat                                          <  2:1049209/148‑1 (MQ=255)
                                                                                                       ccATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAAtt                                         >  2:3223099/1‑148 (MQ=255)
                                                                                                        cATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAAttt                                        >  2:2195727/1‑148 (MQ=255)
                                                                                                        cATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAAttt                                        >  1:188629/1‑148 (MQ=255)
                                                                                                         aTTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTc                                       >  2:581689/1‑148 (MQ=255)
                                                                                                         aTTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTc                                       <  1:398463/148‑1 (MQ=255)
                                                                                                             aTCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGTCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTcc                                   <  2:853095/148‑1 (MQ=255)
                                                                                                             aTCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTcc                                   <  2:1654071/148‑1 (MQ=255)
                                                                                                             aTCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTcc                                   <  2:1239641/148‑1 (MQ=255)
                                                                                                                aaaTGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAg                                <  1:2512065/148‑1 (MQ=255)
                                                                                                                   tGGGGAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGAt                             >  2:2283113/1‑148 (MQ=255)
                                                                                                                    ggggAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGAtt                            <  2:2376309/148‑1 (MQ=255)
                                                                                                                      ggAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTc                          >  1:3483419/1‑148 (MQ=255)
                                                                                                                      ggAATCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTc                          <  2:3310828/148‑1 (MQ=255)
                                                                                                                        aaTCCGGCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGa                        >  2:1011821/1‑148 (MQ=255)
                                                                                                                             ggCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCa                   >  1:2539039/1‑148 (MQ=255)
                                                                                                                             ggCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCa                   <  2:2447977/148‑1 (MQ=255)
                                                                                                                              gCGTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCAt                  <  2:38495/148‑1 (MQ=255)
                                                                                                                                gTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCAttt                <  1:3402479/148‑1 (MQ=255)
                                                                                                                                gTAATTAATGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCAttt                <  1:3402621/148‑1 (MQ=255)
                                                                                                                                    ttaatGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAAt            >  2:2868666/1‑148 (MQ=255)
                                                                                                                                      aatGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCg          <  1:1168087/148‑1 (MQ=255)
                                                                                                                                       atGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCgg         <  1:312001/148‑1 (MQ=255)
                                                                                                                                        tGTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCggg        <  1:1446804/148‑1 (MQ=255)
                                                                                                                                         gTTTTCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGt       >  2:1270547/1‑148 (MQ=255)
                                                                                                                                           tttCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTc     >  1:3025325/1‑148 (MQ=255)
                                                                                                                                            ttCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCa    >  1:63548/1‑148 (MQ=255)
                                                                                                                                            ttCTCTAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCa    <  1:2368098/148‑1 (MQ=255)
                                                                                                                                              ctctAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAAc  >  2:1337611/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCCGAGAATAATACCACCATGACAAGTTGCACATTTATTATCTTTAAATAATTGATAGCCTTTTTTCTGTTGCGCCGTCAGCGCATTTTCATCTCCGCGCAACCATTTATCAAATGGGGAATCCGGCGTAATTAATGTTTTCTCAAATTCAGCAATGGCATCAGTAATATTTTCGCCACTGAAACCTTGCGGATAGACTTCGAGGAACTGCGTTTTAAGCTGCGGATCTTTTTCCAGCTTAGCAATAATTTCGTCCCAGGATTTCGACGCCATTTCAATCGGGTTCAAC  >  CP009273/3661465‑3661754

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: