Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,865,863 T→A 100% T166S (ACC→TCC)  dgoD ← galactonate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,865,8630TA96.7% 736.7 / 9.4 227T166S (ACC→TCC) dgoDgalactonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base A (99/119);  total (101/126)
Fisher's exact test for biased strand distribution p-value = 7.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 228 total aligned reads displayed.

GCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAA  >  CP009273/3865718‑3866005
                                                                                                                                                 |                                                                                                                                              
 ccAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGtt                                                                                                                                             <  2:3108505/148‑1 (MQ=255)
  cAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTa                                                                                                                                            <  1:3000305/148‑1 (MQ=255)
  cAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTa                                                                                                                                            <  1:1227606/148‑1 (MQ=255)
   aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTaa                                                                                                                                           <  2:566725/148‑1 (MQ=255)
   aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTaa                                                                                                                                           <  1:1614649/148‑1 (MQ=255)
    gCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAc                                                                                                                                          >  2:92411/1‑148 (MQ=255)
     cACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAcc                                                                                                                                         >  1:3124881/1‑148 (MQ=255)
         ggCTCCTCAATAAACCGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAAccgccg                                                                                                                                     <  2:1005970/148‑1 (MQ=255)
          gCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACcgccgc                                                                                                                                    >  1:715307/1‑148 (MQ=255)
          gCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACcgccgc                                                                                                                                    >  1:320212/1‑148 (MQ=255)
              ctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCt                                                                                                                                >  2:3513605/1‑148 (MQ=255)
                   tAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACcgc                                                                                                                           >  1:552598/1‑148 (MQ=255)
                    aaaCAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACcgcg                                                                                                                          <  2:535540/148‑1 (MQ=255)
                       cAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCgg                                                                                                                       <  2:2767839/148‑1 (MQ=255)
                       cAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCgg                                                                                                                       >  2:1293777/1‑148 (MQ=255)
                         gCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGa                                                                                                                     >  1:52266/1‑148 (MQ=255)
                             gCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                 >  2:2965156/1‑148 (MQ=255)
                              cGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTAt                                                                                                                >  2:160748/1‑148 (MQ=255)
                                aTACGGCTCCAGGTCTTTATTCAGCACTTTCGCCTTCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCa                                                                                                              <  2:1179028/148‑1 (MQ=255)
                                 tACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCaa                                                                                                             <  1:108821/148‑1 (MQ=255)
                                  aCGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAt                                                                                                            <  2:2656315/148‑1 (MQ=255)
                                  aCGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAt                                                                                                            <  2:1265741/148‑1 (MQ=255)
                                   cGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAtt                                                                                                           >  1:1611997/1‑148 (MQ=255)
                                    ggCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTa                                                                                                          <  1:636807/148‑1 (MQ=255)
                                    ggCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTa                                                                                                          <  1:420715/148‑1 (MQ=255)
                                      cTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGc                                                                                                        <  2:3306866/148‑1 (MQ=255)
                                      cTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGc                                                                                                        >  1:2459385/1‑148 (MQ=255)
                                        ccAGTTATTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGccc                                                                                                      <  2:3139058/148‑1 (MQ=255)
                                         cAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGcccc                                                                                                     >  2:212947/1‑148 (MQ=255)
                                         cAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGcccc                                                                                                     >  2:1059634/1‑148 (MQ=255)
                                             tCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCCATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGtt                                                                                                 <  2:191182/148‑1 (MQ=255)
                                             tCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGtt                                                                                                 <  1:808228/148‑1 (MQ=255)
                                              cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATGTACCGCGCGGGAGTTATCAATTAGCCCCAGttc                                                                                                >  2:1516721/1‑148 (MQ=255)
                                              cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTCCCGCGCGGGAGTTATCAATTAGCCCCAGttc                                                                                                >  2:814223/1‑148 (MQ=255)
                                              cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttc                                                                                                <  2:1155598/148‑1 (MQ=255)
                                                    tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTcaca                                                                                          >  1:1772470/1‑148 (MQ=255)
                                                    tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTcaca                                                                                          <  2:380372/148‑1 (MQ=255)
                                                     cAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACaa                                                                                         <  2:2032717/148‑1 (MQ=255)
                                                      aGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAAc                                                                                        >  1:1931916/1‑148 (MQ=255)
                                                        cACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGTCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCg                                                                                      >  2:475904/1‑148 (MQ=255)
                                                        cACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCg                                                                                      >  1:1253517/1‑148 (MQ=255)
                                                           tttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTc                                                                                   <  1:3075010/148‑1 (MQ=255)
                                                            ttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCa                                                                                  <  1:2116425/148‑1 (MQ=255)
                                                            ttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCa                                                                                  <  1:1188516/148‑1 (MQ=255)
                                                            ttCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACACAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCTCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCa                                                                                  <  1:1023276/148‑1 (MQ=255)
                                                             tCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAg                                                                                 <  1:1328367/148‑1 (MQ=255)
                                                               gCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGtt                                                                               <  1:2397657/148‑1 (MQ=255)
                                                                ccATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGttt                                                                              <  1:986375/148‑1 (MQ=255)
                                                                ccATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGttt                                                                              >  1:897333/1‑148 (MQ=255)
                                                                 cATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTg                                                                             <  2:607757/148‑1 (MQ=255)
                                                                       cgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCGACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTa                                                                       <  2:1601324/148‑1 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGa                                                                   <  2:908025/148‑1 (MQ=255)
                                                                           cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTGTGAAGGTATCGa                                                                   >  2:661871/1‑148 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACACCCGTTCAGTTTGAAGGTATCGAAg                                                                 <  2:2261238/148‑1 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                 <  2:754653/148‑1 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                 <  2:2012519/148‑1 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                 <  2:3401716/148‑1 (MQ=255)
                                                                             gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                 <  2:3155052/148‑1 (MQ=255)
                                                                              aCGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGc                                                                <  1:1450435/148‑1 (MQ=255)
                                                                                gcgACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCg                                                              <  1:1528042/148‑1 (MQ=255)
                                                                                gcgACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCg                                                              <  1:1059634/148‑1 (MQ=255)
                                                                                 cgACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGa                                                             <  1:2546711/148‑1 (MQ=255)
                                                                                   aCCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTCCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAtt                                                           <  2:1039772/148‑1 (MQ=255)
                                                                                     cGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTc                                                         >  1:2949831/1‑148 (MQ=255)
                                                                                       tGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTcgc                                                       <  1:2774712/148‑1 (MQ=255)
                                                                                         gAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGt                                                     >  2:2490913/1‑148 (MQ=255)
                                                                                          aaaTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTa                                                    >  2:1944454/1‑148 (MQ=255)
                                                                                            aTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGc                                                  >  2:384276/1‑148 (MQ=255)
                                                                                               aaGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGtt                                               <  2:2250940/148‑1 (MQ=255)
                                                                                                aGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTCCCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGttt                                              >  2:2437630/1‑148 (MQ=255)
                                                                                                  aCCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTa                                            <  2:298675/148‑1 (MQ=255)
                                                                                                   ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTaa                                           <  1:3188687/148‑1 (MQ=255)
                                                                                                   ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTaa                                           <  1:2007772/148‑1 (MQ=255)
                                                                                                    cAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAAt                                          <  2:2459385/148‑1 (MQ=255)
                                                                                                      aaCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTCCCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGc                                        >  2:70637/1‑148 (MQ=255)
                                                                                                         tCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGt                                     >  1:1799378/1‑148 (MQ=255)
                                                                                                         tCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGt                                     >  1:1492904/1‑148 (MQ=255)
                                                                                                         tCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGt                                     >  2:1419951/1‑148 (MQ=255)
                                                                                                         tCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTAAATGCCGt                                     <  1:1349314/148‑1 (MQ=255)
                                                                                                         gcAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATGTCGCGTAGCGTGTTAATGCCGt                                     >  1:1799675/2‑148 (MQ=255)
                                                                                                              cTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATa                                <  2:262132/148‑1 (MQ=255)
                                                                                                                gATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAAc                              <  1:2762301/148‑1 (MQ=255)
                                                                                                                gATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGATCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAAc                              >  2:2352364/1‑148 (MQ=255)
                                                                                                                 aTTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACa                             >  2:1776799/1‑148 (MQ=255)
                                                                                                                      cAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgc                        <  2:552598/148‑1 (MQ=255)
                                                                                                                      cAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgc                        <  1:92439/148‑1 (MQ=255)
                                                                                                                       aaaaGCTTCACGAATTTGTGCCACGGGGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgcc                       >  1:590649/1‑148 (MQ=255)
                                                                                                                          aGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGa                    <  1:2606243/148‑1 (MQ=255)
                                                                                                                           gCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGAc                   >  1:2222214/1‑148 (MQ=255)
                                                                                                                           gCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGAc                   >  1:983541/1‑148 (MQ=255)
                                                                                                                             ttCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGa                 <  1:1038810/148‑1 (MQ=255)
                                                                                                                             ttCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGa                 <  2:3319827/148‑1 (MQ=255)
                                                                                                                             ttCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGa                 <  2:1557231/148‑1 (MQ=255)
                                                                                                                              tCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGAt                >  1:581811/1‑148 (MQ=255)
                                                                                                                                aCGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcg              >  2:1653931/1‑148 (MQ=255)
                                                                                                                                 cGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgc             <  2:2724817/148‑1 (MQ=255)
                                                                                                                                   aaTTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccg           >  2:157635/1‑148 (MQ=255)
                                                                                                                                     tttGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCgccgcc         <  2:2262979/148‑1 (MQ=255)
                                                                                                                                      ttGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGGTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGcca        <  2:3206207/148‑1 (MQ=255)
                                                                                                                                         tgCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaacc     <  2:1164812/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAA  >  CP009273/3865718‑3866005

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: