Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,148,639 A→G 59.1% Y340C (TAC→TGC)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,148,6390AG59.1% 52.7 / 103.1 93Y340C (TAC→TGC) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base G (36/19);  total (57/36)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

GCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAA  >  NC_000913/1148542‑1148737
                                                                                                 |                                                                                                  
gCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATAcc                                                                                                  <  2:708262/100‑1 (MQ=255)
gCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATAcc                                                                                                  >  1:2037133/1‑100 (MQ=255)
  cAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCa                                                                                                >  1:389169/1‑100 (MQ=255)
  cAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCa                                                                                                >  2:616422/1‑100 (MQ=255)
  cAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCa                                                                                                >  2:2259808/1‑100 (MQ=255)
     tCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCCGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGc                                                                                              <  1:1289319/99‑1 (MQ=255)
      gatctaGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTg                                                                                            <  1:1056969/95‑1 (MQ=255)
      cAGCGAGCTTTTGCGGTCGCGATTGAACGGGCAGTGCAGGCGGTGCAGCGACAAGTTCCGCAGCGAATTGCCGCTCGCCTGGAAGCTGTATGCCCAGCt                                                                                             <  1:520415/99‑1 (MQ=255)
      cAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTg                                                                                            >  2:1909467/1‑100 (MQ=255)
        gCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGt                                                                                          >  1:1829271/1‑100 (MQ=255)
        gCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGt                                                                                          >  2:1361874/1‑100 (MQ=255)
         cGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGt                                                                                          >  1:1011059/1‑99 (MQ=255)
         cGAGCTTTTGCGGACGAGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGTGTCTGTATACCCAGCTGttt                                                                                         >  1:2425683/1‑100 (MQ=255)
           aGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGtttt                                                                                       <  2:389171/100‑1 (MQ=255)
           aGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGtttt                                                                                       <  2:425052/100‑1 (MQ=255)
           aGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTgagatc                                                                                       >  2:1056973/1‑95 (MQ=255)
            gCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTg                                                                                      <  1:1138815/100‑1 (MQ=255)
            gCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTg                                                                                      <  1:1840199/100‑1 (MQ=255)
             cTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGcgtttcgtc                                                                                      >  1:2572440/1‑91 (MQ=255)
               tttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgc                                                                                   >  2:2497557/1‑100 (MQ=255)
               tttGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAgc                                                                                   <  2:1391421/100‑1 (MQ=255)
                ttGCGGTCGCGATTGAACAGGCCGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgct                                                                                  >  1:1503710/1‑100 (MQ=255)
                  gCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAgctgc                                                                                <  2:182566/100‑1 (MQ=255)
                   cGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAgctgct                                                                               <  2:1012844/100‑1 (MQ=255)
                    ggTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAgctgct                                                                               >  1:1129848/1‑99 (MQ=255)
                    ggTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGctgctg                                                                              >  2:639200/1‑100 (MQ=255)
                    ggTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGctgctg                                                                              <  2:1227805/100‑1 (MQ=255)
                     gTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTgg                                                                             >  1:2215797/1‑100 (MQ=255)
                      tCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGa                                                                            >  2:802745/1‑100 (MQ=255)
                      tCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGa                                                                            <  2:856695/100‑1 (MQ=255)
                        gcgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACg                                                                          >  1:1881128/1‑100 (MQ=255)
                        gcgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCGGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACg                                                                          <  1:1602111/100‑1 (MQ=255)
                         cgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTGAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACgg                                                                         >  1:1410071/1‑100 (MQ=255)
                         cgATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACgg                                                                         <  1:2401753/100‑1 (MQ=255)
                            ttGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTgg                                                                      <  2:1694283/100‑1 (MQ=255)
                            ttGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTgg                                                                      <  1:1514720/100‑1 (MQ=255)
                             tGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGc                                                                     >  1:1213240/1‑100 (MQ=255)
                             tGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGc                                                                     <  1:855767/100‑1 (MQ=255)
                              gAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCa                                                                    >  2:2508231/1‑100 (MQ=255)
                              gAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCa                                                                    >  1:1936993/1‑100 (MQ=255)
                              gAACAGGCAGTGCAGGCGGTGCAGCGACAAGTGCCTCAGCGAATTGCCGCGCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCa                                                                    >  1:2508051/1‑100 (MQ=255)
                               aaCAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCa                                                                    <  1:757905/99‑1 (MQ=255)
                                   ggCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGc                                                               >  2:1000513/1‑100 (MQ=255)
                                   ggCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGc                                                               >  1:2006353/1‑100 (MQ=255)
                                    ctaGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCg                                                              <  1:1329228/98‑1 (MQ=255)
                                     cAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCgg                                                             >  1:882962/1‑100 (MQ=255)
                                      aGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCgag                                                            >  2:1329235/1‑98 (MQ=255)
                                        tGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAAc                                                          <  2:630481/100‑1 (MQ=255)
                                            ggCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTg                                                      >  2:2222164/1‑100 (MQ=255)
                                             gCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGc                                                     <  1:807603/100‑1 (MQ=255)
                                               tccgatctgACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCgg                                                   <  2:1376265/92‑1 (MQ=255)
                                               ggTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCgg                                                   <  1:745528/100‑1 (MQ=255)
                                               ggTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCgg                                                   <  2:2378521/100‑1 (MQ=255)
                                                gTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCgg                                                   <  2:1410079/99‑1 (MQ=255)
                                                 tGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGt                                                  >  1:1954909/1‑99 (MQ=255)
                                                 tGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTa                                                 >  1:2412872/1‑100 (MQ=255)
                                                    aGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCa                                              <  1:1685190/100‑1 (MQ=255)
                                                     gCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAg                                             <  1:1259635/100‑1 (MQ=255)
                                                      cGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAgg                                            >  2:2179190/1‑100 (MQ=255)
                                                       gACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGa                                           >  2:2071740/1‑100 (MQ=255)
                                                       gACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCGGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCggagatcgga                                           >  1:1376257/1‑92 (MQ=255)
                                                         cAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACg                                         >  2:164202/1‑100 (MQ=255)
                                                            gTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGACCCGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTg                                      >  1:1602824/1‑100 (MQ=255)
                                                              tCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGcc                                    >  2:1044143/1‑100 (MQ=255)
                                                                cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAg                                  >  1:1844856/1‑100 (MQ=255)
                                                                cTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAg                                  >  2:2360976/1‑100 (MQ=255)
                                                                 tCAGCGAATTGCCGCTCGGCTGGAATCTGGATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACACTGcc                                    <  1:329713/97‑1 (MQ=255)
                                                                 tCAGCGAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAAATCTGCGGTAGCAGGACGCTGCCAGc                                 >  1:2442533/1‑100 (MQ=255)
                                                                   aGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGa                               >  2:1402023/1‑100 (MQ=255)
                                                                   aGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGa                               >  2:129203/1‑100 (MQ=255)
                                                                   aGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGa                               >  2:964014/1‑100 (MQ=255)
                                                                    gCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGaa                              >  1:2395655/1‑100 (MQ=255)
                                                                    gCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTACCAGCGaa                              >  1:2367388/1‑100 (MQ=255)
                                                                      gAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACt                            <  2:213757/100‑1 (MQ=255)
                                                                      gAATTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACt                            >  2:2532171/1‑100 (MQ=255)
                                                                        aTTGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCg                          <  2:2215812/100‑1 (MQ=255)
                                                                         ttGCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGc                         <  2:1552954/100‑1 (MQ=255)
                                                                         ttGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGc                         <  1:724854/100‑1 (MQ=255)
                                                                           gCCGCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAg                       >  1:992208/1‑100 (MQ=255)
                                                                           gCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAg                       >  1:440662/1‑100 (MQ=255)
                                                                             cGCTCGCCTTGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGtt                     <  2:2584070/100‑1 (MQ=255)
                                                                             cGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGtt                     >  2:1691818/1‑100 (MQ=255)
                                                                              gCTCGCCTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGttt                    >  1:2230886/1‑100 (MQ=255)
                                                                              gCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGttt                    >  2:455395/1‑100 (MQ=255)
                                                                               cTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTg                   <  2:497783/100‑1 (MQ=255)
                                                                                 cGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGc                  <  2:577788/99‑1 (MQ=255)
                                                                                 cGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCa                 >  2:1734354/1‑100 (MQ=255)
                                                                                   ccTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAg               >  1:2474739/1‑100 (MQ=255)
                                                                                    cTGGAATCTGTATGCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGt              >  2:2005251/1‑100 (MQ=255)
                                                                                      ggAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa            >  2:1465095/1‑100 (MQ=255)
                                                                                      ggAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGa            >  2:637369/1‑100 (MQ=255)
                                                                                         aTCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACgg         >  2:8221/1‑100 (MQ=255)
                                                                                           cTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGta       >  1:569123/1‑100 (MQ=255)
                                                                                               ataCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGtatata   <  2:225294/100‑1 (MQ=255)
                                                                                                taCCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATaa  >  2:611421/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GCCAATCAGCGAGCTTTTGCGGTCGCGATTGAACAGGCAGTGCAGGCGGTGCAGCGACAAGTTCCTCAGCGAATTGCCGCTCGCCTGGAATCTGTATACCCAGCTGGTTTTGAGCTGCTGGACGGTGGCAAAAGCGGAACTCTGCGGTAGCAGGACGCTGCCAGCGAACTCGCAGTTTGCAAGTGACGGTATATAA  >  NC_000913/1148542‑1148737

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