Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,029,311 G→A 5.4% D707N (GAC→AAC)  htrA → S1C family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,029,3110GA 5.4% 250.6 / 3.1 85D707N (GAC→AAC) htrAS1C family peptidase
Reads supporting (aligned to +/- strand):  ref base G (30/50);  new base A (2/3);  total (32/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

ATTTAGAAACAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAATTACAGCG  >  CP000730/1029174‑1029451
                                                                                                                                         |                                                                                                                                            
gtttAGAAACAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGaca                                                                                                                                           <  2:989917/140‑1 (MQ=255)
aTTTAGAAACAAAAAGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGGCAAGCTGTTAAATTGCCAGGAAAAGTTAAGAACGGTGGTGTTGTAGATCAAGTTGaca                                                                                                                                           <  1:674214/141‑1 (MQ=255)
   tAGAAACAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGacaaca                                                                                                                                        >  2:917555/1‑141 (MQ=255)
       aaCAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATgg                                                                                                                                    >  1:866415/1‑141 (MQ=255)
        aCAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGt                                                                                                                                   >  1:563649/1‑141 (MQ=255)
        aCAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTAACAACAATGGt                                                                                                                                   >  2:964285/1‑141 (MQ=255)
        aCAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTAACAACAATGGt                                                                                                                                   >  2:341703/1‑141 (MQ=255)
          aaaaGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGttt                                                                                                                                 >  1:399753/1‑141 (MQ=255)
               gTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAg                                                                                                                            >  2:943479/1‑141 (MQ=255)
                  aaaTTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGAta                                                                                                                         >  2:821052/1‑140 (MQ=255)
                         cTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTCGCAGATCAATCTgg                                                                                                                  <  1:576510/141‑1 (MQ=255)
                          tATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGt                                                                                                                 <  1:682374/141‑1 (MQ=255)
                             cccGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGACCGGTGTTGTCGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTg                                                                                                              <  1:861133/141‑1 (MQ=255)
                             cccGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTg                                                                                                              <  1:730526/141‑1 (MQ=255)
                             cccGATGGAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTg                                                                                                              <  1:345418/141‑1 (MQ=255)
                                gATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGGTTAGCAGATCAATCTGGTCTGaaa                                                                                                           >  1:408781/1‑141 (MQ=255)
                                 aTGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGaaaa                                                                                                          <  1:233457/141‑1 (MQ=255)
                                 aTGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGaaaa                                                                                                          <  1:780731/141‑1 (MQ=255)
                                 aTGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTAGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGaaaa                                                                                                          <  2:674251/141‑1 (MQ=255)
                                 aTGTAGGTGTTAAAATGAAGAATATTGCCAGTATAAATAGTATAGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCATATCAAtaaaaaaaaaaaa                                                                                                          <  2:509492/141‑13 (MQ=255)
                                   gTAGGTTTTATAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGaaaaaa                                                                                                        <  1:44647/141‑1 (MQ=255)
                                   gTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGaaaaaa                                                                                                        <  1:380259/141‑1 (MQ=255)
                                         gtTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGAt                                                                                                  <  2:120273/141‑1 (MQ=255)
                                            aaaaTGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTa                                                                                               <  1:501215/141‑1 (MQ=255)
                                              aaTGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAAt                                                                                             <  2:839589/141‑1 (MQ=255)
                                               aTGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAAtt                                                                                            >  2:28056/1‑141 (MQ=255)
                                               aTGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAAtt                                                                                            >  1:654588/1‑141 (MQ=255)
                                               aTGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAAtt                                                                                            >  1:678687/1‑141 (MQ=255)
                                               aTGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAAtt                                                                                            >  2:769052/1‑141 (MQ=255)
                                                  aagaaTATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAACAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACt                                                                                         >  2:602061/1‑141 (MQ=255)
                                                   agaaTATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAGGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCGGAAAAAAAGTGAGGTAATTAATg                                                                                        >  2:601582/1‑141 (MQ=255)
                                                     aaTATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGaa                                                                                      <  1:841789/141‑1 (MQ=255)
                                                       tatTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAAtt                                                                                    <  2:1027343/141‑1 (MQ=255)
                                                         ttGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAg                                                                                  <  2:1080661/141‑1 (MQ=255)
                                                           ggcAGTCTGAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGAt                                                                                <  2:252300/139‑1 (MQ=255)
                                                           gcgaGTATAAATAGTTTTAAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGAt                                                                                <  1:75131/138‑1 (MQ=255)
                                                            tcAGTCTAAATAGGTTTGAAAGACAAGAAGTTAAATTGCGAGGAAAAGTTAAGAACGGTGCTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATAAAAAAAAAAAAAAAAAAGGTGATGTAATTACTGAATTAGATg                                                                               <  2:781265/140‑1 (MQ=255)
                                                            ccAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGCTTAGCAGATCAATCTGGCCGGAAAAAAGGTGATGTAATTACTGAATTAGATg                                                                               <  2:63739/141‑1 (MQ=255)
                                                            ccAGTCAAAATAGTTTTGAAAGACAAACAATTAAAGTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGAACAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATg                                                                               <  2:148843/141‑1 (MQ=255)
                                                             cAGTCTAAAAAGTATTAAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGAAAATCAAGTTGACAACAATGGTTTAGCAGATCAAAAAAAAAAAAAAAAAGGTGATGTAATTACTGAATTAGATgg                                                                              <  2:988239/141‑1 (MQ=255)
                                                              aGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGc                                                                             <  1:917555/141‑1 (MQ=255)
                                                                        gTTTTCAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTag                                                                   <  1:359831/141‑1 (MQ=255)
                                                                         tttgGAAAGACAAGCAGTGAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTaga                                                                  <  1:244230/141‑1 (MQ=255)
                                                                              aaaGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAgatg                                                             >  1:107484/1‑141 (MQ=255)
                                                                                aGACAAGCAGTTAAATTGCCAGGATAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAgatgat                                                           <  2:692183/141‑1 (MQ=255)
                                                                                aGACAAGCAGTTAAATTGCCAGGAAAGGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATGTGGTCTGAAAAAACGTGATGTAATTACTGAATTAGATCGCAAACTTTTAGAAgatgat                                                           <  1:673183/141‑1 (MQ=255)
                                                                                aGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCATATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAgatgat                                                           <  2:392553/141‑1 (MQ=255)
                                                                                aGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAgatgat                                                           <  1:801760/141‑1 (MQ=255)
                                                                                aGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAgatgat                                                           <  1:842827/141‑1 (MQ=255)
                                                                                aGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGGGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAgatgat                                                           >  2:108217/1‑141 (MQ=255)
                                                                                  gctgtgtggtcgaTTGCCGTGAAAAGTTAAGGACGGTGTTGTTGCACATCAAGTCGACATACATAGTTTATCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTAGTGAATTAGATGTCAAACTTTTAGAAGATGAttt                                                         <  2:327759/129‑1 (MQ=255)
                                                                                        cgTTAAATTGCCCGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCtt                                                   <  2:1048158/140‑1 (MQ=255)
                                                                                        cgTTAAATTGCCAGGAAAAGTTAAGGACGGTTTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTCTTCGAAGCTCCTTTACGCtt                                                   <  2:165024/140‑1 (MQ=255)
                                                                                        aattAAAGTGCCAGGAAAAGTTAAGATCGGGTTTGTGGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATAAAAAAAAAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCtt                                                   <  2:197850/139‑1 (MQ=255)
                                                                                              aTTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCa                                             <  2:1049038/141‑1 (MQ=255)
                                                                                               ttGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAg                                            <  2:563649/141‑1 (MQ=255)
                                                                                               ttGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTAACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAg                                            <  1:964285/141‑1 (MQ=255)
                                                                                               ttGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTAACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAg                                            <  1:341703/141‑1 (MQ=255)
                                                                                                   cAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATta                                        >  2:648204/1‑141 (MQ=255)
                                                                                                   cAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGGGATGTAATTACTGAATTAGATGGCAAACTCTTAGAAGATGATTTACGATTTAGGCCGATta                                        >  2:866554/1‑141 (MQ=255)
                                                                                                   cAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATAAAAAAAAAAAAAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATta                                        <  2:1054039/141‑1 (MQ=255)
                                                                                                     ggAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATtata                                      <  1:437606/141‑1 (MQ=255)
                                                                                                         aaGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTa                                  <  1:937181/141‑1 (MQ=255)
                                                                                                          aGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAg                                 >  2:355799/1‑141 (MQ=255)
                                                                                                          aGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAg                                 >  2:683316/1‑141 (MQ=255)
                                                                                                          aGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAg                                 >  2:985179/1‑141 (MQ=255)
                                                                                                              aaGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCAt                             >  1:356071/1‑141 (MQ=255)
                                                                                                               aGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATa                            <  1:1057079/141‑1 (MQ=255)
                                                                                                                   cGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGa                        <  1:472762/141‑1 (MQ=255)
                                                                                                                       gttgttgtAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGAc                    >  1:748887/1‑141 (MQ=255)
                                                                                                                        ttgttgtAAATCAAGTTGAAACAAATGGTTTAGCAGCTCAATCTGGTCTGAAAAAAGGTCATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACt                   <  1:182612/141‑1 (MQ=255)
                                                                                                                          gttgtAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTg                 >  2:336097/1‑141 (MQ=255)
                                                                                                                          gttgtAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTg                 >  2:593248/1‑141 (MQ=255)
                                                                                                                            tgtAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGaa               >  2:987032/1‑141 (MQ=255)
                                                                                                                            tgtAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGaa               >  2:579315/1‑141 (MQ=255)
                                                                                                                                gATTAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCa           >  2:1015799/1‑141 (MQ=255)
                                                                                                                                gATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCa           >  2:811976/1‑141 (MQ=255)
                                                                                                                                  tcgaGTTGACAACAATGGTTGAGCAGATCAATCTGGTCTGAAAAAAGGTGATTTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAAt         <  2:200916/138‑1 (MQ=255)
                                                                                                                                  tCAAGGTAACAACAATGGTTTAGCAGATTAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAAt         <  2:283121/141‑1 (MQ=255)
                                                                                                                                    aaGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAATTa       <  1:821052/141‑1 (MQ=255)
                                                                                                                                    aaGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATAATTAGTCATAAAGATGACTTGAAATCAATTa       >  2:700205/1‑141 (MQ=255)
                                                                                                                                      gctGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAACAGGTGATGTAATTATAGAATTTGATAAATAAGAATGGCGAGATTATGTAGGCTTTAGGCAGATGAGATTTAGTCATAAAGATGACTTGAAATCAATTACa     >  1:218262/3‑141 (MQ=255)
                                                                                                                                      gTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAATTACa     <  2:371875/141‑1 (MQ=255)
                                                                                                                                      gTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAATTACa     <  2:489878/141‑1 (MQ=255)
                                                                                                                                       ttGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAATTACAg    <  2:735879/141‑1 (MQ=255)
                                                                                                                                         gACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTCGAAATCAATTACAGCg  <  2:399753/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
ATTTAGAAACAAAAGGTAAAATTGACTATCCCGATGTAGGTGTTAAAATGAAGAATATTGCCAGTCTAAATAGTTTTGAAAGACAAGCAGTTAAATTGCCAGGAAAAGTTAAGAACGGTGTTGTTGTAGATCAAGTTGACAACAATGGTTTAGCAGATCAATCTGGTCTGAAAAAAGGTGATGTAATTACTGAATTAGATGGCAAACTTTTAGAAGATGATTTACGCTTTAGGCAGATTATATTTAGTCATAAAGATGACTTGAAATCAATTACAGCG  >  CP000730/1029174‑1029451

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