Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,367,702 A→C 8.8% I209S (ATC→AGC)  aspA ← aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,367,7020AC 8.8% 157.1 / 4.8 65I209S (ATC→AGC) aspAaspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (24/34);  new base C (3/4);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCA  >  NC_000913/4367563‑4367840
                                                                                                                                           |                                                                                                                                          
ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGc                                                                                                                                         >  1:495129/1‑143 (MQ=255)
ttCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGc                                                                                                                                         >  1:830214/1‑143 (MQ=255)
    gCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaa                                                                                                                                     <  2:655133/143‑1 (MQ=255)
     ccAGTTTGTTCACTGCGAGCGGAGAGTACTCTTTCGGCTTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                                                                                                    <  2:14527/143‑1 (MQ=255)
     ccAGTTCTTTCACTGCCAGCGGAGAGTACGCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGaaa                                                                                                                                    <  1:222957/143‑1 (MQ=255)
     ccAGTGTTTTCAGGGCCGGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTCACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGaaa                                                                                                                                    <  2:389035/143‑1 (MQ=255)
      cAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                                                                                                   <  1:42167/143‑1 (MQ=255)
      cAGTTGTGTCACTGCCAGCGGAGAGTACGCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAg                                                                                                                                   <  2:620251/143‑1 (MQ=255)
         ttttttCACTGCCAGCGGCGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCATCCTGCTGAAAgaag                                                                                                                                <  1:56538/143‑4 (MQ=255)
            tttCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGCTGCTGAAAGCGCGGa                                                                                                                             >  1:639336/1‑143 (MQ=255)
             ctcgcTGCCAGCCGAGCGTACGCTTTCGGCGTGTTCAGGCCAGTACCGATTGCTGGTGCACCAGGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAGGCGCGGaa                                                                                                                            <  2:859613/139‑1 (MQ=255)
                aCTGCCAGCTGTGAGTACTCTTTCGGCGTGCTCAGACCAGTTCCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTc                                                                                                                         <  2:483609/143‑1 (MQ=255)
                aCTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCTGTACGTTGTATGTTGTTGGTTTCTGCCTTCTTCATGATGCTGAAAGCGCGGAATTc                                                                                                                         <  2:66553/143‑1 (MQ=255)
                   gCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTGAGCAGGATGCTGAAAGCGCGGAATTCCTg                                                                                                                      >  1:85128/1‑143 (MQ=255)
                    cgaGCGGAGAGTGCGCTTGCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGa                                                                                                                     <  1:116737/141‑1 (MQ=255)
                     cAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGAc                                                                                                                    >  1:432184/1‑143 (MQ=255)
                       gCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCCCGGTTAACTTCCAGCAGCCGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCg                                                                                                                  >  1:545432/1‑143 (MQ=255)
                               tACTCTTTCGGAGTGTTCATACCAGTACCGAGTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCAt                                                                                                          <  2:140888/143‑1 (MQ=255)
                                     gtCGGCGTGGTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTGCGATGACATTCATTGCGTCTTTCAGCAGGAGGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTAc                                                                                                    <  1:209451/142‑1 (MQ=255)
                                        ggCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTAAAAGCGCGGAATTCCTGACCGAGGGTCATCGATACtac                                                                                                 >  1:306527/1‑141 (MQ=255)
                                           gtgtTCAGCCCAGTACCGGTTGCGGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTc                                                                                              <  2:774561/143‑1 (MQ=255)
                                           gtgtGCAGACCAGTACCGATTGCTGTTGCACGAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTGTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTc                                                                                              <  2:631619/143‑1 (MQ=255)
                                              ttCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGCTGCTGAAAGCGCGGAATTCTTGACCGAGGGTCATCAGTACTGCGTCCTg                                                                                           >  1:524998/1‑143 (MQ=255)
                                              ttCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGCCCGAGGGTCATCGGTACTGCGTCCTg                                                                                           >  1:688090/1‑143 (MQ=255)
                                              tgcAGACCAGTACCGATTGCTGTTGCACCAAGGTGAGCTTCCAGCAGCAGTTCAGCGGTACGTGGGATGTGTTTCACTGCTTCTTTCAGCGGGATGCTGAAGGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTg                                                                                           <  2:212500/141‑1 (MQ=255)
                                               tCAGACCAGTGCCGGTTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGc                                                                                          <  1:396175/143‑1 (MQ=255)
                                               tCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGc                                                                                          <  1:687449/143‑1 (MQ=255)
                                                     cAGTACCTTTTGCTGTTTCCCCAAGTTTTACTTCCAGCAGCAGTTCAGCTGTACGTTGGATGTTTTTCTCTTTTTCTTGCAGCAGGATGCTGAAAGCGCTGATTTCCTGACCGAGGGTTATCTGTACTGCGTCCTGCAGCTgg                                                                                    <  2:73985/143‑1 (MQ=255)
                                                     cAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTgg                                                                                    <  1:429298/143‑1 (MQ=255)
                                                     cAGTACCGAGTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTgg                                                                                    <  1:886617/143‑1 (MQ=255)
                                                     cAGGACCGATTGCTGTTGCACCAAGGTTAACTTCCATCATCAGTTCAGCGGTACGTTTGTTGTTTTTCACTTCTTATCCCGTGTGGATGCTGAAAGCGCGGAATTCCTGTCCGAGGGTCATCGGTACTGCGTCCTGCAGCTgg                                                                                    <  2:599456/143‑1 (MQ=255)
                                                        tACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGTGTa                                                                                 >  1:321862/1‑143 (MQ=255)
                                                        caCCGAGTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTa                                                                                 <  1:656176/142‑1 (MQ=255)
                                                         aCCGATGGCTGCCGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTGTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAc                                                                                <  1:229985/143‑1 (MQ=255)
                                                           cGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCCGGCTGCTGAACGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGa                                                                              >  1:98235/1‑143 (MQ=255)
                                                                      gCACGAAGGTTAACTTCCGGCGGCAGTTCGGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCGGAAAGCGCGGATTTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAg                                                                   <  2:210181/143‑1 (MQ=255)
                                                                          cAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTAGGACCCATTTTCAGGaag                                                               >  1:410743/1‑141 (MQ=255)
                                                                            aGGTTAACTTCCAGCGGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTc                                                             >  2:771801/1‑143 (MQ=255)
                                                                                tAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCCGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCCGCTGGGTACGAGCCCGTTTCAGCATGCCCggg                                                         >  2:731196/1‑143 (MQ=255)
                                                                                gtacctCCAGCAGCCGGTCAGCGGTATGTTGGATGCTTTTCGCTTCTTCTTTCAGCAGGAGGCTGAACGCGCGGAATTCCTGACGGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGGACGACCCATTTTCAGGATGTCCTgg                                                         <  1:361280/138‑1 (MQ=255)
                                                                                 aaCTTCCGGCAGCAGTTCAGCGGGACGTTGGATGTTTTTCACTTCTTGTGTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGa                                                        <  1:654145/143‑1 (MQ=255)
                                                                                 aaCTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGa                                                        <  2:875710/143‑1 (MQ=255)
                                                                                    ttCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAAtt                                                     >  1:184228/1‑143 (MQ=255)
                                                                                    ggcgagcagcagGTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAAtt                                                     <  1:493688/139‑1 (MQ=255)
                                                                                       cagcagcagTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGa                                                  <  1:892523/143‑1 (MQ=255)
                                                                                         gcagtcgTTCAGCGGTACGTTGGATGTTTTTCGCTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACa                                                <  1:897137/137‑1 (MQ=255)
                                                                                            gcagTTCAGCGGTACGTTGGATGCTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCt                                             <  1:379640/143‑1 (MQ=255)
                                                                                                ttCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATGCGACAGCTTTAc                                         >  1:179631/1‑143 (MQ=255)
                                                                                                   aGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGATTAATGCTCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCCTTTTCCGGCTGTCCTGGAATTCGCCCGCTTTCCGtg                                      >  1:468704/1‑142 (MQ=255)
                                                                                                      ggTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGCCCTGGAATTCGACAGCTTTACGTTCac                                   >  1:393141/1‑142 (MQ=255)
                                                                                                       gTACGTTGGATGTTTTTCACTTCTTCTTTCAGATGGATGCTGAAAGCGCTGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCGGGATGTCCTGGAATTCGACAGCTTTACGTTCaaa                                  <  1:909744/143‑1 (MQ=255)
                                                                                                        tACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCCGGAGGTCCTGGAAGTCGACAGCGTGACGGTCAAAg                                 >  2:470351/1‑143 (MQ=255)
                                                                                                            ttGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCtt                             >  2:506323/1‑143 (MQ=255)
                                                                                                            ttGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCtt                             >  2:497460/1‑143 (MQ=255)
                                                                                                                 tGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGCCCCATTTTCAGGATGTCCTGGAATTCGACAGATTTACGTACAAAGCCTTCACGc                        >  1:829416/1‑143 (MQ=255)
                                                                                                                    ttttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGt                     <  2:830214/143‑1 (MQ=255)
                                                                                                                    tgttCACTTCTTCTTTCAGCATGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGt                     <  2:495129/141‑1 (MQ=255)
                                                                                                                      ttCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATATTTTGGAAGTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTg                   >  1:513046/1‑143 (MQ=255)
                                                                                                                         aCTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGtt                >  2:6294/1‑143 (MQ=255)
                                                                                                                         aCTTCTTCCTGCAGCAGGATGCTGAAGGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGtt                <  1:394174/143‑1 (MQ=255)
                                                                                                                           ttcttcttTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTCCACGCCGTTGGTTaa              >  2:576452/1‑143 (MQ=255)
                                                                                                                                      gcagGATGCTGAAAGCGCGGAATTCCTGCCCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTAcc   <  1:312837/143‑1 (MQ=255)
                                                                                                                                       cagGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTAGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCa  >  1:178621/1‑143 (MQ=255)
                                                                                                                                          gATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCa                       >  2:168136/1‑119 (MQ=255)
                                                                                                                                          gATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCa                       <  1:168136/119‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TTCAGCCAGTTTTTTCACTGCCAGCGGAGAGTACTCTTTCGGCGTGTTCAGACCAGTACCGATTGCTGTTGCACCAAGGTTAACTTCCAGCAGCAGTTCAGCGGTACGTTGGATGTTTTTCACTTCTTCTTTCAGCAGGATGCTGAAAGCGCGGAATTCCTGACCGAGGGTCATCGGTACTGCGTCCTGCAGCTGGGTACGACCCATTTTCAGGATGTCCTGGAATTCGACAGCTTTACGTTCAAAGCCTTCACGCAGTTGGTTAATCGCATCTACCA  >  NC_000913/4367563‑4367840

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