Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,511,746 (A)8→9 coding (599/906 nt) yhfZ ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,511,7381.A96.8% 195.5 / ‑1.0 62coding (607/906 nt)yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (32/28);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑A‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGC  >  NC_000913/3511600‑3511871
                                                                                                                                           |                                                                                                                                         
gcggcgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGcc                                                                                                                                            <  2:1286273/139‑1 (MQ=255)
 cggcgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a                                                                                                                                         <  2:401977/139‑1 (MQ=255)
   gcgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aa                                                                                                                                     >  1:1288861/1‑139 (MQ=255)
   gcgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aa                                                                                                                                     >  1:1296221/1‑139 (MQ=255)
    cgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aaa                                                                                                                                    >  2:1291525/1‑139 (MQ=255)
    cgTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCAAGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aaa                                                                                                                                    >  1:1141537/1‑139 (MQ=255)
           tctAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aaaaaaaaca                                                                                                                             <  1:108984/139‑4 (MQ=255)
           tctAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aaaaaaaaca                                                                                                                             <  2:1157589/139‑4 (MQ=255)
                 cccAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTAATAAGATAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTc                                                                                                                       <  1:1291525/139‑1 (MQ=255)
                    aGCATGGTCATTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACCATCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATg                                                                                                                    <  2:1296221/139‑1 (MQ=255)
                    aGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATg                                                                                                                    <  2:1302161/139‑1 (MQ=255)
                         ggTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGAtttt                                                                                                               >  1:1161380/1‑139 (MQ=255)
                          gTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aaaaaaaaca                                                                                                                             <  2:730596/124‑4 (MQ=255)
                          gTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a‑‑aaaaaaaaca                                                                                                                             >  1:730596/1‑121 (MQ=255)
                           tCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCgg                                                                                                                         >  2:574384/1‑127 (MQ=255)
                           tCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCgg                                                                                                                         <  1:574384/127‑1 (MQ=255)
                              gTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGAt                                                                                                          <  2:1095908/139‑1 (MQ=255)
                                                  ttCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑AAGAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAg                                                                                       >  2:960200/1‑139 (MQ=255)
                                                        aTCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                >  1:587689/1‑139 (MQ=255)
                                                           aCCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACgcgc                                                                             >  1:1105177/1‑139 (MQ=255)
                                                            ccGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAATAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCt                                                                            >  1:1353228/1‑139 (MQ=255)
                                                            ccGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGACTTTCTGATCCGCCGAACgg                                                                                               <  2:1074077/120‑1 (MQ=255)
                                                            ccGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑A‑‑AAAAAAAGCATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                               >  1:1074077/1‑120 (MQ=255)
                                                             cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                           >  1:653508/1‑139 (MQ=255)
                                                               cATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCa                                                                         <  2:1141537/139‑1 (MQ=255)
                                                               cATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCa                                                                         <  2:1105177/139‑1 (MQ=255)
                                                                           cTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACt                                                             <  2:386910/139‑1 (MQ=255)
                                                                            tttAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTc                                                            <  2:653508/139‑1 (MQ=255)
                                                                                   aTGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                               >  1:480946/1‑97 (MQ=255)
                                                                                   aTGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                               <  2:480946/97‑1 (MQ=255)
                                                                                    tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                    >  2:1307037/1‑139 (MQ=255)
                                                                                     gCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                           <  2:632165/99‑1 (MQ=255)
                                                                                     gCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                           >  1:632165/1‑99 (MQ=255)
                                                                                      cGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                    <  1:454613/89‑1 (MQ=255)
                                                                                      cGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                    >  2:454613/1‑89 (MQ=255)
                                                                                         tGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTc                                                                                                                       <  1:167315/67‑1 (MQ=255)
                                                                                         tGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTc                                                                                                                       >  2:167315/1‑67 (MQ=255)
                                                                                             aaCTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCa                                           >  2:438598/1‑139 (MQ=255)
                                                                                               ctctCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAAc                                         <  2:14238/139‑1 (MQ=255)
                                                                                                tctcGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                <  2:1355782/99‑1 (MQ=255)
                                                                                                tctcGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                >  1:1355782/1‑99 (MQ=255)
                                                                                                tctcGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACt                                        >  2:1247360/1‑139 (MQ=255)
                                                                                                   cGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTgg                                      >  1:578071/1‑139 (MQ=255)
                                                                                                   cGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTg                                                                   >  2:400949/1‑109 (MQ=255)
                                                                                                   cGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTg                                                                   <  1:400949/109‑1 (MQ=255)
                                                                                                   cGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGt                                     >  2:518974/1‑139 (MQ=255)
                                                                                                    gTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTg                                    <  2:1353228/139‑1 (MQ=255)
                                                                                                     tGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGc                                   >  1:1363468/1‑139 (MQ=255)
                                                                                                         aaGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                               <  1:1247360/139‑1 (MQ=255)
                                                                                                         aaGAGAGATCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                               <  2:587689/139‑1 (MQ=255)
                                                                                                            agagaTCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                <  2:1512626/87‑1 (MQ=255)
                                                                                                            agagaTCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                >  1:1512626/1‑87 (MQ=255)
                                                                                                                aTCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                        >  1:143842/1‑139 (MQ=255)
                                                                                                                 tCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAg                       <  2:1363468/139‑1 (MQ=255)
                                                                                                                 tCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAg                       >  1:1367328/1‑139 (MQ=255)
                                                                                                                 tCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAg                       >  1:1191083/1‑139 (MQ=255)
                                                                                                                 tCGACTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCACTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAg                       >  1:853779/1‑139 (MQ=255)
                                                                                                                    aCTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGgtg                    >  2:1270772/1‑139 (MQ=255)
                                                                                                                     cTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                >  1:557855/1‑78 (MQ=255)
                                                                                                                     cTCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                                <  2:557855/78‑1 (MQ=255)
                                                                                                                      tCGTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCACTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGtgtg                  <  2:853779/139‑1 (MQ=255)
                                                                                                                        gTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATccg                                                                                                       <  1:1097506/52‑1 (MQ=255)
                                                                                                                        gTTCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATccg                                                                                                       >  2:1097506/1‑52 (MQ=255)
                                                                                                                         ttCCACATCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCgg               >  1:1028817/1‑139 (MQ=255)
                                                                                                                               aTCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                        >  1:1322121/1‑124 (MQ=255)
                                                                                                                               aTCACTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                        <  2:1322121/124‑1 (MQ=255)
                                                                                                                                   cTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAATGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACtac     >  1:968990/1‑137 (MQ=255)
                                                                                                                                   cTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTc                                                                          >  2:1181186/1‑70 (MQ=255)
                                                                                                                                   cTACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTc                                                                          <  1:1181186/70‑1 (MQ=255)
                                                                                                                                    tACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAAc                                                                                                 >  1:157829/1‑46 (MQ=255)
                                                                                                                                    tACCGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAAc                                                                                                 <  2:157829/46‑1 (MQ=255)
                                                                                                                                      ccGCCAA‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAgc  >  1:1282311/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                         
GCGGCGTCGCCTCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑A‑‑AAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCAGC  >  NC_000913/3511600‑3511871

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