Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,498,206 C→T 6.7% intergenic (+64/+88) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,498,2060CT 6.7% 660.8 / 25.2 209intergenic (+64/+88)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base C (73/121);  new base T (10/4);  total (84/125)
Fisher's exact test for biased strand distribution p-value = 2.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 210 total aligned reads displayed.

GAATGACAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCG  >  CP009273/3498061‑3498355
                                                                                                                                                 |                                                                                                                                                      
gAATGACAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGt                                                                                                                                                    >  1:4707938/1‑150 (MQ=255)
 aaTGACAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTa                                                                                                                                                   >  2:2320790/1‑150 (MQ=255)
 aaTGACAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTa                                                                                                                                                   >  1:3482025/1‑150 (MQ=255)
  aTGACCCACCGGTGCACGTCTCCGACCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTAc                                                                                                                                                  <  2:5173284/150‑1 (MQ=255)
   tGACAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGGTGAATTTTACATCCCGTACg                                                                                                                                                 >  1:3095486/1‑150 (MQ=255)
      cAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGtgg                                                                                                                                              >  1:2101888/1‑148 (MQ=255)
       aaaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACgggga                                                                                                                                             >  1:522338/1‑146 (MQ=255)
        aaCCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACggtggc                                                                                                                                            >  1:506420/1‑145 (MQ=255)
           cGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTc                                                                                                                                        <  1:3218616/150‑1 (MQ=255)
           cGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTAACATCCCGTACGTTCCCCTc                                                                                                                                        <  2:2775212/150‑1 (MQ=255)
           cGGTGCACGTCACCGAACTCTTCTGCCGCGCCAATCGCATCACCTTACCTAGTGATAATAACAGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTAACATCCCGAACGTTCCCCTc                                                                                                                                        <  2:3329681/150‑1 (MQ=39)
             gTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTCCGTTGCCCGCAc                                                                                                                                      >  1:4755248/1‑150 (MQ=255)
                cACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCcccca                                                                                                                                   >  1:2769452/1‑149 (MQ=255)
                cACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCt                                                                                                                                   <  1:2140210/150‑1 (MQ=255)
                cACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCt                                                                                                                                   <  1:2757023/150‑1 (MQ=255)
                 aCGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTa                                                                                                                                  <  2:4602687/150‑1 (MQ=255)
                 aCGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTa                                                                                                                                  <  2:4602641/150‑1 (MQ=255)
                   gTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCCAAc                                                                                                                                >  1:4086587/1‑150 (MQ=255)
                   gTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTGCCCACCCCCaaaa                                                                                                                                >  1:5043384/1‑149 (MQ=39)
                    tctcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAcc                                                                                                                               >  2:273895/1‑150 (MQ=255)
                      tcCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACCCCCAAAcccc                                                                                                                             >  1:3776165/1‑149 (MQ=255)
                        cGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCCCACCCAAAcccaca                                                                                                                           >  2:2801123/1‑147 (MQ=255)
                        cGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATACCGTACGTTTCCCACACCCaaaaaccaa                                                                                                                           >  2:1766344/1‑144 (MQ=39)
                         gAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCtctc                                                                                                                          <  2:3256145/150‑1 (MQ=255)
                          aaCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCAAAaccacaac                                                                                                                         >  1:2524570/1‑145 (MQ=255)
                          aaCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcc                                                                                                                         <  1:4449467/150‑1 (MQ=255)
                          aaCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCACCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTGGCCCCCCCCCCAAAcccacac                                                                                                                          >  1:968952/1‑146 (MQ=38)
                           aCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGCTACTAAAAGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTccc                                                                                                                        <  2:3110982/150‑1 (MQ=255)
                            ctctTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTcccc                                                                                                                       <  2:3625458/150‑1 (MQ=255)
                            ctctTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAGAcccacaaccc                                                                                                                       >  2:466599/1‑143 (MQ=39)
                            ctctTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCGCACCCAAACCCACTaccc                                                                                                                       >  2:1128773/1‑150 (MQ=39)
                            ctctTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTGCCCGCACCCCAACCCACAcccc                                                                                                                       >  2:2014124/1‑150 (MQ=39)
                                 ctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATTCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgg                                                                                                                  >  2:1307644/1‑150 (MQ=255)
                                  tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                 <  1:2495980/150‑1 (MQ=255)
                                  tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                 <  1:5075275/150‑1 (MQ=255)
                                  tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCGAACCCGCACCCCAAAggg                                                                                                                 >  2:585576/1‑150 (MQ=255)
                                         ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                          <  2:2391522/150‑1 (MQ=255)
                                         ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCTCACCCCCAAGGGGCGCggg                                                                                                          >  2:4031522/1‑150 (MQ=39)
                                           aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCCCCCCAAACCACAACCCCAACGGGCGAggggg                                                                                                        >  1:1743993/1‑149 (MQ=37)
                                           aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCACAACCCCAAAAGGCGCGGGGa                                                                                                        >  1:3348124/1‑150 (MQ=38)
                                              cGCATCACCTTAACTATTGATAATAAAAGACCTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCCCCCGAACCCTCTCCCCAAAGGGGCGAGGGGACCg                                                                                                     <  2:863266/150‑1 (MQ=255)
                                              cGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAACCCACACCCCCCAGGGGCGAGGGGGccc                                                                                                     >  1:767121/1‑149 (MQ=38)
                                              cGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCg                                                                                                     <  1:3362321/150‑1 (MQ=255)
                                                cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCggt                                                                                                   >  2:2944048/1‑148 (MQ=255)
                                                 aTCACCGTAACTATTGAGAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCCCCCTAACCCTCTCCCCACCGGGGCGAGGGGACCGTTc                                                                                                  <  2:4511834/150‑1 (MQ=39)
                                                    aCCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAACCCAAAACCCAAAGGGGCGAGGGGACCggtgcgt                                                                                               >  2:2109023/1‑144 (MQ=38)
                                                     ccTGAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGAGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACCTCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                              <  2:2589518/150‑1 (MQ=255)
                                                      cTTAACTATTGATAATAAAAGACATGAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                             <  2:5086826/150‑1 (MQ=255)
                                                          aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCCAACCCACTCCCCCAAGGGGCGCGGGGACCGGTCACTACGGac                                                                                         >  2:4066287/1‑149 (MQ=38)
                                                            tcttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                       <  1:3536698/148‑1 (MQ=255)
                                                            tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                       <  2:3453799/150‑1 (MQ=255)
                                                             aTTGATAAGAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACAGCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAttt                                                                                      <  2:2288412/150‑1 (MQ=255)
                                                              ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                     <  1:3558442/150‑1 (MQ=255)
                                                               tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATTCCGTACGTTCCCCTCCCCCTAACCCACTCCCCAAAGGGGCGAGGGGACCGTTCCGTACGGAATTTAg                                                                                    >  2:674155/1‑150 (MQ=255)
                                                               tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATTCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAATGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                    >  2:667412/1‑150 (MQ=255)
                                                                  taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCtt                                                                                 <  1:3707320/150‑1 (MQ=255)
                                                                  taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCtt                                                                                 <  1:2336911/150‑1 (MQ=255)
                                                                  taataaAAGACATGAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGTGCGTTCAGTACGGAATTTAGCtt                                                                                 <  2:4404346/150‑1 (MQ=255)
                                                                     taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTCACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                              <  2:378621/150‑1 (MQ=255)
                                                                     taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                              <  2:378368/150‑1 (MQ=255)
                                                                      aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACg                                                                             <  1:1545517/150‑1 (MQ=255)
                                                                        aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                           <  1:561300/150‑1 (MQ=255)
                                                                        aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                           <  1:2377205/150‑1 (MQ=255)
                                                                         aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCCTCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                          <  2:2130038/150‑1 (MQ=255)
                                                                         aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCGAAACCACACCCCCAAGGGGCGGGGGGACCGTTCAGTACCGCATTTTTCTTGACgggt                                                                          >  2:4606791/1‑148 (MQ=37)
                                                                             aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAt                                                                      <  1:2014124/150‑1 (MQ=255)
                                                                             aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACAACCAAAACCACACCCCAAATGGGGGGGGGGCCCGGTCATTAAGGAAATTAGATTGTAGGTGCTaa                                                                      >  1:4325595/1‑149 (MQ=21)
                                                                                aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATTCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                   <  1:1278237/150‑1 (MQ=255)
                                                                                aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                   <  2:4408521/150‑1 (MQ=255)
                                                                                aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCCCACCCAAACCCACAACCCCAACGGGCGCGGGGCCCGGTCACTAAGGAATTTATATTGACGGTGCTATCCa                                                                   >  1:4276247/1‑150 (MQ=25)
                                                                                aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCACACCCAAAACCACAACCCCAAAGGGCCAGGGGACCCTTCACTACCGAAATTATCATGACGGTGCTATACa                                                                   >  1:4734520/1‑150 (MQ=25)
                                                                                 aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATTCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                  <  2:2662854/150‑1 (MQ=255)
                                                                                    gTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAt                                                               <  2:2672764/150‑1 (MQ=255)
                                                                                    gTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAt                                                               <  1:4954/150‑1 (MQ=255)
                                                                                         ccGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAAt                                                          <  2:4458266/150‑1 (MQ=255)
                                                                                         ccGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACCACCAAAACCACAACCCCAAGGGGCGGGGGGTCCGGTCAGTTAGGAAATTATCTTTTCGGTGCTATCCCGCCTTaaaa                                                          >  1:4724055/1‑149 (MQ=255)
                                                                                         ccGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGCCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAACGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAAt                                                          <  2:3248150/150‑1 (MQ=255)
                                                                                           ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAc                                                        <  1:2975946/150‑1 (MQ=255)
                                                                                           ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAc                                                        <  2:2939353/150‑1 (MQ=255)
                                                                                             aTAAGATGCGCCAGCATCGCCTCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                      <  2:4136731/150‑1 (MQ=255)
                                                                                             aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGATGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                      >  1:2091394/1‑150 (MQ=255)
                                                                                              tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                     <  1:3794122/150‑1 (MQ=255)
                                                                                              tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAAAGTAGCTTGACGGGGCTATCCAGCAGTAAAAACggg                                                     >  2:1253413/1‑150 (MQ=38)
                                                                                               aaGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCGCACCCAAAACCGCTACCCAAAGGGGCGAGGGGACCGTTCAGTACCGAATTTTTCTTGACGGTGCGATCCAACATTAAATACGGGc                                                    >  2:1691764/1‑150 (MQ=255)
                                                                                                aGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCa                                                   <  1:3640163/150‑1 (MQ=255)
                                                                                                  aTGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATTCCGTACGTTCCCCTCACCCAAACCCTCTACCCCAAGGGGCGAGGGGACCGGTCAGTACCGAAATTAGCGTGTCGGTGCTATCCCGCATTAAATACTGGCAAt                                                 >  2:1377019/1‑150 (MQ=255)
                                                                                                   tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAATGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATa                                                <  2:9123/150‑1 (MQ=255)
                                                                                                      gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                             <  2:3824639/150‑1 (MQ=255)
                                                                                                       ccAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActct                                            <  2:3294202/150‑1 (MQ=255)
                                                                                                         aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                          <  2:1913919/150‑1 (MQ=255)
                                                                                                          gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                         <  1:2452489/150‑1 (MQ=255)
                                                                                                            atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTACCCACACACAAACCCACACCCCAAATGGGCGAGGGGACCGTTAAGTACGGAATTTAGCTTGACGGTGCGATAAAGCAAAAAAAACAGGCAATACTCTAACGc                                       >  2:3199392/1‑150 (MQ=255)
                                                                                                                         aTGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcagca                          <  1:1724979/150‑1 (MQ=255)
                                                                                                                         aTGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcagca                          <  1:1131324/150‑1 (MQ=255)
                                                                                                                           gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                        <  1:2068948/150‑1 (MQ=255)
                                                                                                                                  cgTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAAACCACACCCCAAAGGGGGGAGGGGGCCGGTCAGTAAGGCATTTATATTGACGGTGGGATACAGCCATAAAAAAATGCAATAAACTCACCCAAAAAAACCACCACCGACTCtt                 >  2:4919382/1‑150 (MQ=255)
                                                                                                                                    ttGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCa               <  1:5066917/150‑1 (MQ=255)
                                                                                                                                        aTTTTACATCCCGTACGTTCCCCTCACCCGAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGc           <  2:869651/150‑1 (MQ=255)
                                                                                                                                          tttACATTCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg         <  1:4577925/150‑1 (MQ=255)
                                                                                                                                                tCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc   <  1:2932449/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
GAATGACAAACCGGTGCACGTCTCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCG  >  CP009273/3498061‑3498355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: