Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,116,964 A→G 100% 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,9640AG99.5% 735.4 / ‑4.8 188intergenic (‑21/+53)dinI/pyrCDNA damage‑inducible protein I/dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (83/103);  total (85/103)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 188 total aligned reads displayed.

ACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTG  >  CP009273/1116819‑1117106
                                                                                                                                                 |                                                                                                                                              
aCCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCt                                                                                                                                              <  1:2340455/148‑1 (MQ=255)
aCCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCt                                                                                                                                              <  2:302616/148‑1 (MQ=255)
aCCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCt                                                                                                                                              <  2:302698/148‑1 (MQ=255)
  cGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTgg                                                                                                                                            >  1:1285086/1‑148 (MQ=255)
   gATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGt                                                                                                                                           <  2:1899769/148‑1 (MQ=255)
     tACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGtta                                                                                                                                         <  2:2559279/148‑1 (MQ=255)
     tACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGtta                                                                                                                                         <  1:2730630/148‑1 (MQ=255)
      aCGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttat                                                                                                                                        <  2:874288/148‑1 (MQ=255)
       cGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttatt                                                                                                                                       <  2:2105370/148‑1 (MQ=255)
          ggCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttatttat                                                                                                                                    <  2:723691/148‑1 (MQ=255)
            ccTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAc                                                                                                                                  <  1:236066/148‑1 (MQ=255)
            ccTTCATTATCAGGAAACGCATACTGAATACGGCGCGAAAGTTCGCCAGCCAGGGAGTCAAAACCCCCAGCTGGCAAAGGAGAAGTTTACGGGATGGTGACATCAATTAGAATAATACACCCCACTTGAAAATTGTGGTTTATtataa                                                                                                                                  >  1:1203736/1‑147 (MQ=255)
             cTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACa                                                                                                                                 <  2:443540/148‑1 (MQ=255)
             cTTCATAATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAAATCGCATAATAGCCACCTGTTGAATATGCTGGTTATTTATACa                                                                                                                                 >  2:1390231/1‑148 (MQ=255)
                cATTATCAGGAAACGCATACTGAATACGGCGGGAAAGGTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGt                                                                                                                              <  2:1546555/148‑1 (MQ=255)
                  ttATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTaa                                                                                                                            >  2:3082535/1‑148 (MQ=255)
                   tATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTaaa                                                                                                                           <  1:3119299/148‑1 (MQ=255)
                   tATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTaaa                                                                                                                           <  1:1020562/148‑1 (MQ=255)
                   tATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAAAGGTaaa                                                                                                                           >  2:3156907/1‑148 (MQ=255)
                    aTCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTaaaa                                                                                                                          <  1:1000170/148‑1 (MQ=255)
                      cAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATa                                                                                                                        >  1:225051/1‑148 (MQ=255)
                        ggAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAg                                                                                                                                >  2:657633/1‑138 (MQ=255)
                             cGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAAt                                                                                                                 <  1:463382/148‑1 (MQ=255)
                               cATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGa                                                                                                               <  1:1114517/148‑1 (MQ=255)
                                aTACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTGTTCGCGATGGTGTCTTCAATTCGCATAATAGCCCCCAGTTGTATATGCTGGTTATTTATACAGGTAAAATAACCTAATGac                                                                                                              >  1:1441859/1‑148 (MQ=255)
                                 tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGaca                                                                                                             >  1:750095/1‑148 (MQ=255)
                                       aTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                                                                                                       >  1:2031903/1‑148 (MQ=255)
                                       aTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                                                                                                       <  1:2595314/148‑1 (MQ=255)
                                        tACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATGCTGGTTATTTATACAGGTAACATAACCTAATGACAACAGGaa                                                                                                      >  2:2777425/1‑148 (MQ=255)
                                         aCGGCGGGAAAGTTCGCCAGCCTGGGCGTCAATAGCACCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAg                                                                                                     <  2:1203736/148‑1 (MQ=255)
                                          cggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACACCAGGAAGc                                                                                                    <  2:692289/148‑1 (MQ=255)
                                          cggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                    <  2:692165/148‑1 (MQ=255)
                                          cggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGc                                                                                                    <  2:1020022/148‑1 (MQ=255)
                                            gcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTa                                                                                                  <  1:1870694/148‑1 (MQ=255)
                                                 aaaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAtt                                                                                             <  1:1664536/148‑1 (MQ=255)
                                                 aaaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAtt                                                                                             <  1:1908924/148‑1 (MQ=255)
                                                  aaGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttt                                                                                            <  1:2275618/148‑1 (MQ=255)
                                                   aGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAtttt                                                                                           <  1:2449957/148‑1 (MQ=255)
                                                    gTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                          <  2:523034/148‑1 (MQ=255)
                                                    gTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                          <  2:1473028/148‑1 (MQ=255)
                                                    gTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt                                                                                          <  2:523642/148‑1 (MQ=255)
                                                      tCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAt                                                                                        <  1:1955128/148‑1 (MQ=255)
                                                        gccagccaGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGACGCTACGATTTTTATTg                                                                                      <  1:2332707/148‑1 (MQ=255)
                                                        gccagccaGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTg                                                                                      <  1:672418/148‑1 (MQ=255)
                                                         ccagccagGGCGTCAATAGCCCCATCTGGCAATGGAGAAGTTTTCGCTTTGGTGTCTTCTATTCTCATAATAGCCCCCTGTTGAATTTGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGt                                                                                     <  1:1167988/148‑1 (MQ=255)
                                                           agccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGttt                                                                                   <  1:1274095/148‑1 (MQ=255)
                                                            gccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                  <  1:1787404/148‑1 (MQ=255)
                                                            gccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTa                                                                                  <  1:641313/148‑1 (MQ=255)
                                                             ccagGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTaa                                                                                 <  1:1576148/148‑1 (MQ=255)
                                                                 ggCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATTGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATAAGCTGGTTATTTATAAAGGTAAAATAACCTAATGACAAAAGGAAGCTACGAATTTTAATGTTTAACGGa                                                                             >  2:2595890/1‑148 (MQ=255)
                                                                 ggCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCAAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGa                                                                             >  1:2181972/1‑148 (MQ=255)
                                                                         tAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTa                                                                     <  2:962530/148‑1 (MQ=255)
                                                                           gCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTAcc                                                                   >  2:1575035/1‑148 (MQ=255)
                                                                            ccccAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCg                                                                  >  1:2414086/1‑148 (MQ=255)
                                                                             cccAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGt                                                                 <  2:320569/148‑1 (MQ=255)
                                                                               cAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAAAGGACCAGCGTACCGttt                                                               >  1:2658327/1‑148 (MQ=255)
                                                                               cAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAAAGGTAAAATAACCAAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGttt                                                               >  1:2658177/1‑148 (MQ=255)
                                                                                 gCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcc                                                             <  1:384527/148‑1 (MQ=255)
                                                                                       aaTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCg                                                       <  2:453111/148‑1 (MQ=255)
                                                                                       aaTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCg                                                       <  2:1508752/148‑1 (MQ=255)
                                                                                        aTGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATATCCCCCAGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGa                                                      >  2:1133652/1‑148 (MQ=255)
                                                                                         tGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAg                                                     >  1:536539/1‑148 (MQ=255)
                                                                                          ggAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    >  2:1274823/1‑148 (MQ=255)
                                                                                           gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   <  2:2911568/148‑1 (MQ=255)
                                                                                           gagaAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   <  1:2480625/148‑1 (MQ=255)
                                                                                                gTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGc                                              >  1:3292109/1‑148 (MQ=255)
                                                                                                  tttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAc                                            >  1:245771/1‑148 (MQ=255)
                                                                                                  tttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAc                                            >  1:375934/1‑148 (MQ=255)
                                                                                                   ttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGACTGGCAcc                                           >  2:2265797/1‑148 (MQ=255)
                                                                                                   ttCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCAcc                                           <  1:1141327/148‑1 (MQ=255)
                                                                                                        tATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgt                                      <  2:1230561/148‑1 (MQ=255)
                                                                                                          tGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgtgt                                    <  2:1128587/148‑1 (MQ=255)
                                                                                                          tGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCgtgt                                    <  1:796655/148‑1 (MQ=255)
                                                                                                           ggTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACTATTTTTATTGTTTAACGGACCAGAGTACAGTTTCACCGGCGAGGAATGGCACCAGCGTGtc                                   >  1:156253/1‑148 (MQ=255)
                                                                                                           ggTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtc                                   >  1:70063/1‑148 (MQ=255)
                                                                                                           ggTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtc                                   >  1:70069/1‑148 (MQ=255)
                                                                                                             tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcat                                 >  2:776491/1‑148 (MQ=255)
                                                                                                             tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcat                                 <  1:1415208/148‑1 (MQ=255)
                                                                                                             tGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcat                                 <  1:2935336/148‑1 (MQ=255)
                                                                                                               aCTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                               <  2:23295/148‑1 (MQ=255)
                                                                                                                cTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcag                              >  2:3209997/1‑148 (MQ=255)
                                                                                                                 ttCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagt                             >  2:2677243/1‑148 (MQ=255)
                                                                                                                   cAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagtca                           >  2:2420554/1‑148 (MQ=255)
                                                                                                                   cAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagtca                           >  2:2420438/1‑148 (MQ=255)
                                                                                                                   cAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCAtcagtca                           >  2:1656248/1‑148 (MQ=255)
                                                                                                                     aTTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATcagtcagt                         >  2:2242287/1‑148 (MQ=255)
                                                                                                                       tCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       >  1:418891/1‑148 (MQ=255)
                                                                                                                         gCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGa                     <  1:1530858/148‑1 (MQ=255)
                                                                                                                          cATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGAt                    <  1:3060358/148‑1 (MQ=255)
                                                                                                                           ataataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATg                   <  1:2560509/148‑1 (MQ=255)
                                                                                                                             aataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCt                 <  1:1820029/148‑1 (MQ=255)
                                                                                                                              ataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCtt                >  1:951169/1‑148 (MQ=255)
                                                                                                                              ataGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCtt                >  1:1575749/1‑148 (MQ=255)
                                                                                                                               taGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCttt               >  1:706345/1‑148 (MQ=255)
                                                                                                                                  cccccTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAg            <  2:2567173/148‑1 (MQ=255)
                                                                                                                                    cccTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCa          >  2:2181810/1‑148 (MQ=255)
                                                                                                                                     ccTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCaa         <  1:472348/148‑1 (MQ=255)
                                                                                                                                        gTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCt      >  1:700826/1‑148 (MQ=255)
                                                                                                                                        gTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCt      >  1:59887/1‑148 (MQ=255)
                                                                                                                                            aaTATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTACAGCAACCTGTTg  <  1:360258/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTG  >  CP009273/1116819‑1117106

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: