Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,323,426 T→C 15.9% K309K (AAA→AAG phnD ← phosphonate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,323,4260TC15.9% 135.4 / 3.0 69K309K (AAA→AAGphnDphosphonate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base C (11/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCAGCCCAGGCGTTCCAGCACC  >  NC_000913/4323285‑4323567
                                                                                                                                             |                                                                                                                                             
cGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGACTCTCGGTGGttt                                                                                                                                              >  2:600773/1‑143 (MQ=255)
  tcggcccaTTGGCGGCTGCGACGCCGGGGCGTCGTGGGCGGCGTGCGCGGTGCGGCACCGCTGTGCTCAGCGAACTCCTCGCGCTTAGCGCGTCGTTCAGGCGGTCCAGGTCGTGCAGCTGCGCCTGAATCTCGGTGGttttt                                                                                                                                            <  2:450089/137‑1 (MQ=255)
  ccggCACATTGCCGGAATCGACGCCGGGACGTCGTATGAGGCCTAGGCGTTACGGCACCGCCTTACTCATCGAGCTCACCGCGCTTAGCGCACCGATCAGGTGGTGCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGttttt                                                                                                                                            <  2:642262/140‑1 (MQ=255)
   cTGCACATTGTCGGATGCGACGCCGGTGCGTCTTATCCGGCCTACGATGTACTGCACCGCTTTGCTCACCGAGCTCCTCCCTCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTg                                                                                                                                           <  1:607933/143‑1 (MQ=255)
      cacaTTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCGCGGTGGTTTTTGCCa                                                                                                                                        >  2:498817/1‑143 (MQ=255)
             tCGGATGCGACGCCGGGGCGTCTTATCCGCCCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                                                                                                                                 <  1:468958/143‑1 (MQ=255)
             ccGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTc                                                                                                                                 <  1:303458/142‑1 (MQ=255)
                gATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTg                                                                                                                              <  2:24496/143‑1 (MQ=255)
                 aTGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGc                                                                                                                             >  1:55203/1‑143 (MQ=255)
                     gACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCAt                                                                                                                         >  2:107716/1‑143 (MQ=255)
                      aCGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCAtt                                                                                                                        <  1:564936/143‑1 (MQ=255)
                                   ttctCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTg                                                                                                           <  1:71084/140‑1 (MQ=255)
                                   gtATCCGACCTACGACTTACTGCACCGCTTTAGCCAGGTAACTCTTCGCGCTTATCTCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTg                                                                                                           <  2:486215/142‑1 (MQ=255)
                                      tCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCTTTGCCAGCTTGTCGTCCTCATTCAGTCGTTTCTTGCtt                                                                                                        >  2:170326/1‑143 (MQ=255)
                                      tCCGACCTAAGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCCTTCAGTCCTTTGTTGCtt                                                                                                        >  1:413496/1‑143 (MQ=255)
                                       ccGCCCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCttt                                                                                                       <  1:370002/143‑1 (MQ=255)
                                        cGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCtttt                                                                                                      >  2:332299/1‑143 (MQ=255)
                                         gACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCATTTATTGCCTTTc                                                                                                     >  2:217370/1‑143 (MQ=255)
                                             tACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGc                                                                                                 <  1:655095/143‑1 (MQ=255)
                                             tACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTCGCCACCTTAACCTGCTCATTCAGTCCTTTATTGCTTTTCgcgc                                                                                                 >  2:28223/1‑143 (MQ=255)
                                             tACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGCCGCAAGCTTATCATGCTCATTCAGTCCTTTCTAGCTTTTCACGc                                                                                                 >  2:63500/1‑143 (MQ=255)
                                              aCGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGcc                                                                                                <  1:170629/143‑1 (MQ=255)
                                                 aTTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTg                                                                                             <  1:175929/143‑1 (MQ=255)
                                                  tttACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTCTTGAATGATGCGTCGCCTCCTTCAATCCCTGGTTGCTTTTCACGCCCTGc                                                                                            >  1:490934/1‑143 (MQ=255)
                                                    ttccgCACCACTTTACACACCGTACTGAACGCGTGTAGATCGTTGTCCAGGCGGTCTAGGTCATCCAGCTGCTCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCAt                                                                                          <  2:193236/139‑1 (MQ=255)
                                                      cTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGCCGTCTGATTGTCCTGCACAGTCAGTCCTTTATTGCTTTCAACTCCCTGCctct                                                                                        >  2:103951/1‑143 (MQ=255)
                                                          aCCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCttt                                                                                    <  1:429860/143‑1 (MQ=255)
                                                          aCCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCttt                                                                                    <  1:173240/143‑1 (MQ=255)
                                                          aCCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTAACGCCCTGCATCTCttt                                                                                    >  1:109082/1‑143 (MQ=255)
                                                            cGCTTTACTCACCGAACTCATCGCGCTGAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTaa                                                                                  <  1:163007/143‑1 (MQ=255)
                                                                      aCCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCTGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTTAACAGCTCGAg                                                                        >  1:135182/1‑143 (MQ=255)
                                                                      aCCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAg                                                                        >  1:454030/1‑143 (MQ=255)
                                                                          aaCTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTgg                                                                    <  1:612326/143‑1 (MQ=255)
                                                                          aaCTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGGCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTgg                                                                    <  1:457990/143‑1 (MQ=255)
                                                                               aTCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAAt                                                               >  1:109737/1‑143 (MQ=255)
                                                                               aTCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAAt                                                               >  1:434745/1‑143 (MQ=255)
                                                                               aTCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAAt                                                               >  1:224576/1‑143 (MQ=255)
                                                                               aTCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCTTTGTCAGCTTGTCATGCTCATTCCGTCCTTTGTTGCGTTTCCCGCCCTGCATCACTTTAAAAAGCGCGAGCTgggcgaa                                                               >  2:338779/1‑138 (MQ=255)
                                                                                    gcTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCa                                                          <  1:196908/143‑1 (MQ=255)
                                                                                     cTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCAc                                                         <  1:72240/143‑1 (MQ=255)
                                                                                     cTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCAc                                                         <  1:258774/143‑1 (MQ=255)
                                                                                      ttAGCGCGTTGTTCAGGCGGTCCAGGTAATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTACAGATCATTCAGTGCTTTATTGCTTTTCTCGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCAAATCGGccac                                                        >  2:536711/1‑140 (MQ=255)
                                                                                          cgcgTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGCTGCCTCCTTGTCCAGCTAATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCCAGCGGTCGACTAGGCAGCAGtc                                                    >  2:24753/1‑142 (MQ=255)
                                                                                          cgcgTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAAGCGCGGTGGTTTTTGCCAGCTTGTCGTGCTCAATGAGTCATTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTAGCGGATCGGCACCAGtg                                                    >  1:482701/1‑142 (MQ=255)
                                                                                           tctttgtgCAGGCGGTGCAGGTCATCCAGCTGCGCCGGCATCTCGGTGGGTTTTGCCAGCTTGTCGAGATGATTCAGGCCTTTCTTGCTTCTCACGCCCTGCATCTCTTTGAGCAGCGCGAGCTGGCGAATCGGCACCAGTTg                                                   <  2:80890/140‑1 (MQ=255)
                                                                                           gcgTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCTTCGTCAGCTTGAACTACACATTCAATCCTCTCTTGCTTTTCACTCCCTGCATCTCTTCAACCACCCCCACCAGGCTACTCGccaccacatg                                                   >  2:521008/1‑139 (MQ=255)
                                                                                              ttgttCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAg                                                <  1:599750/143‑1 (MQ=255)
                                                                                              ttgttCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTGTCGCCAGCTTATACTGCTCATGCAGTCCTTTATTGCTGTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAAGCGGCGCCAGTTGCAg                                                >  2:150989/1‑143 (MQ=255)
                                                                                                  tCAGGCGGTCCAGGTCACCCAGCTGCCCCTTAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCGTTATTGCGTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCg                                            <  1:21101/143‑1 (MQ=255)
                                                                                                   cAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGc                                           >  1:589325/1‑143 (MQ=255)
                                                                                                    cggCGGTGCAGGTCATCCCGCCGCGCCGGAATCTCGGTGGTTGTTGCCGGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCGGGCGAATCGGCACCAGTTGCAGGTCGCt                                          <  2:166182/142‑1 (MQ=255)
                                                                                                     ggCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTg                                         <  1:416569/143‑1 (MQ=255)
                                                                                                         gTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGCTGCCACCTTCAACGACTCATTCAGTCATTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCGGGCGAATCCGCACCAGTACCAGGCCGCTGGAAg                                     >  2:533276/1‑143 (MQ=255)
                                                                                                          tCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCTTTGCCTGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGACTCGGCACCAGTTGCCACTCGATGGAAgc                                    >  2:358192/1‑143 (MQ=255)
                                                                                                          tCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGCCGCCACCTTGACCTGCACATTCAGTCCTTCATTGCTTTTCACGCCCGGCATCTCTCTAAACAGCGCGAGCTGGCAACTCGGCACCAGTTGCACGCCGCTgcacgc                                    >  2:191003/1‑138 (MQ=255)
                                                                                                             aGGTCCTGCAGCCGTGCCTGAATCTCGGTGGTTTTTGCCAGCTCGTCCTGCTCACTCGGTCCTTTCTTGCTTTTCACGCCCTGCATCTCTTTCCACAGCGCGCGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGcgcg                                 <  2:111609/143‑1 (MQ=255)
                                                                                                             aGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTGTTGCAAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAACTGGCGAATCGGCACCACTTGCAGTGCGCTGGAcgcgcg                                 >  2:541034/1‑143 (MQ=255)
                                                                                                             aGGTCATCCAGCTGCGCCTGAATCTCGGTGGTCGTTGCCAGCTTGTCCTGCACATTCAGTCCTGTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCCCCAGTTGCAGGTCGCTGGAcgcgcg                                 >  2:73413/1‑143 (MQ=255)
                                                                                                               gTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGa                               <  1:547003/143‑1 (MQ=255)
                                                                                                                tCATCCAGCTGCGCCTGAATCTCGGTGGTCGCCGCCACCTTGACCTGCACATTCCGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGaa                              >  2:57924/1‑143 (MQ=255)
                                                                                                                  aTCCAGCTGCGCCTGAATCTCGGTGGTTATTGTCACCTTGTCCACAACATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCCACCTAGCTAATTGGCACCAGTTGCAGGACACCGGAACCGCGGAAta                            >  2:373093/1‑142 (MQ=255)
                                                                                                                     cAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCCGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGgcg                         >  1:571454/1‑143 (MQ=255)
                                                                                                                      aGCTGCGCCTGAATCCCGGTGGTTCTTGCCAGCTTGTCCTGCTCATTCAGGCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGgcgc                        <  2:134080/143‑1 (MQ=255)
                                                                                                                        cTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgccc                      >  1:75078/1‑143 (MQ=255)
                                                                                                                         tggcccTGATTCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgccca                     <  2:577053/139‑1 (MQ=255)
                                                                                                                         tGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgccca                     <  1:498817/143‑1 (MQ=255)
                                                                                                                              cTGAATCTCGGTGGTTTTCGCAACCTTAACCAGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCgcacagcaca                >  2:541350/1‑143 (MQ=255)
                                                                                                                                   gctcGGTGGTCTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCAGCCCAGGCGt           <  1:68652/142‑1 (MQ=255)
                                                                                                                                     ttcgtggtCTTCCCCAGCTTGTCCTCCTCCTTCCTCCCTTGATTCCTTGTCACGCCCTGCCTCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCCGGTCTCTGGAAGCGCGGAATGGCGCCCAGCCCAGGCGTTc         <  2:85740/140‑1 (MQ=255)
                                                                                                                                     tcGGTGGTTGTTGCCTGCTTGAACTGCTCATTCAGTCCTTTATTGCTTTTCACGCCTTGCATCTCTTTAAACAGCGCGAGCTGGCGACTCGGCCCCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCCGCCCTGGCgcgc         >  2:351715/1‑140 (MQ=255)
                                                                                                                                        gtggtTTTCTCCAGCTTGGCCTGCTCATTCAGTCCTGTATTGCTCTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCAGCCCAGGCGTTCCAg      <  2:398167/143‑1 (MQ=255)
                                                                                                                                          ggtTCTTGCCAGCTTGTCCTGCCCCGGCAGGCCGTTATCGCCGTTCACGCCCTGCCTCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAAGGGCGCCCAGCCCAGGCGTTCCAGCa    <  2:177950/143‑1 (MQ=255)
                                                                                                                                            tttttGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAACGGCGCCCCGGCCAGCCGTTCCAGcccc  >  2:52076/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                             
CGCCTGCACATTGTCGGATGCGACGCCGGGGCGTCTTATCCGACCTACGATTTACTGCACCGCTTTACTCACCGAACTCATCGCGCTTAGCGCGTTGTTCAGGCGGTCCAGGTCATCCAGCTGCGCCTGAATCTCGGTGGTTTTTGCCAGCTTGTCCTGCTCATTCAGTCCTTTATTGCTTTTCACGCCCTGCATCTCTTTAAACAGCGCGAGCTGGCGAATCGGCACCAGTTGCAGGTCGCTGGAAGCGCGGAATGGCGCCCAGCCCAGGCGTTCCAGCACC  >  NC_000913/4323285‑4323567

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: