Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,116,964 A→G 61.4% 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,9640AG61.4% 182.5 / 262.5 221intergenic (‑21/+53)dinI/pyrCDNA damage‑inducible protein I/dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (37/48);  new base G (63/71);  total (102/119)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 222 total aligned reads displayed.

GTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCT  >  CP009273/1116817‑1117110
                                                                                                                                                   |                                                                                                                                                  
gTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCt                                                                                                                                                  >  2:441116/1‑150 (MQ=255)
 tACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTg                                                                                                                                                 >  1:1488368/1‑150 (MQ=255)
 tACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTg                                                                                                                                                 <  2:458474/150‑1 (MQ=255)
 tACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTAGCTATGGTGACTTCAATTAGCATAATATCCCCCTGTTGAATATGCTg                                                                                                                                                 >  1:1629258/1‑150 (MQ=255)
   ccGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGt                                                                                                                                               >  1:1579167/1‑150 (MQ=255)
     gATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGtta                                                                                                                                             >  1:2308764/1‑150 (MQ=255)
     gATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTAGCATAATAGCCCCCTGTTGAATATACTGGtta                                                                                                                                             >  1:4663230/1‑150 (MQ=255)
             gCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACa                                                                                                                                     <  2:306930/150‑1 (MQ=255)
             gCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACa                                                                                                                                     >  1:2499957/1‑150 (MQ=255)
              ccTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAg                                                                                                                                    <  2:2596724/150‑1 (MQ=255)
               cTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAAAgg                                                                                                                                   >  1:312666/1‑150 (MQ=255)
                     tATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAAt                                                                                                                             <  2:280785/150‑1 (MQ=255)
                        cAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAAc                                                                                                                          <  2:792590/150‑1 (MQ=255)
                         aGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAAcc                                                                                                                         <  2:2272058/150‑1 (MQ=255)
                         aGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTGTTAGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATTCTGGTTATTTATACAGGTAAAATAAcc                                                                                                                         >  1:4094732/1‑150 (MQ=255)
                          ggAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCt                                                                                                                        <  2:2657105/150‑1 (MQ=255)
                            aaaCGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTaa                                                                                                                      >  2:1710777/1‑150 (MQ=255)
                             aaCGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAAt                                                                                                                     >  2:3993864/1‑150 (MQ=255)
                             aaCGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAAt                                                                                                                     >  1:4853637/1‑150 (MQ=255)
                                gCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGac                                                                                                                  <  2:3760006/150‑1 (MQ=255)
                                 cATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGaca                                                                                                                 >  1:693213/1‑150 (MQ=255)
                                 cATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAAAAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTAAGGTGACATCAATTCGCATAATAGCCCCCTGGTGAATATGCTCGTTAGTTATACAGGTAAAATAACCAAATGaca                                                                                                                 >  1:692862/1‑150 (MQ=255)
                                  aTACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGacaa                                                                                                                <  1:3250476/150‑1 (MQ=255)
                                  aTACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGacaa                                                                                                                >  1:2394241/1‑150 (MQ=255)
                                  aTACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATACTGGTTATTTATAAAGGTAAAATAACCTAATTacaa                                                                                                                >  1:4085812/1‑150 (MQ=255)
                                   tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaac                                                                                                               <  2:173806/150‑1 (MQ=255)
                                   tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaac                                                                                                               <  2:1291448/150‑1 (MQ=255)
                                     cTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAAAGGTAAAATAACCTAATGACAACAg                                                                                                             >  1:47902/1‑150 (MQ=255)
                                       gAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                                                                                                           <  1:3490308/150‑1 (MQ=255)
                                       gAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                                                                                                           <  1:4198644/150‑1 (MQ=255)
                                       gAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACGGGTTATTTATACAGGTAAAATAACCTAATGACACCAGGa                                                                                                           <  2:3108324/150‑1 (MQ=255)
                                        aaTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGaa                                                                                                          >  2:3522318/1‑150 (MQ=255)
                                        aaTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTACAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGaa                                                                                                          >  2:2737200/1‑150 (MQ=255)
                                         aTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAg                                                                                                         <  2:3960807/150‑1 (MQ=255)
                                          tACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        >  2:3903610/1‑150 (MQ=255)
                                          tACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                        >  1:1741606/1‑150 (MQ=255)
                                           aCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCt                                                                                                       <  1:1083847/150‑1 (MQ=255)
                                            cggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTa                                                                                                      <  2:789825/150‑1 (MQ=255)
                                                     aGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTAGCTATGGTGACTTAAATTCGCCTAATAGCCCCCTGTTGAATATACTGGTTATTTATAAACGTAAAATAACCGAATTACAACACGAAGCTCCGATTTTTa                                                                                             >  1:1217384/1‑150 (MQ=255)
                                                      gTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAt                                                                                            <  2:3140418/150‑1 (MQ=255)
                                                        tCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTg                                                                                          >  1:4271012/1‑150 (MQ=255)
                                                          gccagccaGGGCGTCAATAGCCCCAGCTGGCACTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGtt                                                                                        <  1:2653487/150‑1 (MQ=255)
                                                           ccagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGttt                                                                                       <  2:2499720/150‑1 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      <  2:924469/150‑1 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      <  1:2976762/150‑1 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      >  1:4757036/1‑150 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      >  2:732021/1‑150 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      <  2:1430681/150‑1 (MQ=255)
                                                            cagccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                      <  1:3886769/150‑1 (MQ=255)
                                                               ccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACg                                                                                   <  1:377695/150‑1 (MQ=255)
                                                                  gggCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCCTAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGAc                                                                                >  1:1476284/1‑150 (MQ=255)
                                                                    gCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCa                                                                              <  1:1344929/150‑1 (MQ=255)
                                                                       tCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCg                                                                           >  2:814246/1‑150 (MQ=255)
                                                                         aaTAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTa                                                                         <  2:1274540/150‑1 (MQ=255)
                                                                           tAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTAcc                                                                       <  1:3826316/150‑1 (MQ=255)
                                                                            aGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATACAGGTTATTTATACAGGTAAAATAACCAAATGACAACAGGAAGCTACGATTTTTATTTTTTAACGGACCAGCGTACCg                                                                      >  1:4590627/1‑150 (MQ=255)
                                                                               cccAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGttt                                                                   <  1:1988091/150‑1 (MQ=255)
                                                                                  aGCTGGCAATGGAGAAGTTTTCGCTATTGTGACTTCAAATCGCATAAAATCCCCCAGTTGAAAATTCTGGTTATTTATACAGGTAAAATAACCTAATTACAACAAGAAGCTACGATTTTTATTTTTTAAAGGAACAGCGTAACCTTTacc                                                                >  1:3443411/1‑150 (MQ=255)
                                                                                  aGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTccc                                                                >  1:4300754/1‑150 (MQ=255)
                                                                                   gCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                               <  1:3243238/150‑1 (MQ=255)
                                                                                   gCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                               <  1:3488288/150‑1 (MQ=255)
                                                                                      ggCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTAGCATAATATCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCAAATGACAACAGGAAGCTCCGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGc                                                            >  2:2937348/1‑150 (MQ=255)
                                                                                           tGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                       <  1:2144061/150‑1 (MQ=255)
                                                                                              agaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATg                                                    <  2:4047867/150‑1 (MQ=255)
                                                                                               gaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                                   <  2:923938/150‑1 (MQ=255)
                                                                                               gaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                                   >  1:294212/1‑150 (MQ=255)
                                                                                                aaGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGc                                                  >  1:2492425/1‑150 (MQ=255)
                                                                                                   ttttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAcc                                               >  2:2852381/1‑150 (MQ=255)
                                                                                                     ttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAg                                             <  1:2815658/150‑1 (MQ=255)
                                                                                                       cGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCg                                           <  1:4242801/150‑1 (MQ=255)
                                                                                                           aTGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtc                                       >  1:660182/1‑150 (MQ=255)
                                                                                                           aTGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtc                                       <  2:3644888/150‑1 (MQ=255)
                                                                                                             ggTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcat                                     <  2:106281/150‑1 (MQ=255)
                                                                                                               tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGACGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                                   >  2:1857417/1‑150 (MQ=255)
                                                                                                                gACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcag                                  <  2:4276774/150‑1 (MQ=255)
                                                                                                                gACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcag                                  <  2:4034517/150‑1 (MQ=255)
                                                                                                                 aCTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagt                                 >  1:3211371/1‑150 (MQ=255)
                                                                                                                      aaTTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGACTGGCACCAGCGTGTCATCAGTCAGTg                            >  2:4684771/1‑150 (MQ=255)
                                                                                                                          cGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGAt                        <  2:1510586/150‑1 (MQ=255)
                                                                                                                            cATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGc                      >  2:1633020/1‑150 (MQ=255)
                                                                                                                             ataataGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCt                     <  2:568905/150‑1 (MQ=255)
                                                                                                                               aataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCttt                   >  2:1697493/1‑150 (MQ=255)
                                                                                                                               aataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCttt                   <  2:2372433/150‑1 (MQ=255)
                                                                                                                               aataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCttt                   <  2:4212333/150‑1 (MQ=255)
                                                                                                                               aataGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCttt                   >  2:4854206/1‑150 (MQ=255)
                                                                                                                                ataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc                  <  1:3916602/150‑1 (MQ=255)
                                                                                                                                 taGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCa                 <  2:369476/150‑1 (MQ=255)
                                                                                                                                  aGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg                >  1:3877786/1‑150 (MQ=255)
                                                                                                                                    cccccTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGACTGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCa              >  1:2775854/1‑150 (MQ=255)
                                                                                                                                         tGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTg         <  2:3320656/150‑1 (MQ=255)
                                                                                                                                         tGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTg         <  2:2396266/150‑1 (MQ=255)
                                                                                                                                         tGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTg         <  1:4449094/150‑1 (MQ=255)
                                                                                                                                          gTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt        <  2:4757036/150‑1 (MQ=255)
                                                                                                                                          gTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt        >  2:726440/1‑150 (MQ=255)
                                                                                                                                          gTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt        >  1:4282642/1‑150 (MQ=255)
                                                                                                                                             gAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGc     >  1:4133662/1‑150 (MQ=255)
                                                                                                                                              aaTATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGACTGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGct    >  2:3510159/1‑150 (MQ=255)
                                                                                                                                              aaTATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGct    >  2:4845660/1‑150 (MQ=255)
                                                                                                                                              aaTATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGct    >  1:3589581/1‑150 (MQ=255)
                                                                                                                                                tatgCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGctct  <  1:814246/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCT  >  CP009273/1116817‑1117110

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: