Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,431 C→A 100% P90T (CCA→ACA)  yiaN → 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,738,4310CA99.5% 607.4 / ‑4.0 186P90T (CCA→ACA) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (103/82);  total (104/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-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.

GATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGG  >  CP009273/3738286‑3738575
                                                                                                                                                 |                                                                                                                                                
 aTCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                               >  2:3273862/1‑148 (MQ=255)
 aTCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAg                                                                                                                                               >  2:1056973/1‑148 (MQ=255)
  tCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAgg                                                                                                                                              >  2:3087278/1‑148 (MQ=255)
            aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggc                                                                                                                                    <  1:1361304/148‑1 (MQ=255)
            aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggc                                                                                                                                    <  2:2627658/148‑1 (MQ=255)
            aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggc                                                                                                                                    <  1:52115/148‑1 (MQ=255)
             aCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCgggctgggct                                                                                                                                   >  1:1065466/1‑148 (MQ=255)
              cGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCGTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCGGGGCTa                                                                                                                                  >  1:886047/1‑148 (MQ=255)
                cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACg                                                                                                                                <  2:3263708/148‑1 (MQ=255)
                cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACg                                                                                                                                <  2:1521946/148‑1 (MQ=255)
                 tGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGt                                                                                                                               <  2:2893862/148‑1 (MQ=255)
                  ggTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTc                                                                                                                              <  2:1447235/148‑1 (MQ=255)
                    tGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCCAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCgg                                                                                                                            >  2:2894307/1‑148 (MQ=255)
                      aaCGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCg                                                                                                                          >  1:2102266/1‑148 (MQ=255)
                        cGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                        <  1:3028425/148‑1 (MQ=255)
                        cGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                        <  1:2347319/148‑1 (MQ=255)
                        cGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                        <  1:25928/148‑1 (MQ=255)
                        cGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTg                                                                                                                        <  1:1777262/148‑1 (MQ=255)
                           aGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTg                                                                                                                     >  1:1090755/1‑148 (MQ=255)
                            gCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTgg                                                                                                                    <  1:1734947/148‑1 (MQ=255)
                               gATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAg                                                                                                                 >  2:2009154/1‑148 (MQ=255)
                                 tAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                               >  2:2107937/1‑148 (MQ=255)
                                  aGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTa                                                                                                              >  2:2863467/1‑148 (MQ=255)
                                    cTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATg                                                                                                            <  1:2242316/148‑1 (MQ=255)
                                     ttCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGa                                                                                                           <  1:936453/148‑1 (MQ=255)
                                      tctcCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGAt                                                                                                          <  1:3281623/148‑1 (MQ=255)
                                        tcCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTa                                                                                                        <  2:884745/148‑1 (MQ=255)
                                         cccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCCTAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTAt                                                                                                       >  1:3168468/1‑148 (MQ=255)
                                          ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATg                                                                                                      >  2:1005018/1‑148 (MQ=255)
                                          ccTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATg                                                                                                      >  2:636632/1‑148 (MQ=255)
                                            tgctgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGc                                                                                                    >  2:1649492/1‑148 (MQ=255)
                                               tgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                 <  1:2252895/148‑1 (MQ=255)
                                               tgGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAg                                                                                                 <  1:591538/148‑1 (MQ=255)
                                                 gCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGcc                                                                                               <  1:2248339/148‑1 (MQ=255)
                                                    aTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCttt                                                                                            <  2:1464671/148‑1 (MQ=255)
                                                    aTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCttt                                                                                            <  2:319222/148‑1 (MQ=255)
                                                    aTTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCttt                                                                                            <  1:1807697/148‑1 (MQ=255)
                                                     ttCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTc                                                                                           <  2:66952/148‑1 (MQ=255)
                                                      tCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCt                                                                                          >  1:626373/1‑148 (MQ=255)
                                                        cGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTgg                                                                                        <  2:2155042/148‑1 (MQ=255)
                                                         gTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGt                                                                                       <  1:603262/148‑1 (MQ=255)
                                                          ttCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                      >  1:1460804/1‑148 (MQ=255)
                                                          ttCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                      >  2:2427551/1‑148 (MQ=255)
                                                          ttCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGtt                                                                                      >  1:2883036/1‑148 (MQ=255)
                                                             tttGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTNGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTg                                                                                   <  1:1678280/148‑1 (MQ=255)
                                                             tttGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTg                                                                                   <  2:1398587/148‑1 (MQ=255)
                                                              ttGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGc                                                                                  >  2:2038712/1‑148 (MQ=255)
                                                              ttGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGc                                                                                  >  2:1296200/1‑148 (MQ=255)
                                                                   cTGGAGGGTGAAATCATGAATGCGGGAGGCTTGTCAAAACGCATTGTTGACCTTCAGATTAAACTGGTGGGGCATAAAACAGGCGGGATGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTg                                                                             <  2:2081348/148‑1 (MQ=255)
                                                                       cGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGa                                                                         <  2:2131874/148‑1 (MQ=255)
                                                                        gggTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGAt                                                                        <  1:335594/148‑1 (MQ=255)
                                                                        gggTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGAt                                                                        <  1:2002712/148‑1 (MQ=255)
                                                                                 cATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgcc                                                               >  2:603436/1‑148 (MQ=255)
                                                                                  aTGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATAccgccgccg                                                              <  2:43250/148‑1 (MQ=255)
                                                                                   tGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGt                                                             <  2:2927034/148‑1 (MQ=255)
                                                                                    gAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGtt                                                            >  1:279771/1‑148 (MQ=255)
                                                                                     aaTGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTCTGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTg                                                           <  2:2175254/148‑1 (MQ=255)
                                                                                     aaTGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTg                                                           >  1:1997160/1‑148 (MQ=255)
                                                                                           ggcggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGctc                                                     >  2:1741546/1‑148 (MQ=255)
                                                                                             cggcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTa                                                   <  1:1850695/148‑1 (MQ=255)
                                                                                              ggcTTGTCAAATCGCATTGTTGACCTGCCGATTAAACTGGTGGGGCATATAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTGTGCCGTTGCGGATACCGCCGCCGTTGCGGctcttc                                                  <  2:2886647/148‑3 (MQ=255)
                                                                                               gcTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACt                                                 >  1:581985/1‑148 (MQ=255)
                                                                                                cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                >  2:2741485/1‑148 (MQ=255)
                                                                                                cTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTg                                                >  1:790764/1‑148 (MQ=255)
                                                                                                 ttGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTgg                                               >  1:3303317/1‑148 (MQ=255)
                                                                                                     cAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGcc                                           >  2:2529576/1‑148 (MQ=255)
                                                                                                       aaaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCga                                         >  1:2809583/1‑148 (MQ=255)
                                                                                                        aaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgat                                        >  2:2099569/1‑148 (MQ=255)
                                                                                                        aaCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgat                                        >  2:1177371/1‑148 (MQ=255)
                                                                                                         aCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                       >  2:2088701/1‑148 (MQ=255)
                                                                                                         aCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                       <  1:822531/148‑1 (MQ=255)
                                                                                                         aCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                       >  1:2307895/1‑148 (MQ=255)
                                                                                                         aCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                       >  1:2288145/1‑148 (MQ=255)
                                                                                                         aCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCgatg                                       >  2:2552117/1‑148 (MQ=255)
                                                                                                            cATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGatgatg                                    <  2:1503067/148‑1 (MQ=255)
                                                                                                            cATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGatgatg                                    <  1:1056973/148‑1 (MQ=255)
                                                                                                              ttgttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcg                                  >  2:478162/1‑148 (MQ=255)
                                                                                                              ttgttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcg                                  >  2:1686623/1‑148 (MQ=255)
                                                                                                               tgttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcgc                                 >  2:1306334/1‑148 (MQ=255)
                                                                                                               tgttgACCTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATgcgc                                 <  1:2952971/148‑1 (MQ=255)
                                                                                                                     ccTGCCGATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGcc                           >  2:253908/1‑148 (MQ=255)
                                                                                                                           gATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTAt                     <  1:2167413/148‑1 (MQ=255)
                                                                                                                           gATGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTAt                     <  2:1084394/148‑1 (MQ=255)
                                                                                                                            aTGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                    >  2:3174855/1‑148 (MQ=255)
                                                                                                                            aTGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                    >  2:1782408/1‑148 (MQ=255)
                                                                                                                            aTGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                    <  2:2315395/148‑1 (MQ=255)
                                                                                                                            aTGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                    >  1:2212294/1‑148 (MQ=255)
                                                                                                                            aTGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                    >  1:2213471/1‑148 (MQ=255)
                                                                                                                            aTGAAACTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATc                    >  1:679555/1‑148 (MQ=255)
                                                                                                                               aaaCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCgg                 >  1:2072549/1‑148 (MQ=255)
                                                                                                                               aaaCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCgg                 >  1:2073035/1‑148 (MQ=255)
                                                                                                                                 aCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTc               >  1:1342274/1‑148 (MQ=255)
                                                                                                                                 aCTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTc               >  1:3263376/1‑148 (MQ=255)
                                                                                                                                  cTGGTGGGGCATAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCa              >  2:3093669/1‑148 (MQ=255)
                                                                                                                                             tAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg   >  2:3124710/1‑148 (MQ=255)
                                                                                                                                             tAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg   <  2:2295630/148‑1 (MQ=255)
                                                                                                                                             tAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg   <  1:1008403/148‑1 (MQ=255)
                                                                                                                                             tAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg   >  1:3294838/1‑148 (MQ=255)
                                                                                                                                             tAAAACAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGgcggcg   <  1:1005018/148‑1 (MQ=255)
                                                                                                                                              aaaaCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGcggcgg  >  2:814498/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTCTTTGTGCTGGCGGGTGAAATCATGAATGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGG  >  CP009273/3738286‑3738575

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: ‑