Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 548,669 A→T 5.7% T105S (ACT→TCT)  ylbE → NO‑induced DUF1116 protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913548,6690AT 5.7% 222.6 / 3.3 70T105S (ACT→TCT) ylbENO‑induced DUF1116 protein
Reads supporting (aligned to +/- strand):  ref base A (41/25);  new base T (2/2);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAT  >  NC_000913/548522‑548815
                                                                                                                                                   |                                                                                                                                                  
ggACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTa                                                                                                                                                    <  2:369396/148‑1 (MQ=255)
  aCCCATGAAAGGGGCGTGCGTTGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt                                                                                                                                                  <  1:299876/148‑1 (MQ=255)
   cccATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt                                                                                                                                                 <  8:44587/148‑1 (MQ=255)
    ccATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc                                                                                                                                                <  6:376607/148‑1 (MQ=255)
     cATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCt                                                                                                                                               >  5:332206/1‑148 (MQ=255)
       tGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGTGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAGGGGCGGTATTACTTCTGc                                                                                                                                             >  7:342374/1‑148 (MQ=255)
         aaaGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGc                                                                                                                                             >  2:2780/1‑146 (MQ=255)
            ggggCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTa                                                                                                                                        <  5:371261/148‑1 (MQ=255)
                  tgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGa                                                                                                                                  <  4:200479/148‑1 (MQ=255)
                       gggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtg                                                                                                                             >  8:132677/1‑148 (MQ=255)
                        ggCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAATCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtg                                                                                                                             <  5:328984/147‑1 (MQ=255)
                          cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtggtg                                                                                                                          >  5:29997/1‑148 (MQ=255)
                                       gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTg                                                                                                              <  5:87328/147‑1 (MQ=255)
                                        aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGAc                                                                                                            <  4:61915/148‑1 (MQ=255)
                                         aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGAcc                                                                                                           >  7:281437/1‑148 (MQ=255)
                                             tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGAccgccg                                                                                                       >  2:338956/1‑148 (MQ=255)
                                             tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACg                                                                                                       >  1:57022/1‑148 (MQ=255)
                                              gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCt                                                                                                                              <  5:272751/124‑1 (MQ=255)
                                              gggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCt                                                                                                                              >  6:272751/1‑124 (MQ=255)
                                                   aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAAcc                                                                                                 >  7:296391/1‑148 (MQ=255)
                                                                aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCaa                                                                                    >  5:180903/1‑148 (MQ=255)
                                                                     cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTc                                                                                >  8:226753/1‑147 (MQ=255)
                                                                             aCTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAgg                                                                        >  2:27501/1‑147 (MQ=255)
                                                                                ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATc                                                                    >  5:48258/1‑148 (MQ=255)
                                                                                  aGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtgg                                                                                                                            <  3:748/90‑1 (MQ=255)
                                                                                  aGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtgg                                                                                                                            >  4:748/1‑90 (MQ=255)
                                                                                    cAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCAAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGc                                                                 >  5:363107/1‑147 (MQ=255)
                                                                                     aGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCa                                                                >  8:229945/1‑147 (MQ=255)
                                                                                              tGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCg                                                      >  7:148370/1‑148 (MQ=255)
                                                                                                aaCTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGt                                                                                                               >  6:17088/1‑89 (MQ=255)
                                                                                                aaCTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGt                                                                                                               <  5:17088/89‑1 (MQ=255)
                                                                                                 aCTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGttt                                                   >  1:261962/1‑148 (MQ=255)
                                                                                                   ttCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTTCTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTg                                                 >  4:31801/1‑148 (MQ=255)
                                                                                                   ttCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTTCTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTg                                                 >  6:230104/1‑148 (MQ=255)
                                                                                                   ttCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGttttc                                                 >  7:218747/1‑147 (MQ=255)
                                                                                                    tCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTgg                                                >  6:28651/1‑148 (MQ=255)
                                                                                                    tCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTgg                                                >  4:291575/1‑148 (MQ=255)
                                                                                                    tCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTgg                                                >  4:21538/1‑148 (MQ=255)
                                                                                                     cATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTg                                                 <  6:48258/146‑1 (MQ=255)
                                                                                                       ttCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGgcgc                                             >  2:289450/1‑148 (MQ=255)
                                                                                                        tCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCt                                            >  7:69631/1‑148 (MQ=255)
                                                                                                            tGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACg                                        >  8:121007/1‑148 (MQ=255)
                                                                                                             gTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTAc                                         <  6:29997/146‑1 (MQ=255)
                                                                                                              tCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACgg                                       <  1:27501/147‑1 (MQ=255)
                                                                                                                aCCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGa                                    <  8:69631/148‑1 (MQ=255)
                                                                                                                      tgAATGCCGTCGGGCCAATGGGCGGTATTTCTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGt                              <  5:230104/148‑1 (MQ=255)
                                                                                                                      tgAATGCCGTCGGGCCAATGGGCGGTATTTCTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGt                              <  3:31801/148‑1 (MQ=255)
                                                                                                                        aaTGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTcc                            <  6:363107/148‑1 (MQ=255)
                                                                                                                          tGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTg                          <  1:2780/148‑1 (MQ=255)
                                                                                                                             cGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACt                       <  7:132677/148‑1 (MQ=255)
                                                                                                                              gTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCGGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTc                      <  6:332206/148‑1 (MQ=255)
                                                                                                                               tCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCg                     >  7:122561/1‑148 (MQ=255)
                                                                                                                               tCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCGTACGGCGAAGATGTCCTGACTCg                     >  1:54048/1‑148 (MQ=255)
                                                                                                                                    ccAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATc                >  2:109754/1‑148 (MQ=255)
                                                                                                                                       aTGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCt             >  5:182973/1‑148 (MQ=255)
                                                                                                                                        tGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTg            <  3:159851/148‑1 (MQ=255)
                                                                                                                                         gggCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTg            >  7:141593/1‑147 (MQ=255)
                                                                                                                                          ggCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTgg           <  3:291575/147‑1 (MQ=255)
                                                                                                                                           gCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTcc                            <  1:251152/129‑1 (MQ=255)
                                                                                                                                           gCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTcc                            >  2:251152/1‑129 (MQ=255)
                                                                                                                                                aTTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCTTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTAGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCCCAGGAGgcgc     >  8:211908/1‑147 (MQ=255)
                                                                                                                                                aTTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCaaa                                                              <  3:10629/90‑1 (MQ=255)
                                                                                                                                                aTTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCaaa                                                              >  4:10629/1‑90 (MQ=255)
                                                                                                                                                aTTACTTCTGCCAGTATGCCAATGCTGGTGGTTGAGAACGTGCCCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATgcgc     >  2:182245/1‑147 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCACCCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAt  >  1:67590/1‑148 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGttt                                                   >  6:109491/1‑99 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGttt                                                   <  5:109491/99‑1 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTgg           <  7:295686/139‑1 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTgg           >  8:295686/1‑139 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAt  >  6:143040/1‑148 (MQ=255)
                                                                                                                                                  tACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAt  >  1:187437/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTGACCGACGGCAACCGGGCGTACTGCAACCTCAACGAAGGTATCGGCAAAGTGATGCGTTTTGGCGCTTACGGCGAAGATGTCCTGACTCGCCATCGCTGGATGCGCGAT  >  NC_000913/548522‑548815

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