Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,129,394 G→C 6.1% intergenic (+20/+20) murJ → / ← flgN putative lipid II flippase/export chaperone for FlgK and FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,129,3940GC 6.1% 305.7 / 9.6 99intergenic (+20/+20)murJ/flgNputative lipid II flippase/export chaperone for FlgK and FlgL
Reads supporting (aligned to +/- strand):  ref base G (51/42);  new base C (3/3);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

Only 100 of 121 total aligned reads displayed.

GTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAGT  >  NC_000913/1129245‑1129543
                                                                                                                                                      |                                                                                                                                                     
gTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑tgcag                                                                                                                                                       >  1:740604/1‑150 (MQ=255)
  cATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑tgcagtg                                                                                                                                                     >  1:340933/1‑150 (MQ=255)
  cATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGACGGACGGTGTAAGAATGCATTCAGGCC‑tgcagtg                                                                                                                                                     >  2:1570871/1‑150 (MQ=255)
  cATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAAAGGCGGTCGTACTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTAAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑tgcaggg                                                                                                                                                     >  1:1506043/1‑148 (MQ=255)
   aTTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑tgcagtgc                                                                                                                                                    >  1:1419121/1‑150 (MQ=255)
   aTTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGGGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑tgcagtgc                                                                                                                                                    <  2:1673775/150‑1 (MQ=255)
     tGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑tgcagtgcag                                                                                                                                                  >  1:560985/1‑150 (MQ=255)
     tGGGTACCATGCCCTGGCGTTTACTGCGTGTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTACGCGGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGCGTAACAATGCATTCCGGCC‑tgcagtgcag                                                                                                                                                  <  2:1140865/150‑1 (MQ=255)
        gTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCACTGCAGGcc                                                                                                                                               <  1:921135/150‑1 (MQ=255)
          aCCATGCCCTGGCGGTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCACTGCAGGCCgg                                                                                                                                             <  2:1109839/150‑1 (MQ=255)
          aCCATGCCCGGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCgg                                                                                                                                             <  1:624638/150‑1 (MQ=255)
            cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGTGGATTGCCTCGTACTTAGATTCACAGGCGGTACTGGGATTCAAAGTTAAAGAATTTGCCCGCCGGACGGAGTAACAATGCATTCCTGCC‑TGCAGAGCAGGCCGgag                                                                                                                                           >  2:1686484/1‑150 (MQ=255)
              tGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCACTGCAGGCCGgaatg                                                                                                                                         <  2:161731/150‑4 (MQ=255)
               gCCCTGGCGTTTACTGCGTTTAATGGCGTTCGTGCTGGCGGGGATTTCAGCGTACTTCGCTGCACTGGCGGTACTTTTCTTCAAAGTTAAAGAATTTGCGCGGCGTACGGTGTAACATTGCATTCCGGCA‑TTCAGTGCAGGCCGGAGATa                                                                                                                                        <  2:793264/150‑1 (MQ=255)
                  cTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCTGCC‑TGCAGTGCAGGCCGGAGATAATc                                                                                                                                     >  2:1389516/1‑150 (MQ=255)
                  cTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATc                                                                                                                                     >  2:1291481/1‑150 (MQ=255)
                    ggCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCtt                                                                                                                                   >  1:718831/1‑150 (MQ=255)
                     gCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTc                                                                                                                                  <  1:1600374/150‑1 (MQ=255)
                     gCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTc                                                                                                                                  <  2:1580152/150‑1 (MQ=255)
                      cGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGTGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCa                                                                                                                                 <  1:841428/150‑1 (MQ=255)
                          tACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATc                                                                                                                             >  1:1125981/1‑150 (MQ=255)
                              gCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGaaa                                                                                                                         >  1:1062116/1‑150 (MQ=255)
                                   tAATGGCGGTCGTTCTGGCGGTGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttt                                                                                                                    <  2:548372/150‑1 (MQ=255)
                                   tAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCCCTGGCGGTACTGGGCTTCAAAGTTAAAGACTTTGCCCGACGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGAAGGCCGGAGATAATCTTCAGATAGAAATCttt                                                                                                                    >  1:97612/1‑150 (MQ=255)
                                    aaTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGgaga                                                                                                                                          <  1:1630249/127‑1 (MQ=255)
                                    aaTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGgaga                                                                                                                                          >  2:1630249/1‑127 (MQ=255)
                                       ggCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAc                                                                                                                >  1:797607/1‑150 (MQ=255)
                                                gCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGAt                                                                                                       >  2:562800/1‑150 (MQ=255)
                                                gCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGAt                                                                                                       <  1:1077509/150‑1 (MQ=255)
                                                   ggCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                    >  1:323287/1‑150 (MQ=255)
                                                    gCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGTATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgt                                                                                                   <  1:887502/150‑1 (MQ=255)
                                                        ggATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttgtt                                                                                               >  2:1658383/1‑150 (MQ=255)
                                                         gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGttgttg                                                                                              <  2:340933/150‑1 (MQ=255)
                                                          aTTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAAGAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGttgttgt                                                                                             >  1:177291/1‑149 (MQ=255)
                                                          aTTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGa                                                                                             >  1:574110/1‑150 (MQ=255)
                                                            tGCCGCGTACTTCGCTGCACTGGCTTTACTGGGATTCAAAGTTAAATAATTTGCCCGCCGGACGGTGTAACATTGCTTTCCGTCC‑TGCAGTGCAGGCCGGAGTTAATCTTCAGATCGATATCTTTTTACCGCCGCTATGTGTTGTTGAgg                                                                                           <  1:807537/150‑1 (MQ=255)
                                                             gCCGCGTTCTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGt                                                                                          >  2:1363970/1‑150 (MQ=255)
                                                             gCCGCGTACTTCGCTGCACTGGCGGAACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGt                                                                                          <  2:15174/150‑1 (MQ=255)
                                                                  gTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGGGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAc                                                                                     <  1:1224973/150‑1 (MQ=255)
                                                                   tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAcc                                                                                    <  2:1269341/150‑1 (MQ=255)
                                                                     cTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                  >  1:1260948/1‑150 (MQ=255)
                                                                      ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAgg                                                                                           <  1:1643646/140‑1 (MQ=255)
                                                                         gCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGc                                                                              <  2:578302/150‑1 (MQ=255)
                                                                          cTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGcc                                                                             >  2:578863/1‑150 (MQ=255)
                                                                            gcacgcgCGGTACTGCGCCTCCTCGTTAATGAATTTGCCCGCTGGACGCTGTAACATTGCATTCCGGCC‑TGCAGCGCAGGCCGGAGTTAATCTTCATATCGAAATCTTTTTACCGCCGCGATGGGTGGTTGAGGTCTGACCGTTCGcccc                                                                           <  2:1075265/144‑1 (MQ=255)
                                                                            gCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCGTCAGCTCGAAATCTTTTTACCGCCGCGATGGGGTGTTGCGGTCTGACCGTTCGcccc                                                                           >  2:959637/1‑150 (MQ=255)
                                                                            gCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑GGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGcccc                                                                           <  2:1542668/150‑1 (MQ=255)
                                                                                    ggTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATATTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCg                                                                   >  2:1204158/1‑150 (MQ=255)
                                                                                        cTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCgg                                                               >  1:442782/1‑150 (MQ=255)
                                                                                        cTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCgg                                                               >  1:365845/1‑150 (MQ=255)
                                                                                             cTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                          <  1:381373/150‑1 (MQ=255)
                                                                                             cTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                          >  1:217389/1‑150 (MQ=255)
                                                                                              ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTg                                                         <  1:809838/150‑1 (MQ=255)
                                                                                              ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTg                                                         <  2:964426/150‑1 (MQ=255)
                                                                                              ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGACAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTg                                                         >  2:579656/1‑150 (MQ=255)
                                                                                                   aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACGCCCGCATTGGGTTGTTATGTGATCACCTTACGCCCAATATAGCTACGTTTCGTTATGCg                                                    >  2:1500035/1‑150 (MQ=255)
                                                                                                   aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTTGTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCg                                                    <  2:206119/150‑1 (MQ=255)
                                                                                                    gTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAAGCgg                                                   >  1:317284/1‑150 (MQ=255)
                                                                                                    gTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCGTTTTACCGCCGCGAGGGGTTGTTGAGGTCTGACCGTTCGCCCCATATTGCGTCGGTTCCTGATGCgg                                                   >  2:894946/1‑150 (MQ=255)
                                                                                                      tAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGGTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                                 >  2:1299623/1‑150 (MQ=255)
                                                                                                       aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttt                                                >  1:988302/1‑150 (MQ=255)
                                                                                                       aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttt                                                >  1:1256958/1‑150 (MQ=255)
                                                                                                        aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                               >  2:1473943/1‑150 (MQ=255)
                                                                                                        aaGAATTTGCCCGCCGGAAGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCATCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttta                                               >  1:1249199/1‑149 (MQ=255)
                                                                                                           aaTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAAc                                            <  1:53476/150‑1 (MQ=255)
                                                                                                             tttCCCCGCCGTAGGGTGTAA‑AATGCCTTCCGGCCGTCCAGTGCTGGCTGTTTTTAA‑CTTTAATTCCTAATTTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACAtt                                         <  1:1500035/150‑1 (MQ=255)
                                                                                                                gCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTc                                       <  2:1715004/150‑1 (MQ=255)
                                                                                                                gCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTc                                       >  2:325781/1‑150 (MQ=255)
                                                                                                                 cccGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTcc                                      >  1:1013448/1‑150 (MQ=255)
                                                                                                                       ggACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                               >  1:946799/1‑135 (MQ=255)
                                                                                                                       ggACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                               <  2:946799/135‑1 (MQ=255)
                                                                                                                         aCGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGATATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGCCCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTg                              >  2:1621600/1‑150 (MQ=255)
                                                                                                                             tgtAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                          >  1:913300/1‑150 (MQ=255)
                                                                                                                             tgtAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCAATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                          <  2:1419121/150‑1 (MQ=255)
                                                                                                                              gtAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTACAGCGCCTGTTGAt                         >  2:1179018/1‑150 (MQ=255)
                                                                                                                              gtAACAAGGCATTCCGGCC‑TGCAGTGCAGGCCTGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAt                         <  1:1516923/150‑1 (MQ=255)
                                                                                                                               aaaCAATGCATTCCGGCC‑TGCAGTTCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCTCCGCGATTGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACTTTTCCAGCGACTGTTGAtt                        <  1:707084/149‑1 (MQ=255)
                                                                                                                                 aCAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGAGCGAAATCTTTTTACCGCCGCGATGGGTTGGTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgc                      <  2:97612/150‑1 (MQ=255)
                                                                                                                                  cAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGATATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcg                     >  2:53016/1‑150 (MQ=255)
                                                                                                                                   aaTGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATTTATGGTCGGGTACTGATGCGGTTTAAACATTTCAAGCGGCTGTTGATTgcgc                    >  2:865192/1‑150 (MQ=255)
                                                                                                                                    aTGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCt                   <  2:574110/150‑1 (MQ=255)
                                                                                                                                      gCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCa                 >  2:907100/1‑150 (MQ=255)
                                                                                                                                        tctCCGGCC‑TGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCaaa               >  2:761378/3‑149 (MQ=255)
                                                                                                                                          tCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCt             >  2:1684178/1‑150 (MQ=255)
                                                                                                                                           ccGGCC‑TGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTg            >  1:1491135/1‑150 (MQ=255)
                                                                                                                                             ggCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGCTCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGCTGCGGTTTCAACATTTCCAGCGCCTGTTGAGTGCGCTCAATCTGTc          >  1:1620103/1‑150 (MQ=255)
                                                                                                                                             ggCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTc          <  1:372599/150‑1 (MQ=255)
                                                                                                                                               cc‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCt        >  1:1013314/1‑150 (MQ=255)
                                                                                                                                                   tcagtgcagGCCGGAGATAAGATTCAGATCGAAATCTTTTTACCGCCGCGATGGCTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc     <  2:282976/149‑1 (MQ=255)
                                                                                                                                                   gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc     >  2:1125990/1‑150 (MQ=255)
                                                                                                                                                   gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc     >  2:1289568/1‑150 (MQ=255)
                                                                                                                                                   gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGAGGGGTTGTTGAGGTCTGACCGTTCGCACCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc     >  2:1514671/1‑150 (MQ=255)
                                                                                                                                                   gcagtgcagGCCGGAGATAAACTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc     <  2:463012/150‑1 (MQ=255)
                                                                                                                                                     agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAg   <  2:365845/150‑1 (MQ=255)
                                                                                                                                                     agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAg   <  1:1291481/150‑1 (MQ=255)
                                                                                                                                                     agtgcagGCCGGAGATAATCTTCAGAGCGAAATCTTATTCCCGCCGCGATGGGGTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTGCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAg   <  1:856719/150‑1 (MQ=255)
                                                                                                                                                      gtgctgGGCGGAGATTATCTTCAGATCGAATTCTTTTTACCGCCGCGATGGGTTGTTGTGGTCTGTTCGTTCGCCCCATTTAGCGTCTGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAGt  <  1:1674156/150‑1 (MQ=255)
                                                                                                                                                      gtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGTGTTCAATCTGTCCTTCCAGt  >  2:461595/1‑150 (MQ=255)
                                                                                                                                                      gtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAGt  >  1:586866/1‑150 (MQ=255)
                                                                                                                                                      gtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAGt  >  2:1408429/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GTCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCC‑TGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAGT  >  NC_000913/1129245‑1129543

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