Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,333,059 G→T 10.5% Q144K (CAA→AAA)  bsaA ← glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,333,0590GT10.5% 171.3 / 3.2 63Q144K (CAA→AAA) bsaAglutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base T (5/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TGTAATAAAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTCCCCTTTTACGT  >  CP000730/1332925‑1333199
                                                                                                                                      |                                                                                                                                            
tGTAATAAAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGc                                                                                                                                        <  1:164762/141‑1 (MQ=255)
  taataaAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCaa                                                                                                                                      <  2:696110/141‑1 (MQ=255)
  taataaAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCaa                                                                                                                                      >  1:521662/1‑141 (MQ=255)
   aataaAGACGCTTCATATAAAGGTCGAGTAAATATATGATACGCCTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTTAGGTGCaaa                                                                                                                                     <  1:369610/141‑1 (MQ=255)
   aataaAAACCCTTTATATTAAGGTCGAGTCAATATATGATCCGACTTTATATGAAGCGTTTTATTTTCTATCAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTTAGGTGCaaa                                                                                                                                     <  1:369537/141‑1 (MQ=255)
    ataaAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAAt                                                                                                                                    <  2:478662/141‑1 (MQ=255)
    ataaAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAAt                                                                                                                                    <  2:478549/141‑1 (MQ=255)
      aaaaaCGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCg                                                                                                                                  >  1:212270/1‑141 (MQ=255)
        aaaCGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGTAAATCGtt                                                                                                                                >  2:760361/1‑141 (MQ=255)
            gCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTaaa                                                                                                                            >  1:527540/1‑140 (MQ=255)
             cTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTaaca                                                                                                                           >  1:66708/1‑141 (MQ=255)
              ttCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTTCAAATCGTTTaacaa                                                                                                                          >  2:3357/1‑141 (MQ=255)
              ttCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTaacaa                                                                                                                          >  2:94594/1‑141 (MQ=255)
              ttCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTaacaa                                                                                                                          >  1:219360/1‑141 (MQ=255)
                 atatAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTTTTTTTCTTTTCGTTTAACAACGt                                                                                                                       >  2:344759/1‑141 (MQ=255)
                   gtaAAGGTCGAGTCGATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTAGAATAGTTTTTCTATTTCTCTTTCGATTTGACCAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTa                                                                                                                     <  2:691249/140‑1 (MQ=255)
                       aGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACACCGTTACCtt                                                                                                                 <  2:405932/141‑1 (MQ=255)
                       aGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCtt                                                                                                                 <  1:764234/141‑1 (MQ=255)
                        ggTTGAGTCTATATATGATTCAAGTTTATATGATGAGTTTTATTTGCTATGGGCTCGTATATTTTCTATTAGTGTTGCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTCTGGTGCAATTCGTTTACCATCGTTACCTTc                                                                                                                <  2:157057/141‑1 (MQ=255)
                           cGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTcgcg                                                                                                             >  2:289921/1‑141 (MQ=255)
                            gAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGa                                                                                                            >  1:325367/1‑141 (MQ=255)
                            gAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGa                                                                                                            >  2:225143/1‑141 (MQ=255)
                            gAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGa                                                                                                            >  2:710946/1‑141 (MQ=255)
                              gTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTTAGTTGCAAATCGTTTAACAACGTTACCTTCGCGATc                                                                                                          >  2:761268/1‑141 (MQ=255)
                               tCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATcc                                                                                                         <  1:494218/141‑1 (MQ=255)
                                 aaTATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCAc                                                                                                       >  1:455200/1‑141 (MQ=255)
                                            cGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAg                                                                                            >  2:479501/1‑141 (MQ=255)
                                              aCTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTg                                                                                          >  2:597191/1‑141 (MQ=255)
                                               cTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGa                                                                                         >  1:146651/1‑141 (MQ=255)
                                               cTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGa                                                                                         >  1:190242/1‑141 (MQ=255)
                                                   atatGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAAtt                                                                                     <  1:289921/141‑1 (MQ=255)
                                                             ttttATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGAtt                                                                           <  1:714262/141‑1 (MQ=255)
                                                             ttttATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGAtt                                                                           <  1:371693/141‑1 (MQ=255)
                                                              tttATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGAttt                                                                          >  2:145586/1‑141 (MQ=255)
                                                               ttATTTGCTATGAGTTAGTATATTTTATAATAATTTTTCTATTTCGCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGAtttt                                                                         <  1:546317/141‑1 (MQ=255)
                                                                 aTTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTtct                                                                       >  1:498042/1‑141 (MQ=255)
                                                                   ttGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTTCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCa                                                                     >  2:263144/1‑141 (MQ=255)
                                                                     gCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCAtt                                                                   <  2:522347/141‑1 (MQ=255)
                                                                      cTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTa                                                                  >  2:101254/1‑141 (MQ=255)
                                                                       tATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTaa                                                                 <  1:548032/141‑1 (MQ=255)
                                                                        aTGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTaaa                                                                <  1:731376/141‑1 (MQ=255)
                                                                           aGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGaa                                                             >  1:723738/1‑141 (MQ=255)
                                                                              tAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATcc                                                          >  1:450351/1‑141 (MQ=255)
                                                                              tAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATcc                                                          >  1:291756/1‑141 (MQ=255)
                                                                                 tataTTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCgtg                                                       >  1:736670/1‑141 (MQ=255)
                                                                                     ttttATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGtgttgt                                                   <  2:450351/141‑1 (MQ=255)
                                                                                         ataataATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCg                                               <  1:503900/141‑1 (MQ=255)
                                                                                           aataatTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTCGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCa                                             <  2:748990/141‑1 (MQ=255)
                                                                                              aatTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTTCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTc                                          >  2:192896/1‑141 (MQ=255)
                                                                                               atgtttCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTTTGTTGTTTATCGTTTATGAACGGGACCTTCGCGATCGACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTTTGGTATACGTGGTGTGcctctgttt                                         >  2:664535/4‑134 (MQ=255)
                                                                                                 ttttCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTTATTTTTATATCGTTTAACAACGTTACCTTCGCGATCTACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATTCGtgtttttctgctgttaaa                                       >  2:121892/1‑127 (MQ=255)
                                                                                                    tCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAatat                                    >  1:108776/1‑141 (MQ=255)
                                                                                                    tCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAatat                                    <  1:447083/141‑1 (MQ=255)
                                                                                                     cTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATc                                   <  2:190242/141‑1 (MQ=255)
                                                                                                             ctTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCtaaataaa                           <  1:94594/141‑1 (MQ=255)
                                                                                                              tttCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAg                          <  2:120185/141‑1 (MQ=255)
                                                                                                                    tttGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAAtt                    <  1:338185/141‑1 (MQ=255)
                                                                                                                     ttGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTg                   >  2:171017/1‑141 (MQ=255)
                                                                                                                     ttGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTg                   >  2:204532/1‑141 (MQ=255)
                                                                                                                           aGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTc             <  1:12015/141‑1 (MQ=255)
                                                                                                                           aGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTc             <  2:486330/141‑1 (MQ=255)
                                                                                                                             gTTTTTTTTTTTTTGTAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTccc           >  2:784137/1‑141 (MQ=255)
                                                                                                                                 tttttGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTCCCCttt       >  2:700315/1‑141 (MQ=255)
                                                                                                                                      gAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTCCCCTTTTACGt  <  2:584672/141‑1 (MQ=255)
                                                                                                                                      gAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTCCCCTTTTACGt  >  2:225347/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                            
TGTAATAAAAACGCTTCATATAAAGGTCGAGTCAATATATGATACGACTTTATATGAAGCGTTTTATTTGCTATGAGCTAGTATATTTTATAATAATTTTTCTATTTCTCTTTCGATTTGAACAGGTTTTTTTTGAGGTGCAAATCGTTTAACAACGTTACCTTCGCGATCCACTAAAAACTTAGTGAAATTCCATTTGATTTTCTCATTAAAGAATCCGTGTTGTGCCGCAGTCAAATATCTAAATAAAGGTAATTGATGTTCCCCTTTTACGT  >  CP000730/1332925‑1333199

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: ‑