Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,897,666 G→A intergenic (‑5444/‑3153) RPPX_08470 ← / → RPPX_08480 tryptophan synthase subunit beta/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,897,6660GA77.8% 2878.8 / inf 2514intergenic (‑5444/‑3153)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (299/256);  new base A (992/952);  total (1298/1216)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-07

Only 100 of 2519 total aligned reads displayed.

CGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACGTCTACCAGGACGCAGGC  >  CP009974/1897516‑1897816
                                                                                                                                                      |                                                                                                                                                      
  ctgATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAcc                                                                                                                                                      >  2:373109/3‑151 (MQ=21)
  cAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAAcc                                                                                                                                                      <  1:1286139/151‑1 (MQ=18)
   tgATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAccg                                                                                                                                                     <  1:2054454/150‑2 (MQ=18)
    gATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAAccga                                                                                                                                                    >  1:4077/1‑149 (MQ=25)
    gATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAAccga                                                                                                                                                    <  1:3286144/151‑3 (MQ=18)
       gCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCt                                                                                                                                                 >  1:2888496/1‑151 (MQ=9)
       gCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCt                                                                                                                                                 <  1:2202836/151‑1 (MQ=21)
        cATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCtt                                                                                                                                                >  2:1025747/1‑151 (MQ=18)
         aTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTc                                                                                                                                               <  1:4730838/151‑1 (MQ=9)
           gcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCAc                                                                                                                                             <  2:5342149‑M2/105‑1 (MQ=255)
            gAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCAcc                                                                                                                                            <  2:4411311/151‑1 (MQ=18)
            gAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCAcc                                                                                                                                            >  2:2913503/1‑151 (MQ=18)
            gAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCAct                                                                                                                                            >  1:5608787/1‑150 (MQ=9)
            gAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCAcc                                                                                                                                            <  1:4371112/151‑1 (MQ=9)
            gAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCAcc                                                                                                                                            >  1:4536253/1‑151 (MQ=18)
             aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCg                                                                                                                                           <  1:2768692/151‑1 (MQ=17)
             aaGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCg                                                                                                                                           >  2:2163755/1‑151 (MQ=17)
               gTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTg                                                                                                                                         >  1:3597665/1‑151 (MQ=14)
               gTTCGAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTg                                                                                                                                         <  2:1233002/151‑1 (MQ=18)
                ttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGc                                                                                                                                        >  1:1802508‑M2/42‑151 (MQ=255)
                ttCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGc                                                                                                                                        >  1:1533124/1‑151 (MQ=18)
                ttCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGc                                                                                                                                        <  2:2525160/151‑1 (MQ=18)
                  cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTCCAc                                                                                                                                      >  2:5969043/1‑151 (MQ=17)
                  cGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCAc                                                                                                                                      >  1:5423286/1‑151 (MQ=18)
                   gAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGTAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACg                                                                                                                                     >  1:4845203/1‑151 (MQ=17)
                   gAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACg                                                                                                                                     <  2:3825957/151‑1 (MQ=2)
                   gAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCAACCTTCACTGTGCACg                                                                                                                                     <  2:6398933/151‑1 (MQ=11)
                    aacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGcacat                                                                                                                                    <  1:6544831‑M2/114‑3 (MQ=255)
                    aaaaCCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGt                                                                                                                                    >  1:6209937/1‑151 (MQ=2)
                       acctgaacctgGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCaa                                                                                                                                 <  2:2688514/151‑1 (MQ=17)
                       acctgaacctgGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTCCACCGTGCACGTCaa                                                                                                                                 >  2:5549959/1‑151 (MQ=17)
                         ctgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAAcc                                                                                                                               >  2:5829431/1‑151 (MQ=17)
                          tgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCACATCaacac                                                                                                                              <  2:2043234‑M2/120‑3 (MQ=255)
                          tgaacctgGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCa                                                                                                                              >  2:2305764/1‑151 (MQ=2)
                           gaacctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAg                                                                                                                             >  2:5713192/1‑151 (MQ=14)
                            acgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCACAtcaacaccg                                                                                                                            >  1:3476945‑M2/30‑147 (MQ=255)
                             acctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGcc                                                                                                                           >  1:6319294/1‑151 (MQ=14)
                             accggGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGcc                                                                                                                           <  2:1542623/151‑1 (MQ=17)
                              cctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCg                                                                                                                          <  2:3477668/151‑1 (MQ=14)
                                 ggACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTg                                                                                                                       <  2:519513/151‑1 (MQ=11)
                                       gaccccggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCCCCGTCCACGTCAACCCGCCGCTCGCCCAc                                                                                                                 >  1:6438247/8‑151 (MQ=17)
                                       aacccccgTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                                                                                                 >  1:1638073/8‑151 (MQ=17)
                                        accccggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGc                                                                                                                            >  2:628970/7‑139 (MQ=17)
                                         cTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                                                                                               >  2:1911098/1‑151 (MQ=18)
                                            cAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGActt                                                                                                            <  1:3379975/151‑3 (MQ=18)
                                              gTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTga                                                                                                          >  2:6535516/1‑150 (MQ=17)
                                              gTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTga                                                                                                          <  2:744105/151‑2 (MQ=17)
                                                caacACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTc                                                                                                        >  2:5104915/4‑151 (MQ=11)
                                                cAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTc                                                                                                        >  2:872751/1‑151 (MQ=11)
                                                    aCCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGa                                                                                                    >  1:5855603/1‑151 (MQ=17)
                                                       gAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGAccc                                                                                                 <  2:3636301/151‑1 (MQ=11)
                                                       gAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGGCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAccc                                                                                                 <  1:5291561/151‑1 (MQ=11)
                                                         ggTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTg                                                                                               >  1:3712296/1‑151 (MQ=14)
                                                         ggTAACCGGCGAGCCCGGCACCCCGGGGAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTg                                                                                               <  2:5618380/151‑1 (MQ=14)
                                                          gTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCCACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGa                                                                                              >  1:3837035/1‑151 (MQ=14)
                                                           tGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAg                                                                                             <  1:3126136/151‑1 (MQ=14)
                                                             aCCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGATCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGca                                                                                           <  1:552548/151‑1 (MQ=14)
                                                             aCCGACGAGCCAGGCACCCCGGGAAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGca                                                                                           <  2:2809634/151‑1 (MQ=14)
                                                              ccGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaa                                                                                          <  1:2298461/151‑1 (MQ=14)
                                                               cgacgaGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaac                                                                                         >  1:1165050/1‑151 (MQ=11)
                                                                gacgaGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGcaaca                                                                                        >  2:4643778/1‑151 (MQ=14)
                                                                gacgaGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaaca                                                                                        <  1:5079469/151‑1 (MQ=255)
                                                                  cgaGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCaacaac                                                                                      >  1:758199/1‑151 (MQ=17)
                                                                  cgaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaat                                                                                      <  1:1006142/151‑2 (MQ=14)
                                                                   gaGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACg                                                                                     <  2:3033609/151‑1 (MQ=17)
                                                                    aGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaacaatgc                                                                                    >  1:461772/1‑148 (MQ=14)
                                                                     gCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCCCCGTGCACGTCACCCAGCCGCCGGCCCACGACCTGGACGTGACCCTGAGcaacaatccc                                                                                   >  2:6181956/1‑147 (MQ=11)
                                                                              cccaggCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCAcc                                                                          >  2:1420223/5‑151 (MQ=11)
                                                                              cccGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCAcc                                                                          <  1:5992418/150‑1 (MQ=14)
                                                                               ccGGGTGACGAAGGCGAGATCATCAAGGCCACCCTCACGGCCGACCAGCCCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCa                                                                         <  2:1166681/151‑1 (MQ=14)
                                                                               ccGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCCCAGGTCACCa                                                                         >  1:5731347/1‑151 (MQ=11)
                                                                                cGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCAt                                                                        <  1:3988819/151‑1 (MQ=14)
                                                                                  ggCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCa                                                                      >  2:4093898/1‑151 (MQ=11)
                                                                                  ggCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                                                                      >  2:4747101/1‑151 (MQ=14)
                                                                                   gCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTg                                                                                                                                         >  2:5009777/1‑83 (MQ=14)
                                                                                   gCAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGCGCAACAATGCTCAGGTCACCATCaa                                                                     <  1:2549577/151‑1 (MQ=11)
                                                                                     aaCGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAgg                                                                   >  2:1784981/1‑151 (MQ=255)
                                                                                        gAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGctg                                                                >  2:3179792/1‑149 (MQ=255)
                                                                                        gAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGctg                                                                <  1:859725/151‑3 (MQ=11)
                                                                                           ggCGATATCATCAAGGTCACCATCACGGCCGCCCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTggcg                                                             <  1:1399822/151‑3 (MQ=11)
                                                                                           ggCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGCTggcg                                                             <  1:4952045/151‑3 (MQ=11)
                                                                                            gCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGa                                                            <  2:5475977/151‑1 (MQ=11)
                                                                                            gCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCa                                                                      <  1:927039/141‑1 (MQ=11)
                                                                                            gCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAATGCTCAGGTCACCATCAAGGctggcga                                                            >  1:1929128/1‑145 (MQ=11)
                                                                                             cGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGaa                                                           >  1:3011526/1‑151 (MQ=14)
                                                                                             cGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGCTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGaa                                                           >  2:4878005/1‑151 (MQ=14)
                                                                                             cGATATCATCAAGGGCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACATTGACCGTGCACGTCAACCAGCCGCTGGCCCACGACCCGATGGTGAGCCTGAGCAACAATGCGCAGGTCACCATCGAAGCCgccgaa                                                           >  2:3341049/1‑146 (MQ=11)
                                                                                              gacttggtcaAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGaaa                                                          >  2:6089264/8‑151 (MQ=11)
                                                                                                     tcaAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAgtgc                                                   <  1:264516/151‑1 (MQ=255)
                                                                                                        aGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGCCCAAACCAGCGCGcc                                                >  2:3999362/1‑151 (MQ=11)
                                                                                                           tCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGCGCAACAATGCTCAGGTCACCATCAACGCTGGCGAAACCAGCGCGCCGTa                                             >  2:4172428/1‑151 (MQ=14)
                                                                                                            cACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTac                                            >  1:6084312/1‑151 (MQ=11)
                                                                                                                aTCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCTGGCGAAACCAGTGCGCCGTACAccc                                        <  2:1643137/151‑1 (MQ=11)
                                                                                                                    ccgCCGACCAGACCTCGGGGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCacga                                    <  2:4286203/149‑1 (MQ=11)
                                                                                                                     cgCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCGTACACCCacgat                                   >  1:423215/2‑150 (MQ=11)
                                                                                                                            cAGACTTCGGTGGCCGAGAACGTCACGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGCCTTGGTCGTGACCCTGCGCACCGACGCCCAGGTCACCAGCAAGGCCGGTGAAACCAGTGCGCCGTCCCCCCACGACGcgccag                            >  2:4515237/1‑146 (MQ=11)
                                                                                                                                           gagaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGt                                                                                                      >  1:2410733/1‑62 (MQ=255)
                                                                                                                                           gagaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGt                                                                                                      <  2:2410733/62‑1 (MQ=255)
                                                                                                                                            agaACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGATCAACAACGCCCAGGTCACCATCAACGCCGGTGAAACCAGCGCGCCGTAAACCCACCACGCGCAAGGCGATGACggcaacaa            >  2:3555894/1‑144 (MQ=11)
                                                                                                                                                 gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTAcaacgatg       >  1:5169137/1‑144 (MQ=11)
                                                                                                                                                      |                                                                                                                                                      
CGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCACAGGGCGATGACGTCTACCAGGACGCAGGC  >  CP009974/1897516‑1897816

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