Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 408,940 G→T 22.1% A757E (GCG→GAG)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273408,9400GT22.1% 709.1 / 164.5 278A757E (GCG→GAG) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (101/115);  new base T (34/28);  total (135/143)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

GCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTTTCTTCGCAGTGCGG  >  CP009273/408791‑409088
                                                                                                                                                     |                                                                                                                                                    
  gTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGcgcg                                                                                                                                                    >  2:1132400/1‑150 (MQ=255)
   tttCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTTTTCATACAATAGCGCCGCAAGGAACGCCTGCTGATCGTCAAATAAGCTGGCCAGTATCTCgg                                                                                                                                                   >  1:3837232/1‑150 (MQ=255)
     tCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGgtg                                                                                                                                                 >  2:4264459/1‑150 (MQ=255)
       tgctgTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTc                                                                                                                                               <  1:1443559/150‑1 (MQ=255)
       tgctgTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTc                                                                                                                                               <  1:2248360/150‑1 (MQ=255)
        gctgTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCa                                                                                                                                              <  1:107724/150‑1 (MQ=255)
        gctgTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCa                                                                                                                                              <  2:3987467/150‑1 (MQ=255)
         ctgTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTTTTCATCCAATAGCGCCGCAAGGAACGCCAGCTGATAGTCAAAGACGCAGGCCTGTAGCTCGGTGTCaa                                                                                                                                             >  1:4014454/1‑150 (MQ=255)
         ctgTCTGAGTGACCAGAGTTTGCGCCAGACGGCGCTGGTTTTACAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTTTTCATCCAATAGCGCCGCAAGGAAACCCAGCTGATCGTCAAAGACGCTGGCCTGTATCGCGGTGTCaa                                                                                                                                             >  1:663712/1‑150 (MQ=255)
          tgTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCaaa                                                                                                                                            >  1:3605707/1‑150 (MQ=255)
             cTGAGTGACGAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGGTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCGGATCGTCAAAGACGCGGGCCTGTAGCTCGGTGTCAAACTg                                                                                                                                         >  2:4609669/1‑150 (MQ=255)
                  tGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTACAGATTCTGCTTTAGCTGTTCCAGCTGCGTTAGTGTTTGTTAATCCATTATCGCCCCAAGGAACGCCTGCTGATCGTCAAAGACGCGGGCCAGTAGCTCGGTGTCAAAAAGCGCCt                                                                                                                                    >  1:2724869/1‑150 (MQ=255)
                     ccAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTATTGTTTGTTCATACATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTggg                                                                                                                                 >  1:2824726/1‑150 (MQ=255)
                     ccAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTATGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTAATACATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAATACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTggg                                                                                                                                 >  2:1868295/1‑150 (MQ=255)
                        gagTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCtt                                                                                                                              <  2:4165145/150‑1 (MQ=255)
                            ttGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGGAGCGCGGTGTCAAACTGCGCCGGGGCGGTTTg                                                                                                                          >  1:2038044/1‑150 (MQ=255)
                             tGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGc                                                                                                                         >  1:395911/1‑150 (MQ=255)
                                   tGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGGTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACtt                                                                                                                   >  2:937201/1‑150 (MQ=255)
                                    gACGGCGCTTGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACttt                                                                                                                  <  2:666548/150‑1 (MQ=255)
                                    gACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACttt                                                                                                                  >  1:4604987/1‑150 (MQ=255)
                                    gACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACttt                                                                                                                  <  2:1276555/150‑1 (MQ=255)
                                    gACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACttt                                                                                                                  >  1:1791435/1‑150 (MQ=255)
                                     aCGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTg                                                                                                                 <  2:2189690/150‑1 (MQ=255)
                                      cGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTgc                                                                                                                <  2:3196828/150‑1 (MQ=255)
                                       ggCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTgcg                                                                                                               <  2:1057272/150‑1 (MQ=255)
                                       ggCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTgcg                                                                                                               <  1:2774506/150‑1 (MQ=255)
                                       ggCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTgcg                                                                                                               >  1:704367/1‑150 (MQ=255)
                                         cgcTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGcgcc                                                                                                             <  2:3273525/150‑1 (MQ=255)
                                            tGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCgccgcc                                                                                                          <  1:4550542/150‑1 (MQ=255)
                                              gTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATGAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAg                                                                                                        <  2:244484/150‑1 (MQ=255)
                                               ttttCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGa                                                                                                       <  2:600583/150‑1 (MQ=255)
                                                  tCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGaccc                                                                                                    >  1:3391742/1‑147 (MQ=255)
                                                  tCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACa                                                                                                    <  2:1484849/150‑1 (MQ=255)
                                                  tCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACa                                                                                                    <  2:3354954/150‑1 (MQ=255)
                                                   ccAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACAt                                                                                                   <  2:2137071/150‑1 (MQ=255)
                                                       aTTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTg                                                                                               <  2:3803624/150‑1 (MQ=255)
                                                       aTTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTg                                                                                               <  2:607644/150‑1 (MQ=255)
                                                         tCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCACCGCCAGAACATCCTGtt                                                                                             <  2:3877470/150‑1 (MQ=255)
                                                         tCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGtt                                                                                             <  1:176673/150‑1 (MQ=255)
                                                         tCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGtt                                                                                             <  1:4444008/150‑1 (MQ=255)
                                                         tCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGtt                                                                                             >  1:4509442/1‑150 (MQ=255)
                                                             cTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTg                                                                                         >  2:4042028/1‑150 (MQ=255)
                                                             cTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTg                                                                                         <  1:2696031/150‑1 (MQ=255)
                                                             cTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTg                                                                                         >  2:4461244/1‑150 (MQ=255)
                                                             cTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTg                                                                                         <  1:2344496/150‑1 (MQ=255)
                                                             cTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTg                                                                                         <  1:786443/150‑1 (MQ=255)
                                                             cTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTg                                                                                         >  2:4128671/1‑150 (MQ=255)
                                                                  gCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACg                                                                                    <  2:1876922/150‑1 (MQ=255)
                                                                    tGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTc                                                                                  >  2:3651258/1‑150 (MQ=255)
                                                                      ttCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTctg                                                                                >  1:122312/1‑150 (MQ=255)
                                                                         cAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTctgctg                                                                             <  1:2071139/150‑1 (MQ=255)
                                                                           gCTGCGTTAGTGTTTGTTCATCCATTAACGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTctgctgct                                                                           <  1:2535741/150‑1 (MQ=255)
                                                                             tGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTgctgg                                                                         >  1:3472008/1‑150 (MQ=255)
                                                                              gCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTgctggc                                                                        >  1:4671764/1‑150 (MQ=255)
                                                                                 ttAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCtgt                                                                     >  2:2464163/1‑150 (MQ=255)
                                                                                  tAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCtgtg                                                                    <  1:3912249/150‑1 (MQ=255)
                                                                                  tAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCtgtg                                                                    >  1:2150409/1‑150 (MQ=255)
                                                                                  tAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCtgtg                                                                    >  1:404220/1‑150 (MQ=255)
                                                                                    gtgtTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTa                                                                  >  1:2493322/1‑150 (MQ=255)
                                                                                    gtgtTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTa                                                                  <  2:2672621/150‑1 (MQ=255)
                                                                                    gtgtTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTa                                                                  >  1:1015327/1‑150 (MQ=255)
                                                                                      gtttgttCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAAt                                                                >  1:3006305/1‑150 (MQ=255)
                                                                                      gtttgttCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAAt                                                                <  2:646188/150‑1 (MQ=255)
                                                                                      gtttgttCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAAt                                                                <  1:3088342/150‑1 (MQ=255)
                                                                                           ttCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCgag                                                           >  2:2768858/1‑150 (MQ=255)
                                                                                            tCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCgaga                                                          <  1:3892065/150‑1 (MQ=255)
                                                                                             cATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGAc                                                         >  2:4856344/1‑150 (MQ=255)
                                                                                              aTCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACa                                                        >  1:1548086/1‑150 (MQ=255)
                                                                                                ccATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACAtt                                                      <  2:2999144/150‑1 (MQ=255)
                                                                                                 cATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTg                                                     <  2:355586/150‑1 (MQ=255)
                                                                                                 cATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTg                                                     >  1:4574079/1‑150 (MQ=255)
                                                                                                  aTTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGt                                                    <  1:848294/150‑1 (MQ=255)
                                                                                                          cgcAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTAc                                            >  1:4146741/1‑150 (MQ=255)
                                                                                                             aaGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTg                                         >  1:2575008/1‑150 (MQ=255)
                                                                                                              aGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTgc                                        >  2:1446121/1‑150 (MQ=255)
                                                                                                                gAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTgccg                                      >  1:2632245/1‑150 (MQ=255)
                                                                                                                gAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTgccg                                      <  2:2583630/150‑1 (MQ=255)
                                                                                                                   cGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCa                                   <  2:1248651/150‑1 (MQ=255)
                                                                                                                      ctgctgATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGtt                                >  2:1617060/1‑150 (MQ=255)
                                                                                                                       tgctgATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGttt                               >  1:3707914/1‑150 (MQ=255)
                                                                                                                        gctgATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGtttt                              >  2:4221049/1‑150 (MQ=255)
                                                                                                                        gctgATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGtttt                              <  1:1605329/150‑1 (MQ=255)
                                                                                                                         ctgATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTc                             <  2:697071/150‑1 (MQ=255)
                                                                                                                         ctgATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTc                             >  2:410658/1‑150 (MQ=255)
                                                                                                                            aTCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCaa                          <  2:1061821/150‑1 (MQ=255)
                                                                                                                            aTCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCaa                          <  1:1471875/150‑1 (MQ=255)
                                                                                                                             tCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAAt                         <  2:4363943/150‑1 (MQ=255)
                                                                                                                              cGTCAAAGACGCTGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATa                        >  1:2403379/1‑150 (MQ=255)
                                                                                                                                   aaGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACcaca                   >  2:2224423/1‑150 (MQ=255)
                                                                                                                                     gACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGATTTTTGCAGAATTTGCACCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTATTCCAATACCACAGt                 >  1:435863/1‑150 (MQ=255)
                                                                                                                                       cGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGttt               <  2:3063094/150‑1 (MQ=255)
                                                                                                                                        gCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTttc              <  2:922146/150‑1 (MQ=255)
                                                                                                                                          tGGCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTttctt            >  2:328301/1‑150 (MQ=255)
                                                                                                                                          tGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTttctt            <  2:3516295/150‑1 (MQ=255)
                                                                                                                                           ggCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTttcttc           <  1:645384/150‑1 (MQ=255)
                                                                                                                                           ggCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTttcttc           <  1:1868321/150‑1 (MQ=255)
                                                                                                                                            gCCTGTAGCTCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACCTTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTTTCTTCg          >  2:2934408/1‑150 (MQ=255)
                                                                                                                                             ccTGGAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACCGTTTCTTCGc         <  2:2038044/150‑1 (MQ=255)
                                                                                                                                                 tAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTTTCTTCGCAGTg     >  2:2847367/1‑150 (MQ=255)
                                                                                                                                                  aGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTTTCTTCGCAGTGc    >  2:1508396/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGTTTCTGCTGTCTGAGTGACCAGAGTTTGCGCCTGACGGCGCTGGTTTTCCAGATTCTGCTTGAGCTGTTCCAGCTGCGTTAGTGTTTGTTCATCCATTAGCGCCGCAAGGAACGCCTGCTGATCGTCAAAGACGCTGGCCTGTAGCGCGGTGTCAAACTGCGCCTGGGCTTTTTGCAGACTTTGCGCCGCCAGAACATCCTGTTGCTGTAACGTCTGCTGCTGGCTGTGTAATGCGAGACATTGTTCATGTACCTGCCGCCAGTTTTCCAATACCACAGTTTCTTCGCAGTGCGG  >  CP009273/408791‑409088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: