Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 543,748 G→C 6.9% T100R (ACA→AGA)  allE ← S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913543,7480GC 6.9% 183.8 / 2.2 64T100R (ACA→AGA) allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (32/27);  new base C (3/2);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTG  >  NC_000913/543607‑543894
                                                                                                                                             |                                                                                                                                                  
aCCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACTAACGTCAGTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGcc                                                                                                                                             >  2:973410/1‑149 (MQ=255)
    gCCACGGGGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAAACATTTGCCttaa                                                                                                                                         <  1:766797/149‑3 (MQ=255)
        cGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGCCTGTCTTCGGCCTGGCCGTTAACAAAAGTCATTAAGTAGCCTGGCGGGCAATAAAGAAAGCCACCTTCGCTTAAGTCAAATGTTTTGCCTTCGGCtt                                                                                                                                     >  1:982972/1‑149 (MQ=255)
         ggTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCttt                                                                                                                                    >  2:84865/1‑149 (MQ=255)
          gTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTTTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATACCCACCTTCGCTTAAGGCAATTCTTTTGCCTTAGGCTTTg                                                                                                                                   >  2:109786/1‑149 (MQ=255)
           tGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTgg                                                                                                                                  <  2:96099/149‑1 (MQ=255)
             cATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCa                                                                                                                                <  1:267848/149‑1 (MQ=255)
               tAGCCCTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGt                                                                                                                              <  2:448385/149‑1 (MQ=255)
                aGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGCTAGCCCCCTTCGCTTAAGGCAAATGTTTAGCCTTCGGCTTTGGCAGTg                                                                                                                             >  2:348379/1‑149 (MQ=255)
                    ttCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAATCGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGatat                                                                                                                         >  2:939275/1‑149 (MQ=255)
                     tCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATAtt                                                                                                                        <  2:5607/149‑1 (MQ=255)
                       tACCGGAACATAGCGGCGCTTATATAAAAAGAGTTGGCGGTCTTCGGCCTGGGCGTTAACAACCGGCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTc                                                                                                                      <  1:585747/149‑1 (MQ=255)
                            gAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCagag                                                                                                                 >  1:50044/1‑149 (MQ=255)
                             aaCATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAgaga                                                                                                                <  1:891067/149‑1 (MQ=255)
                              acctAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCATATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGAt                                                                                                               <  2:773289/146‑1 (MQ=255)
                              aCATAGCGGCTCTTATATAAAAAGATTTGGATGTCTTCGGCCTGGGCGTTAACAAACGTCATTTAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGAt                                                                                                               <  2:715750/149‑1 (MQ=255)
                                aTAGCGGCGCTTATATAATAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATca                                                                                                             <  2:790976/149‑1 (MQ=255)
                                  aGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATcaca                                                                                                           >  1:298012/1‑149 (MQ=255)
                                   gcggcgCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACAt                                                                                                          >  2:415295/1‑149 (MQ=255)
                                     ggggcTTATATAAAAAGATTTGGCTGGCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATCAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCGTTGGCAGTGATATTTCCAGAGATCACATAc                                                                                                        <  2:807799/146‑1 (MQ=255)
                                      gcgcTTAAATAAAAACATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACa                                                                                                       <  2:835252/149‑1 (MQ=255)
                                       cgattATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAg                                                                                                      <  2:101340/146‑1 (MQ=255)
                                          ttATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGaa                                                                                                   <  1:6949/149‑1 (MQ=255)
                                           tatataAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCTGAGATCACATACAGGAAc                                                                                                  >  2:997703/1‑149 (MQ=255)
                                            atatatAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCAGTAAGGAGCCTGGCGGGCAATAAAGATAGCCACGTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACg                                                                                                 >  2:209811/1‑149 (MQ=255)
                                              ataAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGtt                                                                                               >  1:409955/1‑149 (MQ=255)
                                               taAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGCTAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGGTGTGGCAGTGATAATTCCAGAGATCACATACAGGAACGttt                                                                                              >  2:750099/1‑149 (MQ=255)
                                                aaaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGttta                                                                                             >  1:213759/1‑148 (MQ=255)
                                                aaaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTc                                                                                             >  1:76098/1‑149 (MQ=255)
                                                aaaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAACGCAAATGTTTTGCCTTCGGCTTTGGCACTGATATTTCCAGAGATCACATACAGGAACGTTTc                                                                                             >  2:244356/1‑149 (MQ=255)
                                                 aaaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCa                                                                                            >  2:520181/1‑149 (MQ=255)
                                                  aaaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAACCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCaa                                                                                           >  1:805870/1‑149 (MQ=255)
                                                   aaGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAAt                                                                                          <  1:534489/149‑1 (MQ=255)
                                                        ttGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCtt                                                                                     <  1:854151/149‑1 (MQ=255)
                                                                 ttCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAATTGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACAGTCAGGACCGTTTCAATGCCTTCGCCACCgt                                                                            >  1:619454/1‑148 (MQ=255)
                                                                 ttCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCCCCGa                                                                            >  1:781738/1‑149 (MQ=255)
                                                                 ttCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGa                                                                            >  1:979972/1‑149 (MQ=255)
                                                                 ttCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGa                                                                            >  1:940999/1‑149 (MQ=255)
                                                                  tCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCTAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGaa                                                                           >  2:482610/1‑149 (MQ=255)
                                                                           gCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCTAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTgttg                                                                  <  2:878831/149‑1 (MQ=255)
                                                                            cGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTgttgg                                                                 >  2:993972/1‑149 (MQ=255)
                                                                            cGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTCTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTgttgg                                                                 <  1:482489/149‑1 (MQ=255)
                                                                                 aCAAACGTCATTAAGGAGCCTGGAGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGACAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCTCCAGCGAAGCCCTGTTGGTTTcc                                                            >  1:923839/1‑149 (MQ=255)
                                                                                     aCGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTTATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCg                                                        >  1:495734/1‑149 (MQ=255)
                                                                                       gTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGGAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGtt                                                      >  1:934248/1‑149 (MQ=255)
                                                                                        tCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGTTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGttt                                                     <  1:183220/149‑1 (MQ=255)
                                                                                          aTTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTg                                                   >  1:16230/1‑149 (MQ=255)
                                                                                                taGCCTGGCGGGCAATAAAGATAGCCACCTTCGATTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAg                                             <  1:65615/148‑1 (MQ=255)
                                                                                                gAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAg                                             <  2:76098/149‑1 (MQ=255)
                                                                                                gAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACAGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAg                                             >  1:182741/1‑149 (MQ=255)
                                                                                                      ggCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGAc                                       <  1:520181/149‑1 (MQ=255)
                                                                                                       gCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCGTTGGCAGTGATATTTCCAGAGATCACATACAGGACCGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGAcc                                      <  1:49041/149‑1 (MQ=255)
                                                                                                          ggCAATAAAGATAGCCACCTTCTCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGa                                   <  1:948931/149‑1 (MQ=255)
                                                                                                              aTAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACTTACAGGAACGTTTCAATGCCTTCGCCCCCGAAGCCCTGTTGGTTGCCCCCGTGTTGATGCAGTGTGACCAGATAAt                               >  1:155083/1‑149 (MQ=255)
                                                                                                              aTAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAAt                               <  1:543965/149‑1 (MQ=255)
                                                                                                                      cgCCCCCTTCGCTACAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTTAGATTTCCAGAGATCACATACAGGACCGTTTCAATGCCTTCCCCCCCTAAGCCCTGTTGGTTGCCACCGTTTTCATGCAGTGTGACCAGATAATCAACaaaa                       <  2:490123/148‑1 (MQ=255)
                                                                                                                      agagaacTTGGCTTAAGGTAAATCTTTTGCGTTCGGCTTTGGCAGTGATAATTCCTTTTATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTTCGACCGTTTTGATGCAGTGTGACCAGATAATCAACaaaa                       <  2:982972/143‑1 (MQ=255)
                                                                                                                        ccACCTTCGCTTAAGGCAAATGTTTTGAATTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCACTGTTGGTTGCCAGCGTTTTGATGCAGTGTGACCAGATACTCACCAAAaga                     >  2:54755/1‑149 (MQ=255)
                                                                                                                         cACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGCGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAagag                    >  1:568413/1‑149 (MQ=255)
                                                                                                                         cACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAagag                    <  1:522680/149‑1 (MQ=255)
                                                                                                                          aCCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGCAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAgg                   >  2:136504/1‑149 (MQ=255)
                                                                                                                                 ccgaaGGCAAAACATTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGGCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGAGAATCAACTTAAGAGcccccccg            >  2:540577/4‑147 (MQ=255)
                                                                                                                                   gaaGGCAAAACATTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCt          >  1:490462/2‑149 (MQ=255)
                                                                                                                                           aTGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTg  <  1:84865/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
ACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAAAAGATTTGGCTGTCTTCGGCCTGGGCGTTAACAAACGTCATTAAGGAGCCTGGCGGGCAATAAAGATAGCCACCTTCGCTTAAGGCAAATGTTTTGCCTTCGGCTTTGGCAGTGATATTTCCAGAGATCACATACAGGAACGTTTCAATGCCTTCGCCACCGAAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTG  >  NC_000913/543607‑543894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑