Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,862,791 T→C Y153H (TAT→CAT)  ygbK → FliA‑regulated DUF1537 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,862,7910TC98.3% 216.9 / ‑4.3 60Y153H (TAT→CAT) ygbKFliA‑regulated DUF1537 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (27/32);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGC  >  NC_000913/2862646‑2862935
                                                                                                                                                 |                                                                                                                                                
tGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCtt                                                                                                                                             >  2:177606/1‑151 (MQ=255)
  ctctCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTcc                                                                                                                                           >  1:682177/1‑151 (MQ=255)
   tctcGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTccc                                                                                                                                          >  1:794987/1‑151 (MQ=255)
         cacGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTg                                                                                                                                    <  2:380731/151‑1 (MQ=255)
         cacGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTg                                                                                                                                    >  2:194528/1‑151 (MQ=255)
             ccGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTg                                                                                                                                <  2:272949/151‑1 (MQ=255)
               gTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGaa                                                                                                                              <  2:71201/151‑1 (MQ=255)
                 ttACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAgc                                                                                                                            >  2:2632/1‑151 (MQ=255)
                   aCGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAgcgc                                                                                                                          <  1:347763/151‑1 (MQ=255)
                       tcttctCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATc                                                                                                                      >  1:164681/1‑151 (MQ=255)
                             ctcCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCcaca                                                                                                                  <  2:159962/149‑1 (MQ=255)
                              tcCGGCCCTGCCGGTCAACGGACGAACGGTCATTCAGGGGTTTTTGTTCGTAATGAACCAACGGCGGCCCGATTCCGGGAGGGCCCATCACCGGGAAATTCCCATGACGGACAGCCACCTTCCCGGCCTGGTTGAGCCCCAATCCACAggg                                                                                                               >  1:176573/1‑151 (MQ=255)
                               ccttcccTTCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATTAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGCCCGACAGCTATCTTCCCCTTCTGGTTGAAGCGCAATCcactt                                                                                                                 <  1:383310/144‑3 (MQ=255)
                               ccGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAg                                                                                                                 <  2:77974/148‑1 (MQ=255)
                                  gCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCt                                                                                                           <  1:166332/151‑1 (MQ=255)
                                    ccTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGtgctgc                                                                                                         <  2:697845/151‑1 (MQ=255)
                                     cTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTgcg                                                                                                        >  2:336110/1‑151 (MQ=255)
                                     cTGCCGGTCAACGAACGTACGGTTTATCAAGGGTATTTGTTCGTAATGAATCAACTGCTTGCCGAATCCCGGATGCGCCATCACCCGGTAAATCCCATCACCGACAGCCATCTTCCCCATCTGTTTGAAGCGCAAt                                                                                                                       <  1:25174/136‑1 (MQ=255)
                                        ccGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTgcggcg                                                                                                     <  2:840391/151‑1 (MQ=255)
                                          ggTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGcgtc                                                                                                   <  1:29437/151‑1 (MQ=255)
                                            tCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGcgtcgt                                                                                                 >  2:531209/1‑151 (MQ=255)
                                            tCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGcgtcgt                                                                                                 >  2:752314/1‑151 (MQ=255)
                                                cGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCg                                                                                             <  1:702219/151‑1 (MQ=255)
                                                   aCGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGc                                                                                           >  2:846961/1‑150 (MQ=255)
                                                      tACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACAt                                                                                       >  2:412489/1‑151 (MQ=255)
                                                      tACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACAt                                                                                       >  2:647726/1‑151 (MQ=255)
                                                       aCGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACAt                                                                                       <  1:299407/150‑1 (MQ=255)
                                                         ggTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGtt                                                                                    >  2:700704/1‑151 (MQ=255)
                                                          gTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGttt                                                                                   <  1:646618/151‑1 (MQ=255)
                                                           tttATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGtttt                                                                                  >  2:784639/1‑151 (MQ=255)
                                                             tATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTc                                                                                 <  2:164586/150‑1 (MQ=255)
                                                             tATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTc                                                                                 <  1:453355/150‑1 (MQ=255)
                                                                     gTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgt                                                                        <  1:470421/151‑1 (MQ=255)
                                                                     gTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgt                                                                        <  2:723202/151‑1 (MQ=255)
                                                                       aTTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGgtgt                                                                      >  2:693621/1‑151 (MQ=255)
                                                                           gTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGa                                                                   <  2:317505/150‑1 (MQ=255)
                                                                           gTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGAt                                                                  >  2:81219/1‑151 (MQ=255)
                                                                             tCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATg                                                                 <  2:525770/150‑1 (MQ=255)
                                                                              cGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGcc                                                               <  2:828141/151‑1 (MQ=255)
                                                                               gTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCg                                                              >  2:86463/1‑151 (MQ=255)
                                                                                 aaTGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGt                                                             <  2:519726/150‑1 (MQ=255)
                                                                                    gAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGt                                                         <  1:214932/151‑1 (MQ=255)
                                                                                        cAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCa                                                       >  1:801285/1‑149 (MQ=255)
                                                                                          aCTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAagag                                                   <  2:787807/151‑1 (MQ=255)
                                                                                              ctgGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTgg                                               >  2:730925/1‑151 (MQ=255)
                                                                                                 gCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTc                                            >  1:476197/1‑151 (MQ=255)
                                                                                                         cGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTAc                                      >  1:115829/1‑149 (MQ=255)
                                                                                                           ggATGCGCCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTGCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCa                                  <  1:540060/151‑1 (MQ=255)
                                                                                                               gcgcCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAagag                              <  1:412457/151‑1 (MQ=255)
                                                                                                               gcgcCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAaga                               <  2:54307/150‑1 (MQ=255)
                                                                                                                cgcCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAgg                             <  1:251830/151‑1 (MQ=255)
                                                                                                                 gcCATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAggg                            >  2:196386/1‑151 (MQ=255)
                                                                                                                  acATCACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGCTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGc                           >  1:536053/2‑151 (MQ=255)
                                                                                                                      cACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTAc                        <  2:197245/150‑1 (MQ=255)
                                                                                                                      cACCCGGTAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTAc                        <  1:733201/150‑1 (MQ=255)
                                                                                                                             tAAATCCCATGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTAcgc                >  2:357207/1‑151 (MQ=255)
                                                                                                                                      tGACCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTg        >  1:632035/1‑150 (MQ=255)
                                                                                                                                        aCCGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATg     <  2:739555/151‑1 (MQ=255)
                                                                                                                                         ccGACAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATg     <  1:671139/150‑1 (MQ=255)
                                                                                                                                            aCAGCCATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGTGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATgcgc  >  1:716188/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGCTCTCGACACGCCGTTTACGGTCTTCTCTCCGGCCCTGCCGGTCAACGGACGTACGGTTTATCAGGGGTATTTGTTCGTAATGAATCAACTGCTGGCCGAATCCGGGATGCGCCATCACCCGGTAAATCCCATGACCGACAGCTATCTTCCCCGTCTGGTTGAAGCGCAATCCACAGGGCGCTGCGGCGTCGTTTCGGCACATGTTTTCGAACAAGGTGTGGATGCCGTTCGTCAAGAGCTGGCTCGCTTACAGCAAGAGGGCTACCGCTACGCGGTGCTTGATGCGC  >  NC_000913/2862646‑2862935

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