Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,897,594 G→C intergenic (‑5372/‑3225) 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,5940GC76.6% 5402.0 / inf 4542intergenic (‑5372/‑3225)RPPX_08470/RPPX_08480tryptophan synthase subunit beta/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (430/471);  new base C (1519/1438);  total (2315/2227)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 4561 total aligned reads displayed.

CAACCGCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGT  >  CP009974/1897446‑1897740
                                                                                                                                                    |                                                                                                                                                  
    tgcaccagACAACGTCTATGTCGGCCCCAATGCGCCGGTGGTCAATGCAATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACATGGACAAGACCCCGGTCAGCACCGAGGGGACCGACGAGCCAGGAACGCCGGGc                                                                                                                                              >  1:812147/8‑151 (MQ=255)
    tgCCCCGGACAACGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGc                                                                                                                                              <  1:3659764/150‑1 (MQ=255)
     gCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCa                                                                                                                                             >  1:2786164/1‑151 (MQ=255)
        ccGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACg                                                                                                                                          >  1:720110/1‑151 (MQ=255)
        ccGGACAACGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAATCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACg                                                                                                                                          >  2:5201018/1‑151 (MQ=25)
          ggACAACGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGaa                                                                                                                                        >  2:2161267/1‑151 (MQ=37)
          agACAACGTCTATGTCGGCACCAACGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGaa                                                                                                                                        <  2:5307284/150‑1 (MQ=255)
           gACAACGTCTATGTCGGCACCAATGCGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAg                                                                                                                                       <  1:3199750/151‑1 (MQ=255)
            aCAACGTCTATGTCGGTACCAACGCGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAgg                                                                                                                                      >  2:4448668/1‑151 (MQ=255)
                cGTCTATGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGat                                                                                                                                  <  2:884298/151‑1 (MQ=255)
                 gTCGATGTGGGCACTCACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGaca                                                                                                                                 <  2:3153293/151‑3 (MQ=255)
                     aTGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGATGAGCCTGGCACCCCAGGCAACGAAGGCGACatcgt                                                                                                                             >  2:4642698/1‑149 (MQ=255)
                     aTGTCGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATatcat                                                                                                                             >  1:4229742/1‑151 (MQ=25)
                      cgtcggtagcaacgacccggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGaaggcgatctggtc                                                                                                                            >  1:1802863‑M2/106‑145 (MQ=255)
                           gCACTAACGCGCCGGTGGTCAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGt                                                                                                                       >  1:1000626/1‑151 (MQ=37)
                            tagcaacgacccggttgtgatgtcgatcgcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTc                                                                                                                      >  1:2089076‑M2/100‑151 (MQ=255)
                             aCTAACCAGCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCa                                                                                                                     >  1:4560919/1‑151 (MQ=255)
                              ccaacgcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCAc                                                                                                                    >  1:5240375/8‑151 (MQ=255)
                               caatgcgCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCAcc                                                                                                                   >  1:1542756/7‑151 (MQ=37)
                                 atgcgCCGGTGGTCAATGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCAt                                                                                                                 <  2:1642605/147‑1 (MQ=25)
                                   gcgCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCa                                                                                                               <  1:5334958/149‑1 (MQ=38)
                                     gCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCActattact                                                                                                             <  2:4181347/151‑8 (MQ=255)
                                       cGGTGGGCAGTGGCATGGTCGCGGTCAGCGGCGCGGATGCGTGGGAGTTGGTAAACCTGGACCTGGACAAGACCCCGGTCAGCACCGGGGTGACCGGCTAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGc                                                                                                           <  2:5109567/151‑1 (MQ=255)
                                           ggtCAATGCCATCGACGCGGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGAc                                                                                                       <  2:4064949/151‑1 (MQ=37)
                                            gtCAATGCCATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGAcc                                                                                                      <  1:3283389/151‑1 (MQ=38)
                                             tCAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGCCGACATCGTCAAAGTCACTATTACTGCTGACCa                                                                                                     >  2:2519643/1‑151 (MQ=37)
                                              cAATGCCATCGATGCCGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAg                                                                                                    >  1:2309395/1‑151 (MQ=25)
                                                tgtcgatcgctaccgttgaaggtgcggacgttggcaagttcgagcaactgacactcgacaaaaccccggtcagcactaccGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGAc                                                                                                  <  2:348327‑M2/71‑1 (MQ=255)
                                                aTGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGAc                                                                                                  >  1:4432302/1‑151 (MQ=25)
                                                aTGCCATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGAc                                                                                                  >  2:2048803/1‑151 (MQ=21)
                                                 tGCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTTAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGAcc                                                                                                 <  2:3248623/151‑1 (MQ=37)
                                                  gCCATCGACGCCGTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGAcgt                                                                                                <  1:1903166/151‑3 (MQ=25)
                                                   ccATCGGCGCGGTCAGCGGCGCAGATGCGTGTAAGTTCGAAAACCTGCACCTGGACAAGACCCCGGTCAGCACCGGGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTc                                                                                               <  2:2317214/151‑1 (MQ=25)
                                                   ccATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGCCAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGActtt                                                                                               >  2:1544252/1‑148 (MQ=25)
                                                   ccATCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTc                                                                                               <  2:3701129/151‑1 (MQ=21)
                                                   ccAGCGACGCCGTTAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTc                                                                                               <  2:4992671/151‑1 (MQ=21)
                                                    cATCGACGCGGTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACAGCCGACCAGACTTCg                                                                                              <  2:1747534/151‑1 (MQ=25)
                                                      tCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGt                                                                                            >  1:3881493/1‑151 (MQ=21)
                                                        gcaaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTgg                                                                                          <  2:1838251‑M2/80‑1 (MQ=255)
                                                          cgcgGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGcc                                                                                        >  2:2840701/1‑150 (MQ=21)
                                                          aaccgttgaaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGgcc                                                                                        <  1:1729340‑M2/82‑2 (MQ=255)
                                                           gcgGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGtagccg                                                                                       >  1:1661529/1‑146 (MQ=21)
                                                            cgGTCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGccga                                                                                      >  1:1340925/1‑148 (MQ=25)
                                                            ccgTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCga                                                                                      >  1:3080890/3‑151 (MQ=25)
                                                             ggTCAGCGGCGCAGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCgag                                                                                     >  2:4271909/1‑151 (MQ=25)
                                                             cgTCAGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCgag                                                                                     <  2:3905574/150‑1 (MQ=21)
                                                              gttaGCGGCGCTGATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGCCACCAGCACGGCCGACCACCCCTCGGTGCCCgaga                                                                                    >  2:4805964/4‑151 (MQ=18)
                                                              gTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCgaga                                                                                    >  1:199651/1‑151 (MQ=18)
                                                               tCAGCGGCGCTGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGaa                                                                                   >  1:4905632/1‑151 (MQ=25)
                                                                 aGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACg                                                                                 >  2:5406305/1‑151 (MQ=18)
                                                                  aaggcgcggatgtcggcaagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGt                                                                                >  1:2984831‑M2/62‑151 (MQ=255)
                                                                      cgcAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCaaa                                                                            >  1:2491113/1‑150 (MQ=17)
                                                                       gctgATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCaaac                                                                           >  2:1952358/4‑149 (MQ=21)
                                                                       gcggacgttggcaagttcgagcaactgacactcgacaaaaccccggtcagcactaccGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGtcaaac                                                                           <  1:3953290‑M2/94‑3 (MQ=255)
                                                                        ctgATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGcc                                                                          <  1:2197204/149‑1 (MQ=255)
                                                                         tgATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAAccg                                                                         >  2:1317951/2‑150 (MQ=18)
                                                                          gATGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAccga                                                                        <  1:5066346/151‑3 (MQ=17)
                                                                          gATGCGTGGAAGTTCGAAAACCTGAACCGGGACAAGACCCCGGGCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAAccga                                                                        <  1:3377337/151‑3 (MQ=17)
                                                                            tGCGTGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGAcc                                                                      >  2:3807338/1‑151 (MQ=17)
                                                                             gcgtGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCt                                                                     >  1:4830799/4‑151 (MQ=25)
                                                                              cATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCTGGCACCCCAGGCAACGAAGGCGACATCGTCAAAGTCACTATTACTGCTGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCtt                                                                    >  1:3779792/1‑151 (MQ=9)
                                                                               gtGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTc                                                                   <  1:1216431/150‑1 (MQ=21)
                                                                                tGGAAGTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCAACCTTCa                                                                  >  1:4976894/1‑151 (MQ=14)
                                                                                  gAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCAcc                                                                <  2:3378003/151‑1 (MQ=25)
                                                                                  caagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCAcc                                                                <  2:1448651‑M2/106‑1 (MQ=255)
                                                                                   aagttcgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCg                                                               >  1:61840‑M2/45‑151 (MQ=255)
                                                                                     gTTCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTg                                                             >  2:4885466/1‑151 (MQ=18)
                                                                                       tCGAAAACCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCa                                                           <  1:431625/151‑1 (MQ=18)
                                                                                        cgaacagctgacgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCAc                                                          >  1:911918‑M2/40‑151 (MQ=255)
                                                                                         gAAAACCTGAACCTGGACAAAACCCCCGTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACg                                                         <  2:699229/151‑1 (MQ=18)
                                                                                           aaaCCTGAACCTGGACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCCAACCGACCTTCACCGTGCACGTc                                                       >  2:2727947/1‑151 (MQ=17)
                                                                                            aacctgaacctgGACAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCa                                                      <  1:2677671/151‑1 (MQ=18)
                                                                                            aacctgaacctgCACAAGACCCCCGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCa                                                      <  1:4033683/151‑1 (MQ=14)
                                                                                                tgaacctgGACAAGCCTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCa                                                  <  1:2407984/151‑1 (MQ=14)
                                                                                                  acgctggacaagacgccggtcagcacgtcGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGACGTTCACCGTGCACAtcaacaccg                                                >  1:3561963‑M2/30‑147 (MQ=255)
                                                                                                      gggACAAGACCCCGGTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCt                                            <  1:5144449/150‑1 (MQ=14)
                                                                                                         aCAAGACTCCAGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGc                                         <  1:5422756/151‑1 (MQ=18)
                                                                                                          cAAGACTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGcc                                        <  2:4519234/151‑1 (MQ=14)
                                                                                                            aGACCCCGGTCAGCACCGAGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCa                                      >  1:413770/1‑151 (MQ=11)
                                                                                                             gaccccggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCAc                                     <  2:360924/144‑1 (MQ=17)
                                                                                                               ccccggTCAGCACCGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                   <  2:590989/146‑1 (MQ=17)
                                                                                                               cTCCAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGa                                   <  1:4656646/151‑1 (MQ=14)
                                                                                                                  cAGTCAACACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                >  2:2615218/1‑151 (MQ=11)
                                                                                                                   ggTCAGCAACGAGGTGTCCGACCAGCCAGGTACGCCGGGCAACGAGGGCGATCTGGTCCAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCATGTCAACCCGCCGCTGGCCCACGACTTg                               <  2:3040243/150‑1 (MQ=11)
                                                                                                                    gTCAGCACCGAGGTAACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTga                              <  2:2238251/151‑2 (MQ=17)
                                                                                                                        acACCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGt                          <  2:2917746/150‑1 (MQ=11)
                                                                                                                         cactaccGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTg                         <  1:347100‑M2/144‑1 (MQ=255)
                                                                                                                          actaccGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACTTGGTCAAAGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                        >  1:4736014‑M2/7‑151 (MQ=255)
                                                                                                                          aCCGAGGTGACCGACGAGCCAGGCACTCCAGGCAGCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                        <  1:2584331/151‑1 (MQ=11)
                                                                                                                            cGAGGTGACCGACGAGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGAcc                      >  2:4716682/1‑151 (MQ=17)
                                                                                                                             gAGGTGACCGACGAGCGAGGTACGGCGCGCAACGAAGGCGATCTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGAccc                     <  2:3977486/151‑1 (MQ=11)
                                                                                                                              aGGTCACCGACGAGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCt                    >  2:1200848/1‑151 (MQ=11)
                                                                                                                                gtcactgacgaGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGa                  >  2:797580/7‑151 (MQ=17)
                                                                                                                                gtcaCCGACGAGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGa                  >  1:1105557/4‑151 (MQ=14)
                                                                                                                                     cgacgaGCCAGGCACCCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaac             >  2:241833/1‑151 (MQ=14)
                                                                                                                                      gacgaGCCAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGcaaca            >  1:2882113/1‑151 (MQ=11)
                                                                                                                                       acgaGCCAGGCACTCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaacaa           >  2:1020411/1‑151 (MQ=17)
                                                                                                                                         gaGCCAGGTACCCCGGGTAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACTTGGTCGTGACCCTGAGCAACAACg         >  2:1090873/1‑151 (MQ=17)
                                                                                                                                         gaGCCAGGTACCCCAGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaacaatg         >  2:436712/1‑149 (MQ=11)
                                                                                                                                            ccAGGCACCCCGGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCATACGTCGGTGGCCGAGAACGTCAAACCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTCACCCTGGGCAACAACGccc      >  1:1434933/1‑151 (MQ=11)
                                                                                                                                                    |                                                                                                                                                  
CAACCGCCCCGGACAACGTCTATGTGGGCACTAACCAGCCGGTGGTCAATGCCATCGACGCGGTCAGCGGCGCAGATGCATGGAAGTTCGAAAACCTGAACCTGGACAAGACTCCAGTCAGCACCGAGGTGACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGATATCATCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCTGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGT  >  CP009974/1897446‑1897740

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: