Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,615,947 G→T 100% E73D (GAG→GAT radA → DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,615,9470GT97.5% 627.8 / 4.5 202E73D (GAG→GATradADNA repair protein
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (100/97);  total (103/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 203 total aligned reads displayed.

CAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAG  >  CP009273/4615803‑4616091
                                                                                                                                                |                                                                                                                                                
cagtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTg                                                                                                                                               <  2:2453204/148‑1 (MQ=255)
  gtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGcc                                                                                                                                             >  2:2310806/1‑148 (MQ=255)
     ccTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcgcg                                                                                                                                          >  2:3408003/1‑148 (MQ=255)
      cTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCGGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGt                                                                                                                                         <  2:220273/148‑1 (MQ=255)
        gTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttt                                                                                                                                       <  1:251571/148‑1 (MQ=255)
        gTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttt                                                                                                                                       >  1:1999110/1‑148 (MQ=255)
         tCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGtttt                                                                                                                                      >  1:2063709/1‑148 (MQ=255)
          cATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAATTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttttt                                                                                                                                     <  2:529874/148‑1 (MQ=255)
          cATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttttt                                                                                                                                     <  2:1038123/148‑1 (MQ=255)
          cATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAAGGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttttt                                                                                                                                     <  2:81484/148‑1 (MQ=255)
             gCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCa                                                                                                                                  <  2:2680883/148‑1 (MQ=255)
             gCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCa                                                                                                                                  >  1:1813747/1‑148 (MQ=255)
                 ggAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCgg                                                                                                                              >  2:2234870/1‑148 (MQ=255)
                 ggAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCgg                                                                                                                              <  1:545341/148‑1 (MQ=255)
                    acacCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGAtt                                                                                                                           >  1:970569/1‑148 (MQ=255)
                       ccATCACCGCGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTAGGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                        <  2:914721/148‑1 (MQ=255)
                        cATCACCGAGGTGCGTCTTGCAGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaaa                                                                                                                       >  2:84511/1‑148 (MQ=255)
                         aTCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGACAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAag                                                                                                                      >  1:2231857/1‑148 (MQ=255)
                           cACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                    >  2:2130162/1‑148 (MQ=255)
                           cACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                    >  2:1715538/1‑148 (MQ=255)
                            aCCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                   >  1:3507031/1‑148 (MQ=255)
                            aCCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                   <  1:1480292/148‑1 (MQ=255)
                            aCCGAGGTGCGTCTTGCTGCGTCGCCAATGGGGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                   <  2:1673245/148‑1 (MQ=255)
                                   tGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACcg                                                                                                            >  2:1620213/1‑148 (MQ=255)
                                   tGCGTCTTGCTGCGTCACCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACcg                                                                                                            >  1:1096807/1‑148 (MQ=255)
                                    gCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACcgc                                                                                                           >  1:1637720/1‑148 (MQ=255)
                                       tCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGTGTCTATGCCGGGAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTa                                                                                                        <  2:2598749/148‑1 (MQ=255)
                                       tCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTa                                                                                                        >  1:1740045/1‑148 (MQ=255)
                                            ctgcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAgg                                                                                                   >  1:1191596/1‑148 (MQ=255)
                                                        tggtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtggt                                                                                       >  1:1207807/1‑148 (MQ=255)
                                                          gtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGc                                                                                     <  1:148211/148‑1 (MQ=255)
                                                             gcgcGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAg                                                                                  >  1:2659609/1‑148 (MQ=255)
                                                              cgcgTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAgg                                                                                 <  1:290570/148‑1 (MQ=255)
                                                               gcgTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGa                                                                                <  2:1904190/148‑1 (MQ=255)
                                                                  tAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGt                                                                             >  1:420860/1‑148 (MQ=255)
                                                                    aCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGc                                                                           >  2:851927/1‑148 (MQ=255)
                                                                      gAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCa                                                                         >  1:3019602/1‑148 (MQ=255)
                                                                       aGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAt                                                                        <  2:2332546/148‑1 (MQ=255)
                                                                        gCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                       <  2:2860317/148‑1 (MQ=255)
                                                                        gCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCCTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                       <  1:1537189/148‑1 (MQ=255)
                                                                         cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                      <  1:3393755/148‑1 (MQ=255)
                                                                         cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                      >  1:3330825/1‑148 (MQ=255)
                                                                         cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                      >  1:1267145/1‑148 (MQ=255)
                                                                            ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                   >  1:3487653/1‑148 (MQ=255)
                                                                            ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                   <  1:1047354/148‑1 (MQ=255)
                                                                            ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                   <  1:1364089/148‑1 (MQ=255)
                                                                               aGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                                                                <  2:1769157/148‑1 (MQ=255)
                                                                                gCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTgg                                                               >  2:610443/1‑148 (MQ=255)
                                                                                gCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTgg                                                               <  2:3337097/148‑1 (MQ=255)
                                                                                 cGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGc                                                              >  2:1662361/1‑148 (MQ=255)
                                                                                   gCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCgg                                                            <  1:703814/148‑1 (MQ=255)
                                                                                       tGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAc                                                        >  1:2966054/1‑148 (MQ=255)
                                                                                       tGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAc                                                        >  2:3045083/1‑148 (MQ=255)
                                                                                          cGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt                                                     <  1:3242787/148‑1 (MQ=255)
                                                                                             tAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGt                                                  >  1:2560871/1‑148 (MQ=255)
                                                                                             tAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGt                                                  <  2:1593383/148‑1 (MQ=255)
                                                                                               gcgcCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGc                                                >  2:163378/1‑148 (MQ=255)
                                                                                                cgcCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCg                                               <  1:1675513/148‑1 (MQ=255)
                                                                                                cgcCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCg                                               >  1:196533/1‑148 (MQ=255)
                                                                                                    ggggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGa                                           <  1:47938/148‑1 (MQ=255)
                                                                                                     gggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaa                                          <  1:1026646/148‑1 (MQ=255)
                                                                                                         ggCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATcc                                      <  2:767829/148‑1 (MQ=255)
                                                                                                           cAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCAc                                    >  2:1369863/1‑148 (MQ=255)
                                                                                                             aaaGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgc                                  >  2:3446721/1‑148 (MQ=255)
                                                                                                             aaaGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgc                                  <  1:814736/148‑1 (MQ=255)
                                                                                                               actCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctg                                <  1:1825642/146‑1 (MQ=255)
                                                                                                               aGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAACTCCACgctg                                >  2:1198544/1‑148 (MQ=255)
                                                                                                               aGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctg                                <  1:2801825/148‑1 (MQ=255)
                                                                                                               aGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctg                                <  1:2801154/148‑1 (MQ=255)
                                                                                                                 tCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgct                              >  2:977740/1‑148 (MQ=255)
                                                                                                                  ccAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTa                             >  2:3134455/1‑148 (MQ=255)
                                                                                                                    aGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACt                           >  2:289486/1‑148 (MQ=255)
                                                                                                                      aaaCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGc                         >  1:641257/1‑148 (MQ=255)
                                                                                                                       aaCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCa                        >  1:561802/1‑148 (MQ=255)
                                                                                                                         ctctCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCaaa                      >  1:921918/1‑148 (MQ=255)
                                                                                                                          tctcCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAAc                     <  2:2058373/148‑1 (MQ=255)
                                                                                                                            gccGATATCAGCCTTGAAGATCTGCCGCGTTGTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGc                   <  1:2700472/147‑1 (MQ=255)
                                                                                                                             ccGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCt                  <  1:2603579/148‑1 (MQ=255)
                                                                                                                             ccGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCt                  >  1:2343250/1‑148 (MQ=255)
                                                                                                                                atatCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCtgtg               >  1:1037810/1‑148 (MQ=255)
                                                                                                                                atatCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCtgtg               >  2:2432376/1‑148 (MQ=255)
                                                                                                                                  atCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCa             >  1:2691211/1‑148 (MQ=255)
                                                                                                                                   tCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaa            >  2:1890749/1‑148 (MQ=255)
                                                                                                                                   tCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaa            <  2:556747/148‑1 (MQ=255)
                                                                                                                                   tCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaa            <  2:552341/148‑1 (MQ=255)
                                                                                                                                    cAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaaa           >  2:3367598/1‑148 (MQ=255)
                                                                                                                                    cAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaaa           >  2:1628366/1‑148 (MQ=255)
                                                                                                                                     aGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATACACGCTGCTACTGCAAACGCTGTGCAAAc          <  1:851927/148‑1 (MQ=255)
                                                                                                                                      gCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACt         >  1:45962/1‑148 (MQ=255)
                                                                                                                                       ccTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg        <  1:1620213/148‑1 (MQ=255)
                                                                                                                                       ccTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg        >  2:1494440/1‑148 (MQ=255)
                                                                                                                                       ccTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg        <  2:2819640/148‑1 (MQ=255)
                                                                                                                                       ccTGGAAGAGCTGCCGCGTTTGTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg        <  2:3250976/148‑1 (MQ=255)
                                                                                                                                         ttGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGc      >  1:912578/1‑148 (MQ=255)
                                                                                                                                          tGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc     <  1:632416/148‑1 (MQ=255)
                                                                                                                                          tGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc     <  2:2427739/148‑1 (MQ=255)
                                                                                                                                          tGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc     <  1:2394959/148‑1 (MQ=255)
                                                                                                                                            aaGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCca   >  2:510528/1‑148 (MQ=255)
                                                                                                                                            aaGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCca   <  2:3218366/148‑1 (MQ=255)
                                                                                                                                             agatCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCCGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCcaa  >  1:754016/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAG  >  CP009273/4615803‑4616091

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: