Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 455,731 A→G D200G (GAC→GGC)  tig → peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913455,7310AG96.7% 197.9 / ‑2.0 61D200G (GAC→GGC) tigpeptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (26/33);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

AAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGCCGGAACTG  >  NC_000913/455587‑455876
                                                                                                                                                |                                                                                                                                                 
aagaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGGCTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATGTCGTACTGGCGATGGGCCAGGGTCGGATGATCCCGGGCTTTGAAGGCGGTAt                                                                                                                                             <  1:79921/151‑1 (MQ=255)
  gaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATc                                                                                                                                            >  1:849100/1‑150 (MQ=255)
  gaaaAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCa                                                                                                                                           <  1:553200/151‑1 (MQ=255)
     aaaGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAg                                                                                                                                        <  2:674504/151‑1 (MQ=255)
     aaaGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGGTCGGTGGACGGCGAAGAGTTCGAAGGCGGTAAAGAGACTGATTTCGTAATGGCGATGGGCCAGGGTCGAATGATCCCGGGCTTTGAAGACGGTATCaa                                                                                                                                          >  1:365470/1‑149 (MQ=255)
      aaGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAgg                                                                                                                                       <  2:199940/151‑1 (MQ=255)
       aGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAgg                                                                                                                                       >  1:138137/1‑150 (MQ=255)
         aCGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGcc                                                                                                                                     <  1:701035/150‑1 (MQ=255)
          cGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCcac                                                                                                                                   >  2:857889/1‑151 (MQ=255)
           ggCGCTGTTGAAGCAGAAGTCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCcaca                                                                                                                                  <  2:617042/151‑1 (MQ=255)
             cgcTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACaaa                                                                                                                                <  2:430237/151‑1 (MQ=255)
             cgcTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACaaa                                                                                                                                <  1:580792/151‑1 (MQ=255)
               cTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGc                                                                                                                              >  1:384253/1‑151 (MQ=255)
                tGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCt                                                                                                                             >  1:101692/1‑151 (MQ=255)
                tGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCt                                                                                                                             >  2:472068/1‑151 (MQ=255)
                     aaGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCg                                                                                                                         >  2:460096/1‑150 (MQ=255)
                      aGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGCAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGa                                                                                                                        >  1:820762/1‑150 (MQ=255)
                           aagaCCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGt                                                                                                                   <  1:555776/150‑1 (MQ=255)
                                      aCCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACg                                                                                                       <  2:220897/151‑1 (MQ=255)
                                                         cTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAAt                                                                                     >  1:569338/1‑150 (MQ=255)
                                                                 ggCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGc                                                                              >  2:316689/1‑149 (MQ=255)
                                                                 ggCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGGCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCAAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCa                                                                                 >  1:543854/1‑146 (MQ=255)
                                                                    gAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGaaa                                                                         <  2:573530/151‑1 (MQ=255)
                                                                           tCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGa                                                                   >  1:523837/1‑150 (MQ=255)
                                                                             gAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGaa                                                                  >  2:445892/1‑149 (MQ=255)
                                                                              aaGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAgg                                                               <  1:116254/151‑1 (MQ=255)
                                                                                ggcggcAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTa                                                             <  2:820818/151‑1 (MQ=255)
                                                                                ggCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTa                                                             <  1:264443/151‑1 (MQ=255)
                                                                                 gCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTaa                                                            <  2:849159/151‑1 (MQ=255)
                                                                                  cGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTaa                                                            >  1:694869/1‑150 (MQ=255)
                                                                                        aGCGTTTGATTTCGTACTGGCGATGGGCCAGGGCCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAagca                                                        >  1:162727/1‑148 (MQ=255)
                                                                                        aGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCg                                                     >  2:785185/1‑151 (MQ=255)
                                                                                         gCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGa                                                    <  1:408251/151‑1 (MQ=255)
                                                                                           gTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGCGACCTTCCCGGAAGAATACCACGCAGAAAACCTAAAAGGTAAAGCAGCGa                                                    >  1:606590/1‑149 (MQ=255)
                                                                                               gATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAAtt                                                <  2:354744/149‑1 (MQ=255)
                                                                                               gATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCg                                              >  1:530207/1‑151 (MQ=255)
                                                                                                aTTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGc                                             <  2:276605/151‑1 (MQ=255)
                                                                                                   tCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTa                                           <  2:144842/150‑1 (MQ=255)
                                                                                                     gTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCa                                        <  2:230003/151‑1 (MQ=255)
                                                                                                         tGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCt                                    >  2:543003/1‑151 (MQ=255)
                                                                                                          ggCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAAcc                                     <  2:458764/149‑1 (MQ=255)
                                                                                                          ggCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCt                                    <  1:594925/150‑1 (MQ=255)
                                                                                                          ggCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCt                                    <  2:126882/150‑1 (MQ=255)
                                                                                                           gCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTga                                  <  2:784908/151‑1 (MQ=255)
                                                                                                                gggCCAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTgaagaa                              >  2:434885/1‑150 (MQ=255)
                                                                                                                   ccAGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGtt                          <  2:677006/151‑1 (MQ=255)
                                                                                                                     aGGGTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGa                        <  1:375898/151‑1 (MQ=255)
                                                                                                                      gggTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGaa                       <  1:539869/151‑1 (MQ=255)
                                                                                                                      gggTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTACCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGa                        >  2:635014/1‑150 (MQ=255)
                                                                                                                       ggTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGa                        >  1:254536/1‑149 (MQ=255)
                                                                                                                        gTCGTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGa                        >  1:318152/1‑148 (MQ=255)
                                                                                                                           gTATGATCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACAACGCAGAAAACCTGAAAGGTAAAGCAGCGACATTCGCTATCAACCTGAAGAAAGTTGAAGAGCg                  >  1:48579/1‑151 (MQ=255)
                                                                                                                             aTGATCCCGGACTTTGAATGCGGTATCAAAGGCCACAAAGCTGGCGGAGAGTTCACCATGGACATGATCTTCCCGGAAGAATACCACGCAAAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAagtg                    <  1:310660/147‑3 (MQ=255)
                                                                                                                                aTCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAAc             <  1:113774/151‑1 (MQ=255)
                                                                                                                                aTCCCGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAAc             >  2:309936/1‑151 (MQ=255)
                                                                                                                                  cccGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACt            <  2:586350/150‑1 (MQ=255)
                                                                                                                                   ccGGGCTTTGAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGc          <  1:93578/151‑1 (MQ=255)
                                                                                                                                            gAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGCCGGAactg  >  2:602187/1‑150 (MQ=255)
                                                                                                                                            gAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGCCGGAactg  <  2:561178/150‑1 (MQ=255)
                                                                                                                                            gAAGGCGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGCCGGAactg  <  1:530255/150‑1 (MQ=255)
                                                                                                                                            gAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGCCGGAactg  >  2:451743/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AAGAAAAAGACGGCGCTGTTGAAGCAGAAGACCGCGTAACCATCGACTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATCGACGTGACCTTCCCGGAAGAATACCACGCAGAAAACCTGAAAGGTAAAGCAGCGAAATTCGCTATCAACCTGAAGAAAGTTGAAGAGCGTGAACTGCCGGAACTG  >  NC_000913/455587‑455876

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