Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,116,964 A→G 54.8% intergenic (‑21/+53) dinI ← / ← pyrC DNA damage‑inducible protein I/dihydro‑orotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,116,9640AG54.8% 106.3 / inf 277intergenic (‑21/+53)dinI/pyrCDNA damage‑inducible protein I/dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (55/70);  new base G (58/93);  total (114/163)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 279 total aligned reads displayed.

GTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCT  >  CP009273/1116817‑1117110
                                                                                                                                                   |                                                                                                                                                  
 tACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACt                                                                                                                                                  <  1:3866650/149‑1 (MQ=255)
  aCCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTg                                                                                                                                                 <  2:1555760/149‑1 (MQ=255)
   ccGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATATCCCCCTGTTGAATATGCTgg                                                                                                                                                >  1:3251156/1‑149 (MQ=255)
   ccGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTgg                                                                                                                                                <  2:1192275/149‑1 (MQ=255)
   ccGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCACATAATAGCCCCCTGTTGAATATGCTgg                                                                                                                                                <  2:3073521/149‑1 (MQ=255)
     gATACGTGGCCTTCATTATCAGGAAACGCATACGGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGtt                                                                                                                                              <  2:3322924/149‑1 (MQ=255)
        aCGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttatt                                                                                                                                           >  1:3398514/1‑149 (MQ=255)
        aCGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGttatt                                                                                                                                           <  2:2814773/149‑1 (MQ=255)
          ggggCCTTCATTATCAGGAAACGCAGACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttattta                                                                                                                                         <  1:2990877/147‑1 (MQ=255)
          gTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttattta                                                                                                                                         <  1:4641677/149‑1 (MQ=255)
          gTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttattta                                                                                                                                         >  2:2306765/1‑149 (MQ=255)
             gCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAc                                                                                                                                      >  2:1492682/1‑149 (MQ=255)
               cTTCATTATCAGGAAACGCATACTGACTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAg                                                                                                                                    <  1:936684/149‑1 (MQ=255)
                  cATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTa                                                                                                                                 >  1:1646814/1‑149 (MQ=255)
                    ttATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTaaa                                                                                                                               >  2:3948120/1‑149 (MQ=255)
                     tATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTaaaa                                                                                                                              <  1:863303/149‑1 (MQ=255)
                         aGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAAc                                                                                                                          <  1:3331792/149‑1 (MQ=255)
                         aGGAAACGCAGACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAAc                                                                                                                          <  2:2122414/149‑1 (MQ=255)
                          ggAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAAcc                                                                                                                         <  1:3003082/149‑1 (MQ=255)
                           gAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCt                                                                                                                        <  2:1976594/149‑1 (MQ=255)
                             aaCGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTaa                                                                                                                      <  2:693926/149‑1 (MQ=255)
                              aCGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTGTTCGCTATGGTGACTTCAATTCGCAAAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAAt                                                                                                                     >  1:113568/1‑149 (MQ=255)
                               cGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATg                                                                                                                    <  2:4964533/149‑1 (MQ=255)
                                gCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGa                                                                                                                   <  1:4490194/149‑1 (MQ=255)
                                gCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGa                                                                                                                   <  1:4775909/149‑1 (MQ=255)
                                    aCTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaac                                                                                                               <  1:4135598/149‑1 (MQ=255)
                                    aCTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaac                                                                                                               >  1:2551902/1‑149 (MQ=255)
                                    aCTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGacaac                                                                                                               <  2:1689858/149‑1 (MQ=255)
                                      tGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAg                                                                                                             >  1:5194731/1‑149 (MQ=255)
                                      tGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAg                                                                                                             <  2:1007735/149‑1 (MQ=255)
                                      tGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATGCTGGTTATTTATACAGGTAAAATAAACTAATGACAACAg                                                                                                             >  1:1798145/1‑149 (MQ=255)
                                       gAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAgg                                                                                                            >  1:2610391/1‑149 (MQ=255)
                                        aaTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                                                                                                           >  1:262073/1‑149 (MQ=255)
                                         aTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGaa                                                                                                          >  2:2456938/1‑149 (MQ=255)
                                           aCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        <  1:4306495/149‑1 (MQ=255)
                                           aCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        <  2:4499807/149‑1 (MQ=255)
                                           aCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        <  1:505576/149‑1 (MQ=255)
                                           aCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        <  2:1890440/149‑1 (MQ=255)
                                              gcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTAc                                                                                                     <  1:2519722/149‑1 (MQ=255)
                                              gcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTAc                                                                                                     <  1:2284172/149‑1 (MQ=255)
                                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAt                                                                                                  <  1:420087/149‑1 (MQ=255)
                                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAt                                                                                                  >  2:269952/1‑149 (MQ=255)
                                                 ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAt                                                                                                  <  1:2086/149‑1 (MQ=255)
                                                   aaaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttt                                                                                                <  1:3484482/149‑1 (MQ=255)
                                                      gTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTAGCATAATATCCCCCTGTTGAATATACTGGTTATTTATAAAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTa                                                                                             >  1:2430357/1‑149 (MQ=255)
                                                        tCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAtt                                                                                           <  2:1344114/149‑1 (MQ=255)
                                                        tCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAtt                                                                                           <  2:1344022/149‑1 (MQ=255)
                                                           ccagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGtt                                                                                        <  1:1097232/149‑1 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGttt                                                                                       <  2:87013/149‑1 (MQ=255)
                                                             agccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      <  2:3419418/149‑1 (MQ=255)
                                                              gccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTaa                                                                                     <  2:1209450/149‑1 (MQ=255)
                                                              gccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTaa                                                                                     >  1:1726534/1‑149 (MQ=255)
                                                                  gggCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGa                                                                                 <  2:4027673/149‑1 (MQ=255)
                                                                    gCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGAcc                                                                               <  1:50389/149‑1 (MQ=255)
                                                                     cGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCa                                                                              <  2:4311065/149‑1 (MQ=255)
                                                                     cGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCa                                                                              <  1:2160059/149‑1 (MQ=255)
                                                                           tAGCCCCAGCCGGCAATGGAGAAGTTTTCGCTATGGTGCCTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTAc                                                                        >  2:2201371/1‑149 (MQ=255)
                                                                               cccAGCTGGCAATGGAGTAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGtt                                                                    <  1:1892953/149‑1 (MQ=255)
                                                                                  aGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcc                                                                 >  2:5154640/1‑149 (MQ=255)
                                                                                   gCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTccc                                                                <  1:5005932/149‑1 (MQ=255)
                                                                                    cTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGACTATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                               <  1:2743418/149‑1 (MQ=255)
                                                                                          aTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAg                                                         <  1:577969/149‑1 (MQ=255)
                                                                                             gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGaa                                                      >  2:4740667/1‑149 (MQ=255)
                                                                                             gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGaa                                                      <  2:1795168/149‑1 (MQ=255)
                                                                                             gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGag                                                      <  1:3967096/149‑2 (MQ=255)
                                                                                             gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGaa                                                      <  1:4094380/149‑1 (MQ=255)
                                                                                              agaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAAt                                                     <  1:1077028/149‑1 (MQ=255)
                                                                                              agaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAAt                                                     <  2:3565405/149‑1 (MQ=255)
                                                                                                  gTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCa                                                 >  2:3649548/1‑149 (MQ=255)
                                                                                                   ttttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAc                                                >  1:1491065/1‑149 (MQ=255)
                                                                                                   ttttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAc                                                <  2:1552224/149‑1 (MQ=255)
                                                                                                        gCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCg                                           >  1:2938613/1‑149 (MQ=255)
                                                                                                         cTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgt                                          <  2:1938569/149‑1 (MQ=255)
                                                                                                         cTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgt                                          >  2:4829611/1‑149 (MQ=255)
                                                                                                         cTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgt                                          >  2:3246475/1‑149 (MQ=255)
                                                                                                           aTGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTGATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgtgt                                        >  1:1455724/1‑149 (MQ=255)
                                                                                                           aTGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgtgt                                        >  1:3364223/1‑149 (MQ=255)
                                                                                                               tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatc                                    <  2:2551902/149‑1 (MQ=255)
                                                                                                               tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatc                                    <  2:2685483/149‑1 (MQ=255)
                                                                                                                   ttCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagtc                                >  1:2033463/1‑149 (MQ=255)
                                                                                                                   ttCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagtc                                >  2:3006375/1‑149 (MQ=255)
                                                                                                                    tCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGACTGGCACCAGCGTGTCAgcagtca                               >  2:4561449/1‑149 (MQ=255)
                                                                                                                      aaTTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATcagtcagt                             <  1:817051/149‑1 (MQ=255)
                                                                                                                       aTTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTg                            >  2:3158568/1‑149 (MQ=255)
                                                                                                                          cGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGa                         <  2:722701/149‑1 (MQ=255)
                                                                                                                             ataataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGc                      >  2:4688782/1‑149 (MQ=255)
                                                                                                                             ataataGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGc                      <  2:1894911/149‑1 (MQ=255)
                                                                                                                             ataataGCCCCCTGTTGAATAGACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGc                      <  2:4898543/149‑1 (MQ=255)
                                                                                                                               aataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCtt                    >  1:4826190/1‑149 (MQ=255)
                                                                                                                                ataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCttt                   >  2:789186/1‑149 (MQ=255)
                                                                                                                                 taGCCCCCTGTTGAATAGACTGGTTATTTATACAGGTAAAATCACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc                  <  1:4213960/149‑1 (MQ=255)
                                                                                                                                  aGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCa                 >  2:1911098/1‑149 (MQ=255)
                                                                                                                                   gCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg                >  1:462243/1‑149 (MQ=255)
                                                                                                                                   gCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg                >  1:4533436/1‑149 (MQ=255)
                                                                                                                                    cccccTGGTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGc               <  2:113180/149‑1 (MQ=255)
                                                                                                                                        cTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAAcc           >  2:669919/1‑149 (MQ=255)
                                                                                                                                          gTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTg         <  1:3645543/149‑1 (MQ=255)
                                                                                                                                          gTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTg         <  1:3645525/149‑1 (MQ=255)
                                                                                                                                            tGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGtt       <  2:3226686/149‑1 (MQ=255)
                                                                                                                                              aaTATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGc     >  2:4409277/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCT  >  CP009273/1116817‑1117110

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: