Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,777,919 G→A 5.5% G174D (GGC→GAC)  lldP → L‑lactate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,777,9190GA 5.5% 235.8 / 2.3 73G174D (GGC→GAC) lldPL‑lactate permease
Reads supporting (aligned to +/- strand):  ref base G (29/40);  new base A (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

GGCGCGTTCCTTGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACGGCTGGCGCGGTATCAAAGAGACG  >  NC_000913/3777774‑3778067
                                                                                                                                                 |                                                                                                                                                    
ggCGCGTTCCTTGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCa                                                                                                                                                    <  5:200725/148‑1 (MQ=255)
 gcggGTTCCTTGAAGGAGCCGCAGGCTTTGGCGCACCAGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTGTTAACCCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGGGGCATTTGGTGCGCTGGGCAtt                                                                                                                                                  >  3:79406/1‑149 (MQ=255)
 gcgcGTTCCTTGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCAtt                                                                                                                                                  <  5:112851/149‑1 (MQ=255)
         cTTGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCt                                                                                                                                          >  6:275438/1‑149 (MQ=255)
          ttGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTg                                                                                                                                         >  1:356519/1‑149 (MQ=255)
          ttGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTg                                                                                                                                         >  1:44901/1‑149 (MQ=255)
            gAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTg                                                                                                                                         <  2:356519/147‑1 (MQ=255)
               ggAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGc                                                                                                                                    <  1:78641/149‑1 (MQ=255)
                    cGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGac                                                                                                                               <  8:173069/149‑1 (MQ=255)
                     gCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGaca                                                                                                                              <  2:347560/149‑1 (MQ=255)
                     gCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGa                                                                                                                                >  2:138976/1‑147 (MQ=255)
                      cAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGacag                                                                                                                             <  2:44901/149‑1 (MQ=255)
                          cTTTGGCGCACCGGTAGCAATTACCGACGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGacaggtaa                                                                                                                         <  3:42139/149‑1 (MQ=255)
                             tGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGacaggtaacag                                                                                                                      >  1:338354/1‑149 (MQ=255)
                             tGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGacaggtaaca                                                                                                                       <  7:204922/148‑1 (MQ=255)
                               gcgcACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGacaggtaacaggt                                                                                                                    <  2:323317/149‑1 (MQ=255)
                                 gcACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCa                                                                                                                                                    <  2:192030/115‑1 (MQ=255)
                                 gcACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCa                                                                                                                                                    >  1:192030/1‑115 (MQ=255)
                                  cACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATc                                                                                                                 >  6:226362/1‑149 (MQ=255)
                                           cAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCttt                                                                                                        <  5:275438/149‑1 (MQ=255)
                                             aTTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCttt                                                                                                        >  4:84090/1‑147 (MQ=255)
                                             aTTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTga                                                                                                      <  4:307095/149‑1 (MQ=255)
                                                aCCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTgaga                                                                                                    <  5:331428/148‑1 (MQ=255)
                                                aCCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTgaga                                                                                                    <  2:76682/148‑1 (MQ=255)
                                                aCCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTgaga                                                                                                    <  5:226362/148‑1 (MQ=255)
                                                  cgccgcATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTg                                                                                                 <  5:326736/149‑1 (MQ=255)
                                                  cgccgcATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTg                                                                                                 <  1:138976/149‑1 (MQ=255)
                                                        aTTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGa                                                                                           <  5:123957/149‑1 (MQ=255)
                                                         ttGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAt                                                                                          >  8:235156/1‑149 (MQ=255)
                                                         ttGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAt                                                                                          >  7:293943/1‑149 (MQ=255)
                                                         ttGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAt                                                                                          >  7:288845/1‑149 (MQ=255)
                                                         ttGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAt                                                                                          >  5:4687/1‑149 (MQ=255)
                                                               gTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAtggtg                                                                                      <  7:65163/147‑1 (MQ=255)
                                                                tCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTg                                                                                                       >  2:81651/1‑129 (MQ=255)
                                                                              aaaCCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGt                                                                      <  7:17505/148‑1 (MQ=255)
                                                                                aCCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTa                                                                   <  8:293943/149‑1 (MQ=255)
                                                                                aCCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGGGGGGCGGCAGCTACCGTTTa                                                                   <  8:288845/149‑1 (MQ=255)
                                                                                 ccGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGggc                                                                                  <  4:258941/133‑1 (MQ=255)
                                                                                 ccGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGggc                                                                                  >  3:258941/1‑133 (MQ=255)
                                                                                   gCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGa                                                                <  1:193265/149‑1 (MQ=255)
                                                                                         cgccgcCGGGCGGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCAtt                                                           <  4:79406/148‑1 (MQ=255)
                                                                                          gccgccGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAtggt                                                                                       >  4:12710/1‑119 (MQ=255)
                                                                                          gccgccGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGAtggt                                                                                       <  3:12710/119‑1 (MQ=255)
                                                                                          gccgccGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTAt                                                         >  8:129826/1‑149 (MQ=255)
                                                                                          gccgccGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGACATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTa                                                                          <  7:246541/132‑1 (MQ=255)
                                                                                          gccgccGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGACATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTa                                                                          >  8:246541/1‑132 (MQ=255)
                                                                                           ccgccgGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATc                                                        <  4:24636/149‑1 (MQ=255)
                                                                                                gggCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCt                                                   <  3:84090/149‑1 (MQ=255)
                                                                                                 ggCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTg                                                  <  8:102850/149‑1 (MQ=255)
                                                                                                  gCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGACATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGt                                                 >  6:98377/1‑149 (MQ=255)
                                                                                                    tgtgCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTc                                               >  8:47009/1‑149 (MQ=255)
                                                                                                           gATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCa                                        >  2:290818/1‑149 (MQ=255)
                                                                                                               gTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATgg                                     <  7:203696/148‑1 (MQ=255)
                                                                                                                ttAACACCGCGCCAGTGGCATTTGTTGCGATGGGCATTCCAATCCTGCTTGCCGGACAGGTAAGAGGTATCGACTGCTTTGAGATTGGTCAGATGGTGGGGCTGCAGCTACCGTTTAt                                                                  <  6:100805/118‑1 (MQ=255)
                                                                                                                ttAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTACAATCCTGGTTGACGGTCAGGTAACAGGTATCTACAGCTTTGAGATTGGTGAGATGGTGTGGCGGCAGATACCGTTTAt                                                                  >  5:100805/1‑118 (MQ=255)
                                                                                                                    cacCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATc                                                        >  8:315796/1‑124 (MQ=255)
                                                                                                                    cacCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATc                                                        <  7:315796/124‑1 (MQ=255)
                                                                                                                     acCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTAt                              >  7:259961/1‑149 (MQ=255)
                                                                                                                      ccGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATg                             <  7:25086/149‑1 (MQ=255)
                                                                                                                      ccGCGCCAGTGGCATTTGGTGCGATGGACATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATg                             <  5:98377/149‑1 (MQ=255)
                                                                                                                          gcCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACg                         >  3:325860/1‑149 (MQ=255)
                                                                                                                              gaggCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCt                                                   >  4:18922/3‑119 (MQ=255)
                                                                                                                              gaggCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCt                                                   <  3:18922/117‑1 (MQ=255)
                                                                                                                               tGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTgaga                                                                                                    >  2:289795/1‑69 (MQ=255)
                                                                                                                                   aTTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGAc                          >  5:156289/1‑139 (MQ=255)
                                                                                                                                   aTTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGAc                          <  6:156289/139‑1 (MQ=255)
                                                                                                                                    tttGGTGCGAAGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACGGCTGGCGCgg               >  4:13668/1‑149 (MQ=255)
                                                                                                                                      tGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACGGCTGGCGCGGt              <  1:290818/148‑1 (MQ=255)
                                                                                                                                       ggTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGGCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGAGCATGGCGATTATGGACGGCTGGCGCGGTa             <  7:235156/148‑1 (MQ=255)
                                                                                                                                         tGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATg                             <  1:176379/130‑1 (MQ=255)
                                                                                                                                         tGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATg                             >  2:176379/1‑130 (MQ=255)
                                                                                                                                         tGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACACGTATCGACAGCTTTGAGATTCGTCAGATGGTGGCGCGGCAACTACCGTTTATGACCATTATCGTGCTGTTATGGATCATGGCGATTATGGACGGCTGGCGCGGTATc           >  1:133319/1‑148 (MQ=255)
                                                                                                                                             aTGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACGGCTGGCGCGGTATCAAaga      <  7:129826/149‑1 (MQ=255)
                                                                                                                                                 gCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACGGCTGGCGCGGTATCAAAGAGACg  >  8:64089/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGCGCGTTCCTTGAAGGAGCCGCAGGCTTTGGCGCACCGGTAGCAATTACCGCCGCATTGCTGGTCGGCCTGGGTTTTAAACCGCTGTACGCCGCCGGGCTGTGCCTGATTGTTAACACCGCGCCAGTGGCATTTGGTGCGATGGGCATTCCAATCCTGGTTGCCGGACAGGTAACAGGTATCGACAGCTTTGAGATTGGTCAGATGGTGGGGCGGCAGCTACCGTTTATGACCATTATCGTGCTGTTCTGGATCATGGCGATTATGGACGGCTGGCGCGGTATCAAAGAGACG  >  NC_000913/3777774‑3778067

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: ‑