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.A93.7% 240.7 / 7.8 79coding (607/906 nt)yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base A (41/33);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCA  >  NC_000913/3511611‑3511869
                                                                                                                                  |                                                                                                                                    
tctAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑aaaaaaaaaca                                                                                                                          <  1:978340/139‑4 (MQ=255)
         aGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATg                                                                                                                  >  2:845381/1‑139 (MQ=255)
          gCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGa                                                                                                                 >  2:495512/1‑139 (MQ=255)
           cATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCa                                                                                                                    <  2:309491/135‑1 (MQ=255)
           cATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCa                                                                                                                    >  1:309491/1‑135 (MQ=255)
                   gTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGt                                                                                                                      >  2:1331244/1‑125 (MQ=255)
                   gTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGt                                                                                                                      <  1:1331244/125‑1 (MQ=255)
                   gTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a                                                                                                                                    <  2:800744/110‑1 (MQ=255)
                   gTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑a                                                                                                                                    >  1:800744/1‑110 (MQ=255)
                         tttCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGTTCATGATTTTCTGATccgccg                                                                                                  >  1:277860/1‑139 (MQ=255)
                         tttCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  >  1:118652/1‑139 (MQ=255)
                              gCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             <  1:495512/139‑1 (MQ=255)
                                  ccacGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTAt                                                                                         >  1:1452940/1‑139 (MQ=255)
                                            gATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             <  2:941385/125‑1 (MQ=255)
                                            gATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             >  1:941385/1‑125 (MQ=255)
                                                  cGCATCGTCATCGCCTTTAACAATACGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGACA‑‑AAAAAAAACATCGGCCATGATTTTCTGATCCGCCGAACGGTTATCGAGCCCCACCCGCtt                                                                         >  2:1122002/1‑139 (MQ=255)
                                                  cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCgg                                                                                                                       <  2:148390/93‑1 (MQ=255)
                                                  cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCgg                                                                                                                       >  1:148390/1‑93 (MQ=255)
                                                  cGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCtt                                                                         >  2:1639473/1‑139 (MQ=255)
                                                    cATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCa                                                                       >  2:543766/1‑139 (MQ=255)
                                                                 tttAACCAGGCGTTGTAAACTCTCGTGATAAGGGAAATCGACTCGTTCCACATCACTACCTCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTc                                                          <  1:1122002/139‑1 (MQ=255)
                                                                    aaCAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAc                                                                               >  1:1131610/1‑115 (MQ=255)
                                                                    aaCAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAc                                                                               <  2:1131610/115‑1 (MQ=255)
                                                                     aCAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGAt                                                                                                        <  1:1790772/89‑1 (MQ=255)
                                                                     aCAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGAt                                                                                                        >  2:1790772/1‑89 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  >  2:39074/1‑91 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  <  1:39074/91‑1 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             >  2:313467/1‑96 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             <  1:313467/96‑1 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                  <  2:277860/139‑1 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                  >  1:76016/1‑139 (MQ=255)
                                                                         tGCGTTGTAAACTCTCGTCATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                  <  2:1235117/139‑1 (MQ=255)
                                                                          gCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATAACTACCGCAA‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCTGAACGTTTATCCAGCCCCACGCGCTTCACGTTTGGGGACTCGCATTTACg                                                 >  1:314198/1‑139 (MQ=255)
                                                                              tGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGaa                                                                                                <  1:1034466/88‑1 (MQ=255)
                                                                              tGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGaa                                                                                                >  2:1034466/1‑88 (MQ=255)
                                                                              tGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCaa                                             >  2:1810063/1‑139 (MQ=255)
                                                                              tGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCaa                                             >  2:77066/1‑139 (MQ=255)
                                                                                        cGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  <  1:800725/76‑1 (MQ=255)
                                                                                        cGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  >  2:800725/1‑76 (MQ=255)
                                                                                            aTAAGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                               <  1:845381/139‑1 (MQ=255)
                                                                                              aaGAGAGATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                             <  1:77066/139‑1 (MQ=255)
                                                                                               agagagATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATg                                                                                                                  >  2:1262820/1‑53 (MQ=255)
                                                                                               agagagATCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATg                                                                                                                  <  1:1262820/53‑1 (MQ=255)
                                                                                                 agagaTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAATGGTTATCCAGCCCCAcg                                                                              >  2:481376/1‑87 (MQ=255)
                                                                                                 agagaTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAATGGTTATCCAGCCCCAcg                                                                              <  1:481376/87‑1 (MQ=255)
                                                                                                     aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGt                       >  2:1016432/1‑139 (MQ=255)
                                                                                                     aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAACTCAACTGGTGCTCGCCAACGTa                      >  2:1296818/1‑139 (MQ=255)
                                                                                                     aTCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTa                      >  1:997584/1‑139 (MQ=255)
                                                                                                      tCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCa                                                                                                                    <  2:497483/44‑1 (MQ=255)
                                                                                                      tCGACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCa                                                                                                                    >  1:497483/1‑44 (MQ=38)
                                                                                                        gACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                              <  1:737762/80‑1 (MQ=255)
                                                                                                        gACTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCAcg                                                                              >  2:737762/1‑80 (MQ=255)
                                                                                                         aCTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGgtg                  >  2:1069697/1‑139 (MQ=255)
                                                                                                         aCTCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGgtg                  >  1:709843/1‑139 (MQ=255)
                                                                                                           tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCa                                                                                                                    <  1:744644/39‑1 (MQ=255)
                                                                                                           tCGTTCCACATCACTACCGCCAAAAAAAAACATCGGTCa                                                                                                                    >  2:744644/1‑39 (MQ=255)
                                                                                                             gTTCCACATCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCg              <  2:118652/139‑1 (MQ=255)
                                                                                                              ttCCACATCACTACCGCCAAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCg              <  1:1016432/139‑1 (MQ=255)
                                                                                                                  acaTCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGcccc         <  1:543766/139‑1 (MQ=255)
                                                                                                                   caTCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCa        <  1:1639473/139‑1 (MQ=255)
                                                                                                                   caTCACTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCa        >  1:268487/1‑139 (MQ=255)
                                                                                                                      cACTACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTc    >  1:1438932/1‑139 (MQ=255)
                                                                                                                       aCTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  >  2:1511419/1‑45 (MQ=255)
                                                                                                                       aCTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  <  1:1511419/45‑1 (MQ=255)
                                                                                                                       aCTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTc    >  2:1810911/1‑139 (MQ=255)
                                                                                                                        cTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  >  2:1338721/1‑44 (MQ=38)
                                                                                                                        cTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATccgccg                                                                                                  <  1:1338721/44‑1 (MQ=255)
                                                                                                                        cTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCa                                                                                <  2:1094034/62‑1 (MQ=255)
                                                                                                                        cTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCa                                                                                >  1:1094034/1‑62 (MQ=255)
                                                                                                                        cTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                               >  2:1661791/1‑111 (MQ=255)
                                                                                                                        cTACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTc                               <  1:1661791/111‑1 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAATGGTTATCCAGCCCCATGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                  <  1:1057941/91‑1 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             >  1:197425/1‑48 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACgg                                                                                             <  2:197425/48‑1 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTAc                                                  >  2:1057941/1‑91 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                             >  2:1039226/1‑112 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGc                             <  1:1039226/112‑1 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGccc          <  2:1614716/131‑1 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGccc          >  1:1614716/1‑131 (MQ=255)
                                                                                                                         tACCGCCAAAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCa  >  1:295792/1‑139 (MQ=255)
                                                                                                                         tACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCa                                                                                <  1:932333/60‑1 (MQ=255)
                                                                                                                         tACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCa                                                                                >  2:932333/1‑60 (MQ=255)
                                                                                                                                  |                                                                                                                                    
TCTAATCCCAGCATGGTCAGTTCGTTTTCCGCCACCACGTTCCAGATCACCGCATCGACATCGCCTTTAACAATGCGTTGTAAACTCTCGTGATAAGAGAGATCGACTCGTTCCACATCACTACCGCC‑‑AAAAAAAACATCGGTCATGATTTTCTGATCCGCCGAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTGCTCGCCAACGTAGGTGTGCGGCCCCAACTCCA  >  NC_000913/3511611‑3511869

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: