Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,521,860 T→A 12.6% H559L (CAC→CTC)  ligA ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,521,8600TA12.6% 922.3 / 103.7 309H559L (CAC→CTC) ligADNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base T (140/129);  new base A (21/18);  total (162/147)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 311 total aligned reads displayed.

GCTGCCCGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGA  >  CP009273/2521712‑2522009
                                                                                                                                                    |                                                                                                                                                     
gCTGCCCGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTg                                                                                                                                                      >  1:1379985/1‑150 (MQ=255)
  tGCCCGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTggg                                                                                                                                                    <  1:2181864/150‑1 (MQ=255)
  tGCCCGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTggg                                                                                                                                                    <  1:3316505/150‑1 (MQ=255)
    cccGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGAt                                                                                                                                                  >  2:2843895/1‑150 (MQ=255)
    cccGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGAt                                                                                                                                                  <  2:1303384/150‑1 (MQ=255)
    cccGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGAt                                                                                                                                                  >  1:2477887/1‑150 (MQ=255)
     ccGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATg                                                                                                                                                 <  2:1731164/150‑1 (MQ=255)
      cGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGc                                                                                                                                                >  1:737004/1‑150 (MQ=255)
      cGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGc                                                                                                                                                <  1:2071491/150‑1 (MQ=255)
         aaGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCTCTGATGACATTGCGGTTGCTTTCTTCGTCAATGAAGTTGTGAACGAGGGATGAAAc                                                                                                                                             >  1:970496/1‑150 (MQ=255)
            caccacGGTTTTACCAGCAAATGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGAc                                                                                                                                          >  1:413650/1‑150 (MQ=255)
            caccacGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGAc                                                                                                                                          >  2:1082528/1‑150 (MQ=255)
                  ggTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGcc                                                                                                                                    >  1:2091980/1‑150 (MQ=255)
                  ggTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGcc                                                                                                                                    <  1:2306879/150‑1 (MQ=255)
                   gTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCa                                                                                                                                   <  2:4373762/150‑1 (MQ=255)
                   gTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCa                                                                                                                                   >  1:1221792/1‑150 (MQ=255)
                        aCCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATc                                                                                                                              >  1:3749001/1‑150 (MQ=255)
                        aCCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATc                                                                                                                              >  2:4596517/1‑150 (MQ=255)
                        aCCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATc                                                                                                                              >  1:910091/1‑150 (MQ=255)
                         ccAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCa                                                                                                                             <  1:1734925/150‑1 (MQ=255)
                           aGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAgg                                                                                                                           <  1:4469232/150‑1 (MQ=255)
                            gCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGc                                                                                                                          >  2:1111876/1‑150 (MQ=255)
                             cAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCa                                                                                                                         >  2:830722/1‑150 (MQ=255)
                              aaaCGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCAc                                                                                                                        >  1:2654754/1‑150 (MQ=255)
                              aaaCGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCAc                                                                                                                        <  2:3686211/150‑1 (MQ=255)
                                 cGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCtt                                                                                                                     >  2:2021916/1‑150 (MQ=255)
                                  gggCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCttt                                                                                                                    <  2:4750/150‑1 (MQ=255)
                                   ggCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGTAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCtttt                                                                                                                   >  2:4558017/1‑150 (MQ=255)
                                      tGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCa                                                                                                                >  2:4636638/1‑150 (MQ=255)
                                          aaTCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGctc                                                                                                            >  1:2732572/1‑150 (MQ=255)
                                          aaTCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGctc                                                                                                            <  1:2362754/150‑1 (MQ=255)
                                          aaTCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGctc                                                                                                            <  1:2362609/150‑1 (MQ=255)
                                              tctTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCa                                                                                                        >  2:2435079/1‑150 (MQ=255)
                                              tctTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCa                                                                                                        >  2:1561932/1‑150 (MQ=255)
                                               ctTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCaa                                                                                                       <  1:153760/150‑1 (MQ=255)
                                                ttCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAAt                                                                                                      >  1:2713185/1‑150 (MQ=255)
                                                 tCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGTACCTTTTGCAGCTCTTCAATc                                                                                                     >  2:1655016/1‑150 (MQ=255)
                                                 tCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATc                                                                                                     <  1:462523/150‑1 (MQ=255)
                                                 tCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATc                                                                                                     <  2:2443327/150‑1 (MQ=255)
                                                   cgcgTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGa                                                                                                   <  2:1760531/150‑1 (MQ=255)
                                                     cgTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAg                                                                                                 <  2:2765507/150‑1 (MQ=255)
                                                     cgTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAg                                                                                                 >  2:541167/1‑150 (MQ=255)
                                                      gTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAgg                                                                                                <  2:3140519/150‑2 (MQ=255)
                                                            aaCGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttc                                                                                          >  2:913535/1‑150 (MQ=255)
                                                             aCGCTCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACACTGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttcc                                                                                         <  2:232953/150‑1 (MQ=255)
                                                               gATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccac                                                                                       <  1:138884/150‑2 (MQ=255)
                                                               gATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTACAATCGAAGCGgcttccag                                                                                       >  1:1009037/1‑150 (MQ=255)
                                                                 tCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcg                                                                                     <  2:1378053/150‑1 (MQ=255)
                                                                   ggCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcgct                                                                                   >  2:1146144/1‑150 (MQ=255)
                                                                    gcgcAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcgctt                                                                                  <  2:2020361/150‑1 (MQ=255)
                                                                       cAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcgcttcca                                                                               <  2:3387655/150‑1 (MQ=255)
                                                                       cAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcgcttcca                                                                               >  2:3620321/1‑150 (MQ=255)
                                                                       cAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACAGTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCGTTTGCAGCTCTTCAATCGAAGCGgattccagcgcttcca                                                                               >  2:302281/1‑150 (MQ=255)
                                                                        aGGCCAGTGAACTCCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcgcttccag                                                                              <  1:3460907/150‑1 (MQ=255)
                                                                        aGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGgcttccagcgcttccag                                                                              <  2:1955400/150‑1 (MQ=255)
                                                                            cAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTg                                                                          >  2:450250/1‑150 (MQ=255)
                                                                              gTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGcc                                                                        <  1:3918865/150‑1 (MQ=255)
                                                                                gAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGa                                                                      <  1:3865692/150‑1 (MQ=255)
                                                                                   cacCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAat                                                                   >  1:2641679/1‑150 (MQ=255)
                                                                                    acCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAata                                                                  >  1:3623984/1‑150 (MQ=255)
                                                                                          cGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATgccgc                                                            >  1:775545/1‑150 (MQ=255)
                                                                                           gCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATgccgcc                                                           >  2:3897481/1‑150 (MQ=255)
                                                                                              aaCAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGACCGGGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGa                                                        >  1:3180035/1‑150 (MQ=255)
                                                                                                cAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGAcc                                                      >  2:1353981/1‑150 (MQ=255)
                                                                                                 aGCTCGCTTATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAAAGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGTCct                                                     <  1:3346696/150‑1 (MQ=255)
                                                                                                 aGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACct                                                     <  1:3254930/150‑1 (MQ=255)
                                                                                                  gCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACctg                                                    <  2:4357302/150‑1 (MQ=255)
                                                                                                      gCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGCACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCtgctgc                                                <  2:3217523/150‑1 (MQ=255)
                                                                                                       cTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCg                                               <  1:808185/150‑1 (MQ=255)
                                                                                                       cTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCg                                               >  2:891828/1‑150 (MQ=255)
                                                                                                        tgatgaCATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCgg                                              >  2:3558867/1‑150 (MQ=255)
                                                                                                        tgatgaCATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCgg                                              >  1:1095599/1‑150 (MQ=255)
                                                                                                         gatgaCATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGt                                             <  2:4430508/150‑1 (MQ=255)
                                                                                                         gatgaCATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGt                                             <  1:1422897/150‑1 (MQ=255)
                                                                                                           tgaCATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTgg                                           >  2:3445840/1‑150 (MQ=255)
                                                                                                           tgaCATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTgg                                           >  1:1838720/1‑150 (MQ=255)
                                                                                                              cATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCt                                        <  1:3519112/150‑1 (MQ=255)
                                                                                                                ttGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCg                                      <  2:2090435/150‑1 (MQ=255)
                                                                                                                ttGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCg                                      >  2:3961255/1‑150 (MQ=255)
                                                                                                                  gCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGcc                                    <  1:2829640/150‑1 (MQ=255)
                                                                                                                       tGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACtt                               <  1:2850009/150‑1 (MQ=255)
                                                                                                                        gCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTc                              <  1:1954532/150‑1 (MQ=255)
                                                                                                                          tttCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCAc                            <  1:215341/150‑1 (MQ=255)
                                                                                                                          tttCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCAc                            <  2:2645076/150‑1 (MQ=255)
                                                                                                                              ttcGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGAt                        <  1:543769/150‑1 (MQ=255)
                                                                                                                              ttcGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGAt                        >  1:4329418/1‑150 (MQ=255)
                                                                                                                              ttcGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGAt                        <  1:1664076/150‑1 (MQ=255)
                                                                                                                                cGGCAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGc                      >  2:141462/1‑150 (MQ=255)
                                                                                                                                  gCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCa                    <  1:1455364/150‑1 (MQ=255)
                                                                                                                                   cAAAGAAGTTGTGAACGAGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCaa                   >  1:4000373/1‑150 (MQ=255)
                                                                                                                                     aagaagTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGt                 >  1:3071493/1‑150 (MQ=255)
                                                                                                                                     aagaagTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGt                 >  2:219394/1‑150 (MQ=255)
                                                                                                                                     aagaagTTGTGAACGAGGGATGCAACGACAATGCCACCATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGt                 <  1:289569/150‑1 (MQ=255)
                                                                                                                                      agaagTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTg                >  1:4385122/1‑150 (MQ=255)
                                                                                                                                          gTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATa            <  2:2887183/150‑1 (MQ=255)
                                                                                                                                           ttGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATag           <  2:1233300/150‑1 (MQ=255)
                                                                                                                                            tgtgAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATaga          >  1:4537021/1‑150 (MQ=255)
                                                                                                                                                   gTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCg   >  2:3281679/1‑150 (MQ=255)
                                                                                                                                                   gTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCg   <  2:3141181/150‑1 (MQ=255)
                                                                                                                                                    tGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGa  <  2:4130278/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTGCCCGTAAGCACCACGGTTTTACCAGCAAACGGGCTGTCAATCTCTTCCGCGTTGATAACGATCGGCGCAGGCCAGTGAACACCTTCCGCCAACAGCTCGCTGATGACATTGCGGTTGCTTTCTTCGGCAAAGAAGTTGTGAACGTGGGATGCAACGACAATGCCAACATCAGGCACCTTTTGCAGCTCTTCAATCGAAGCGGCTTCCAGCGCTTCCAGCGTGCCGAAATATGCCGCCAGACCTGCTGCGGTGGCCTCGCCGACTTCACGGATGCCAAGTGCATAGAGGAAGCGA  >  CP009273/2521712‑2522009

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: