Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,318,228 C→A 100% S251* (TCA→TAA)  yciV → PHP domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,318,2280CA99.4% 574.8 / ‑4.0 174S251* (TCA→TAA) yciVPHP domain protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (83/90);  total (84/90)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 177 total aligned reads displayed.

AATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTT  >  CP009273/1318083‑1318370
                                                                                                                                                 |                                                                                                                                              
aaTCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATaac                                                                                                                                              >  1:2666999/1‑148 (MQ=255)
 aTCTTAGTGCTAAATGGCTGAAAAGATTGGGTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATaaca                                                                                                                                             <  2:2469699/148‑1 (MQ=255)
   cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                           >  1:2502477/1‑148 (MQ=255)
   cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                           >  2:5697/1‑148 (MQ=255)
   cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                           >  1:403445/1‑148 (MQ=255)
   cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                           >  2:1254953/1‑148 (MQ=255)
   cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                           >  2:1255005/1‑148 (MQ=255)
   cTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAg                                                                                                                                           >  2:1842881/1‑148 (MQ=255)
    ttAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAgg                                                                                                                                          <  1:89762/148‑1 (MQ=255)
    ttAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAgg                                                                                                                                          <  2:1621947/148‑1 (MQ=255)
     tAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGa                                                                                                                                         >  2:3141409/1‑148 (MQ=255)
      aGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGAt                                                                                                                                        <  2:1974950/148‑1 (MQ=255)
       gTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATc                                                                                                                                       <  1:1428696/148‑1 (MQ=255)
         gCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTg                                                                                                                                     <  2:920845/148‑1 (MQ=255)
         gCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTg                                                                                                                                     <  2:2322602/148‑1 (MQ=255)
            aaaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGAtt                                                                                                                                  >  2:199791/1‑148 (MQ=255)
             aaTGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGAttt                                                                                                                                 >  1:3048337/1‑148 (MQ=255)
              aTGGCTGAATAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGAtttt                                                                                                                                <  2:2077618/148‑1 (MQ=255)
               tGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTtc                                                                                                                               >  2:727605/1‑148 (MQ=255)
                ggCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTtca                                                                                                                              >  2:492777/1‑148 (MQ=255)
                 gCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTtcat                                                                                                                             <  2:2653801/148‑1 (MQ=255)
                    gAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAg                                                                                                                          >  1:1289698/1‑148 (MQ=255)
                        aGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCAt                                                                                                                      <  1:1743920/148‑1 (MQ=255)
                          aTTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGc                                                                                                                    >  1:2440487/1‑148 (MQ=255)
                           ttGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGcc                                                                                                                   >  1:1058103/1‑148 (MQ=255)
                            tGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGccc                                                                                                                  >  2:2932642/1‑148 (MQ=255)
                               ttGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTg                                                                                                               >  2:3133570/1‑148 (MQ=255)
                               ttGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTg                                                                                                               >  1:1316015/1‑148 (MQ=255)
                                 gcgcATTGTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGa                                                                                                             <  2:46910/148‑1 (MQ=255)
                                  cgcATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGAt                                                                                                            <  2:1056467/148‑1 (MQ=255)
                                  cgcATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGAt                                                                                                            <  2:2292238/148‑1 (MQ=255)
                                  cgcATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGAt                                                                                                            >  2:342536/1‑148 (MQ=255)
                                    cATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCg                                                                                                          >  1:18197/1‑148 (MQ=255)
                                    cATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCg                                                                                                          >  1:32388/1‑148 (MQ=255)
                                     aTTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGa                                                                                                         >  1:2480481/1‑148 (MQ=255)
                                      ttttGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGaa                                                                                                        >  1:2816908/1‑148 (MQ=255)
                                          gCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTgg                                                                                                    >  1:1742829/1‑148 (MQ=255)
                                            cGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGt                                                                                                  <  2:2316616/148‑1 (MQ=255)
                                              aaCACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCg                                                                                                <  1:824226/148‑1 (MQ=255)
                                              aaCACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCg                                                                                                <  2:2405137/148‑1 (MQ=255)
                                              aaCACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCg                                                                                                >  1:3354648/1‑148 (MQ=255)
                                              aaCACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCg                                                                                                >  1:3504372/1‑148 (MQ=255)
                                               acacCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGt                                                                                               >  1:2292111/1‑148 (MQ=255)
                                                 accacGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTaa                                                                                             >  1:1192510/1‑148 (MQ=255)
                                                   cacGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAAc                                                                                           >  2:872125/1‑148 (MQ=255)
                                                      ggTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAActct                                                                                        <  1:208602/148‑1 (MQ=255)
                                                       gTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTg                                                                                       <  2:1578344/148‑1 (MQ=255)
                                                            gcgATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGc                                                                                  <  1:3471529/148‑1 (MQ=255)
                                                               aTGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCg                                                                               >  2:260306/1‑148 (MQ=255)
                                                                 ggAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCa                                                                             <  2:2001916/148‑1 (MQ=255)
                                                                 ggAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCa                                                                             >  1:1940536/1‑148 (MQ=255)
                                                                 ggAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCa                                                                             >  1:631732/1‑148 (MQ=255)
                                                                  gAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAg                                                                            >  1:1004881/1‑148 (MQ=255)
                                                                   aaGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCCTCAGTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAgg                                                                           >  1:920357/1‑148 (MQ=255)
                                                                    agacgcgCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGCTCGTAAACTCTGGCTGCCCGCAGGc                                                                          <  2:2668024/145‑1 (MQ=255)
                                                                     gTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCg                                                                         <  1:2259837/148‑1 (MQ=255)
                                                                      tCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGt                                                                        >  2:2493613/1‑148 (MQ=255)
                                                                         cgcAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGa                                                                     <  2:2369043/148‑1 (MQ=255)
                                                                              tGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCg                                                                >  2:1973172/1‑148 (MQ=255)
                                                                                ccAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTc                                                              <  1:2267582/148‑1 (MQ=255)
                                                                                    cAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGc                                                          >  1:56513/1‑148 (MQ=255)
                                                                                         gCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTa                                                     <  1:254/148‑1 (MQ=255)
                                                                                          cccAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTAt                                                    >  2:1544211/1‑148 (MQ=255)
                                                                                          cccAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTAt                                                    >  1:407896/1‑148 (MQ=255)
                                                                                                          aGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGaa                                    >  1:2860440/1‑148 (MQ=255)
                                                                                                            cTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAaca                                  >  2:239923/1‑148 (MQ=255)
                                                                                                                cggcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcaccac                              <  1:303233/148‑1 (MQ=255)
                                                                                                                cggcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAAcaccac                              <  1:2749968/148‑1 (MQ=255)
                                                                                                                  gcCCTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACag                            <  1:3511110/148‑1 (MQ=255)
                                                                                                                   cccTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACaga                           >  1:2634919/1‑148 (MQ=255)
                                                                                                                   cccTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGACGGGGTAGTAAAAGCAGGCGTCCCCCCCGCCTTTAAGGCCGCATGCAGCCAAGGGAACCGCCCCCGAAcaccccccag                          >  1:2025570/1‑142 (MQ=255)
                                                                                                                     cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                         <  1:18693/148‑1 (MQ=255)
                                                                                                                     cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                         <  1:17403/148‑1 (MQ=255)
                                                                                                                     cTTGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagaga                         <  1:1388100/148‑1 (MQ=255)
                                                                                                                      ttGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACagagag                        <  1:2148990/148‑1 (MQ=255)
                                                                                                                       tGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAgg                       <  1:426014/148‑1 (MQ=255)
                                                                                                                       tGCGCGTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAgg                       <  2:2502477/148‑1 (MQ=255)
                                                                                                                       tGCGCGTCAGCATCATTGATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAgg                       <  2:1322434/148‑1 (MQ=255)
                                                                                                                         cgcgTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGa                     >  2:1366660/1‑148 (MQ=255)
                                                                                                                         cgcgTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGa                     <  2:3048337/148‑1 (MQ=255)
                                                                                                                         cgcgTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGa                     >  2:1366555/1‑148 (MQ=255)
                                                                                                                          gcgTCAGCATCATTTATGGTCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGaa                    <  1:3298851/148‑1 (MQ=255)
                                                                                                                          gcgTCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGaa                    <  2:110530/148‑1 (MQ=255)
                                                                                                                            gTCAGCATCATTTATGGGCATAACAAGGACCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAAtt                  <  2:1767473/148‑1 (MQ=255)
                                                                                                                             tCAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTa                 <  1:868818/148‑1 (MQ=255)
                                                                                                                              cAGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTAt                >  2:813950/1‑148 (MQ=255)
                                                                                                                               aGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATg               <  1:1254953/148‑1 (MQ=255)
                                                                                                                               aGCATCATTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATg               <  1:5697/148‑1 (MQ=255)
                                                                                                                                  atcatTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGc            <  1:2699320/148‑1 (MQ=255)
                                                                                                                                   tcatTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGcc           >  1:1499424/1‑148 (MQ=255)
                                                                                                                                    catTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCa          <  1:2884640/148‑1 (MQ=255)
                                                                                                                                    catTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCa          <  1:1367509/148‑1 (MQ=255)
                                                                                                                                     atTTATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAg         <  2:1331317/148‑1 (MQ=255)
                                                                                                                                      tttATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGt        <  2:3236703/148‑1 (MQ=255)
                                                                                                                                      tttATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGt        <  2:3501091/148‑1 (MQ=255)
                                                                                                                                       ttATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGACCAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtt       >  1:3326337/1‑148 (MQ=255)
                                                                                                                                       ttATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGtt       <  2:2666024/148‑1 (MQ=255)
                                                                                                                                        tATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttt      <  1:2003667/148‑1 (MQ=255)
                                                                                                                                        tATGGGCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttt      <  1:161287/148‑1 (MQ=255)
                                                                                                                                            ggCATAACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGttttttt  <  2:32388/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AATCTTAGTGCTAAATGGCTGAAAAGATTGGTTGCGCATTTTGCCGAACACCACGGTGACGCGATGGAAGTCGCGCAGTGCCAGCAATCGCCCAATGAACGTACCCAGCTGGCGGCCCTTGCGCGTCAGCATCATTTATGGGCATCACAAGGATCTGATTTTCATCAGCCATGCCCGTGGATCGAACTGGGTCGTAAACTCTGGCTGCCCGCAGGCGTTGAAGGCGTCTGGCAGCTATGGGAACAGCCGCAGAACACCACAGAGAGGGAATTATGAGCCAGTTTTTTT  >  CP009273/1318083‑1318370

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: