Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,191,159 C→A 7.8% pseudogene (1207/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,191,1590CA 7.8% 900.8 / 34.7 243pseudogene (1207/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base C (104/120);  new base A (10/9);  total (114/129)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 243 total aligned reads displayed.

TAGAAAAATATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCAAAAAGAAGA  >  CP009273/2191013‑2191307
                                                                                                                                                  |                                                                                                                                                    
tAGAAAAATATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTc                                                                                                                                                   >  1:2632761/1‑150 (MQ=255)
  gAAAAATATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCat                                                                                                                                                 <  2:2190181/150‑1 (MQ=255)
  gAAAAATATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCat                                                                                                                                                 <  2:326020/150‑1 (MQ=255)
   aaaaaTATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCata                                                                                                                                                <  2:3695056/150‑1 (MQ=255)
    aaaaTATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCatat                                                                                                                                               >  2:3869519/1‑150 (MQ=255)
        tatCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGAGTCATatcaat                                                                                                                                           <  2:2787608/150‑1 (MQ=255)
                aCCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGtt                                                                                                                                   <  2:2756394/150‑1 (MQ=255)
                aCCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGtt                                                                                                                                   <  2:3741691/150‑1 (MQ=255)
                 ccATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTc                                                                                                                                  <  2:2394902/150‑1 (MQ=255)
                    tatGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCa                                                                                                                               <  1:2717958/150‑1 (MQ=255)
                     atGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCAt                                                                                                                              <  1:1028906/150‑1 (MQ=255)
                       gCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTg                                                                                                                            >  2:2134431/1‑150 (MQ=255)
                       gCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTg                                                                                                                            <  2:176080/150‑1 (MQ=255)
                       gCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTg                                                                                                                            <  2:176313/150‑1 (MQ=255)
                         cATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGac                                                                                                                          <  1:947697/150‑1 (MQ=255)
                          aTTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGaca                                                                                                                         <  1:605409/150‑1 (MQ=255)
                           tttttGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGacac                                                                                                                        <  1:1818246/150‑1 (MQ=255)
                           tttttGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGacac                                                                                                                        <  1:1180379/150‑1 (MQ=255)
                             tttGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTg                                                                                                                      <  2:1456519/150‑1 (MQ=255)
                             tttGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGAGTCATATCAATCATCCGTTCGCATTGACACTg                                                                                                                      <  2:673281/150‑1 (MQ=255)
                              ttGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGc                                                                                                                     >  2:2214481/1‑150 (MQ=255)
                               tgatgaTTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCt                                                                                                                    <  1:136034/150‑1 (MQ=255)
                                gatgatTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCtt                                                                                                                   >  1:2294321/1‑150 (MQ=255)
                                gatgatTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCtt                                                                                                                   <  2:502930/150‑1 (MQ=255)
                                gatgatTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCtt                                                                                                                   >  1:1226755/1‑150 (MQ=255)
                                gatgatTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCtt                                                                                                                   <  2:3007151/150‑1 (MQ=255)
                                gatgatTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGAGTCATATCAATCATCCGTTCGCATTGACACTGCtt                                                                                                                   <  2:138431/150‑1 (MQ=255)
                                 atgatTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTta                                                                                                                  <  2:3518432/150‑1 (MQ=255)
                                     ttACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACg                                                                                                              <  2:133087/150‑1 (MQ=255)
                                     ttACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACg                                                                                                              >  2:3790262/1‑150 (MQ=255)
                                          aTTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAg                                                                                                         >  2:31958/1‑150 (MQ=255)
                                          aTTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAg                                                                                                         >  2:31794/1‑150 (MQ=255)
                                          aTTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAg                                                                                                         >  1:3517157/1‑150 (MQ=255)
                                            tgtgGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGcc                                                                                                       <  2:13707/150‑1 (MQ=255)
                                            tgtgGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGcc                                                                                                       <  2:13576/150‑1 (MQ=255)
                                             gtggtgCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCg                                                                                                      <  1:3472434/150‑1 (MQ=255)
                                              tggtgCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCgg                                                                                                     >  1:3954105/1‑150 (MQ=255)
                                              tggtgCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCgg                                                                                                     <  2:1480525/150‑1 (MQ=255)
                                               ggtgCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCggg                                                                                                    <  2:368668/150‑1 (MQ=255)
                                                 tgCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCa                                                                                                  <  1:1165756/150‑1 (MQ=255)
                                                   cGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCaaa                                                                                                >  2:153937/1‑150 (MQ=255)
                                                   cGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCaaa                                                                                                <  1:605570/150‑1 (MQ=255)
                                                    gATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAAc                                                                                               <  1:3247384/150‑1 (MQ=255)
                                                      tatGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTa                                                                                             <  2:3281811/150‑1 (MQ=255)
                                                      tatGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTa                                                                                             <  2:3636343/150‑1 (MQ=255)
                                                         ggTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAAc                                                                                          >  1:805805/1‑150 (MQ=255)
                                                         ggTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAAc                                                                                          <  2:210854/150‑1 (MQ=255)
                                                          gTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACg                                                                                         >  2:3619554/1‑150 (MQ=255)
                                                          gTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACg                                                                                         <  1:1514996/150‑1 (MQ=255)
                                                          gTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACg                                                                                         <  1:801706/150‑1 (MQ=255)
                                                                 aCCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCAc                                                                                  <  2:1136590/150‑1 (MQ=255)
                                                                    gTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGAt                                                                               >  2:1970352/1‑150 (MQ=255)
                                                                       ttACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGAGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCgg                                                                            <  2:2678682/150‑1 (MQ=255)
                                                                        tACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGa                                                                           <  1:1619621/150‑1 (MQ=255)
                                                                        tACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGa                                                                           <  1:1620691/150‑1 (MQ=255)
                                                                          cAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATg                                                                         >  2:2492993/1‑150 (MQ=255)
                                                                             gAGCAAGGTGTTGCAGCCCAGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCACACTAAACGCTGTCACGATCGGATGACg                                                                      >  1:970316/1‑150 (MQ=255)
                                                                              aGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGa                                                                     >  1:1858269/1‑150 (MQ=255)
                                                                                 aaGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAg                                                                  >  1:2071701/1‑150 (MQ=255)
                                                                                  aGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGc                                                                 >  2:624583/1‑150 (MQ=255)
                                                                                      gtTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCAtt                                                             >  1:4293940/1‑150 (MQ=255)
                                                                                      gtTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCAtt                                                             >  1:1198871/1‑150 (MQ=255)
                                                                                      gtTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCAtt                                                             <  1:4350027/150‑1 (MQ=255)
                                                                                      gtTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCAtt                                                             <  2:2142573/150‑1 (MQ=255)
                                                                                      gtTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCAtt                                                             <  2:3155328/150‑1 (MQ=255)
                                                                                         gCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATtgct                                                          <  2:2542953/150‑1 (MQ=255)
                                                                                         gCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGAGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATtgct                                                          >  1:1161402/1‑150 (MQ=255)
                                                                                                     cGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACAt                                              <  1:984777/150‑1 (MQ=255)
                                                                                                     cGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACAt                                              <  2:3957926/150‑1 (MQ=255)
                                                                                                     cGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGAGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACAt                                              <  1:3859492/150‑1 (MQ=255)
                                                                                                       aTTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTGACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATgc                                            <  1:3692469/150‑1 (MQ=255)
                                                                                                       aTTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATgc                                            >  1:701668/1‑150 (MQ=255)
                                                                                                         ttGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATgcag                                          <  1:3679428/150‑1 (MQ=255)
                                                                                                         ttGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATgcag                                          >  1:4260655/1‑150 (MQ=255)
                                                                                                            cTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCaa                                       <  2:2467615/150‑1 (MQ=255)
                                                                                                              cccTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCTATg                                     >  1:2906195/1‑150 (MQ=255)
                                                                                                               ccTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATgg                                    >  2:4354247/1‑150 (MQ=255)
                                                                                                                cTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATggc                                   >  1:3739461/1‑150 (MQ=255)
                                                                                                                cTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATggc                                   >  1:3380161/1‑150 (MQ=255)
                                                                                                                 tATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATggcg                                  >  1:926843/1‑150 (MQ=255)
                                                                                                                      cGCAAGTGACTACTGCGCCGATGTTTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACt                             >  1:4448478/1‑150 (MQ=255)
                                                                                                                      cGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACt                             >  1:1389546/1‑150 (MQ=255)
                                                                                                                      cGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACt                             >  1:3658548/1‑150 (MQ=255)
                                                                                                                         aaGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGAc                          <  1:4068353/150‑1 (MQ=255)
                                                                                                                         aaGTGACTACTGCGCCGATGTGTTGAGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGAc                          >  2:1874817/1‑150 (MQ=255)
                                                                                                                           gTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACgg                        <  1:3162229/150‑1 (MQ=255)
                                                                                                                            tGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGa                       <  1:1009136/150‑1 (MQ=255)
                                                                                                                             gACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGaa                      >  1:4014716/1‑150 (MQ=255)
                                                                                                                             gACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGaa                      >  2:4079621/1‑150 (MQ=255)
                                                                                                                                tactGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAActt                   <  1:2134431/150‑1 (MQ=255)
                                                                                                                                     cgcCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTgg              <  1:2886441/150‑1 (MQ=255)
                                                                                                                                      gcCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTggg             <  2:4341569/150‑1 (MQ=255)
                                                                                                                                      gcCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTggg             >  1:2322626/1‑150 (MQ=255)
                                                                                                                                         gATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCaa          <  1:3867051/150‑1 (MQ=255)
                                                                                                                                          aTGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCaaa         >  2:3518200/1‑150 (MQ=255)
                                                                                                                                          aTGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCaaa         >  2:3518081/1‑150 (MQ=255)
                                                                                                                                          aTGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCaaa         >  2:2593672/1‑150 (MQ=255)
                                                                                                                                            gtgtTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCaaaaa       <  1:2243145/150‑1 (MQ=255)
                                                                                                                                             tgtTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCAAAaag      <  1:2873874/150‑1 (MQ=255)
                                                                                                                                               ttGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCAAAaagaa    >  2:4255748/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAGAAAAATATCACGAACCATATGCCATTTTTGATGATTACTATTGTGGTGCGATATGGTCAGCCACCGTATTACAGGAGCAAGGTGTTGCAGCCCTGCCCCGATTTGCTCCCTATGCCGCAAGTGACTACTGCGCCGATGTGTTGCGTCATATCAATCATCCGTTCGCATTGACACTGCTTATACGTGTAGCCGGGCAAACTAAACGCTGTCACGATCGGATGACGAAAGCCATTGCTGCGTTCCCACATGCAGCAATGGCGGCACTGACGGAACTTCTTGGGCAAAAAGAAGA  >  CP009273/2191013‑2191307

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