Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,129,394 G→C 6.8% 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.8% 254.1 / 9.5 73intergenic (+20/+20)murJ/flgNputative lipid II flippase/export chaperone for FlgK and FlgL
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base C (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

ATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCA  >  NC_000913/1129248‑1129541
                                                                                                                                                  |                                                                                                                                                   
atcgggTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCCCGTACTTCGCTGCACTGGCGGAACTGGGCTTCAAAGGTAAAGAATTTGCCCGTCGGACGGAGTAACCACGCATTCCGGCCtgcagtg                                                                                                                                                   >  2:477915/4‑149 (MQ=255)
aTTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtg                                                                                                                                                   <  4:300527/149‑1 (MQ=255)
aTTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCagcagtg                                                                                                                                                   >  4:528042/1‑149 (MQ=255)
 ttGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCGAAGTTTAAGAATTTGCCCGCCGGGCGGTGTAACAATGCATTCCGGCCtgcagtgc                                                                                                                                                  >  2:286324/1‑149 (MQ=255)
   gggTACCATGCCCTGGCGTTTATTGCGTTTAATGGCGTTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTACAGAATTTTCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgcag                                                                                                                                                <  4:848524/149‑1 (MQ=255)
         cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                          >  3:655175/1‑149 (MQ=255)
          aTGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGgag                                                                                                                                         >  3:146400/1‑149 (MQ=255)
          aTGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGgag                                                                                                                                         >  2:757245/1‑149 (MQ=255)
          aTGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGgaa                                                                                                                                         <  1:491393/149‑2 (MQ=255)
               cTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAAt                                                                                                                                    >  4:575043/1‑149 (MQ=255)
                tGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATc                                                                                                                                   <  1:329415/149‑1 (MQ=255)
                tGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATc                                                                                                                                   <  2:245507/149‑1 (MQ=255)
                    gTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCa                                                                                                                               <  2:63565/149‑1 (MQ=255)
                       tACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGAt                                                                                                                            >  3:262876/1‑149 (MQ=255)
                         cTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCg                                                                                                                          >  4:537372/1‑149 (MQ=255)
                          tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGa                                                                                                                         >  3:646584/1‑149 (MQ=255)
                               ttAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCt                                                                                                                    >  1:813510/1‑149 (MQ=255)
                                  aTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                 <  4:676884/149‑1 (MQ=255)
                                  aTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                 <  4:646584/149‑1 (MQ=255)
                                   tGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttttt                                                                                                                <  2:692925/149‑1 (MQ=255)
                                   tGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttttt                                                                                                                <  2:835692/149‑1 (MQ=255)
                                         tCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGAGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAccgcc                                                                                                          <  4:581430/149‑1 (MQ=255)
                                                ggCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATgg                                                                                                   <  3:11913/149‑1 (MQ=255)
                                                  cGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgt                                                                                                 <  1:757245/149‑1 (MQ=255)
                                                    gggATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttg                                                                                               <  4:480421/149‑1 (MQ=255)
                                                               gTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGa                                                                                    <  1:640354/149‑1 (MQ=255)
                                                                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGa                                                                                                                             <  3:725089/107‑1 (MQ=255)
                                                                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGa                                                                                                                             <  1:596455/107‑1 (MQ=255)
                                                                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGa                                                                                                                             >  4:725089/1‑107 (MQ=255)
                                                                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGa                                                                                                                             >  2:596455/1‑107 (MQ=255)
                                                                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAc                                                                                   >  4:264401/1‑149 (MQ=255)
                                                                   ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  1:180747/1‑149 (MQ=255)
                                                                   ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  3:524482/1‑149 (MQ=255)
                                                                   ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  1:445772/1‑149 (MQ=255)
                                                                   ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCGCCGGACGGTGTAACCATGCATTCCGGCCTGCAGAGCAGGCCGGAGATAATCTTTAGAACGAAATTTTTTTTCCGCCGCGGTGGGGTGTTTAGGTCTGGCCGtt                                                                               >  2:648096/1‑149 (MQ=255)
                                                                     cGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTc                                                                              <  1:708995/149‑1 (MQ=255)
                                                                        tGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGcc                                                                           <  4:721361/149‑1 (MQ=255)
                                                                             tGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCata                                                                      >  2:101843/1‑149 (MQ=255)
                                                                                 ggTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGc                                                                  <  1:637521/149‑1 (MQ=255)
                                                                                      tGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCgg                                                             >  1:522867/1‑149 (MQ=255)
                                                                                      tGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCgg                                                             >  3:506754/1‑149 (MQ=255)
                                                                                       gggCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGt                                                            <  3:114845/149‑1 (MQ=255)
                                                                                           ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                        >  2:370241/1‑149 (MQ=255)
                                                                                           ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                        <  1:101843/149‑1 (MQ=255)
                                                                                             cAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGa                                                      >  4:154928/1‑149 (MQ=255)
                                                                                              aaaGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt                                                     >  4:298180/1‑149 (MQ=255)
                                                                                               aaGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATg                                                    >  1:445351/1‑149 (MQ=255)
                                                                                                 gTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCg                                                  >  3:635286/1‑149 (MQ=255)
                                                                                                  ttAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCgg                                                 >  2:751818/1‑149 (MQ=255)
                                                                                                    aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                               >  1:581931/1‑149 (MQ=255)
                                                                                                    aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                               >  2:370545/1‑149 (MQ=255)
                                                                                                    aaaGAATTTGCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttt                                              <  1:648096/149‑1 (MQ=255)
                                                                                                     aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttt                                              <  1:790576/149‑1 (MQ=255)
                                                                                                     aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCg                                                                                                                                            <  2:298049/55‑1 (MQ=255)
                                                                                                     aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCg                                                                                                                                            >  1:298049/1‑55 (MQ=255)
                                                                                                     aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCAGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttt                                              <  3:528042/149‑1 (MQ=255)
                                                                                                      aGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                             >  4:30952/1‑149 (MQ=255)
                                                                                                        aaTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTCTCaa                                           <  4:506754/149‑1 (MQ=255)
                                                                                                         aTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAAc                                          >  2:662242/1‑149 (MQ=255)
                                                                                                            tGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACAtt                                       <  1:494920/149‑1 (MQ=255)
                                                                                                             gCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACAttt                                      <  3:432158/149‑1 (MQ=255)
                                                                                                               ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGa                                                                                           <  1:812297/94‑1 (MQ=255)
                                                                                                               ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGa                                                                                           >  2:812297/1‑94 (MQ=255)
                                                                                                                         gtgtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGtt                          <  2:802183/149‑1 (MQ=255)
                                                                                                                          tgtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTg                         <  4:797499/149‑1 (MQ=255)
                                                                                                                           gtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                        <  3:861130/149‑1 (MQ=255)
                                                                                                                           gtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                        <  2:813510/149‑1 (MQ=255)
                                                                                                                             aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAtt                      >  1:319714/1‑149 (MQ=255)
                                                                                                                             aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCAATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAtt                      >  3:86784/1‑149 (MQ=255)
                                                                                                                                aaTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcg                   >  4:498442/1‑149 (MQ=255)
                                                                                                                                  tGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCt                 >  3:340130/1‑149 (MQ=255)
                                                                                                                                   gCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTc                >  1:401125/1‑149 (MQ=255)
                                                                                                                                     aTTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGCCCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCaa              <  3:129568/149‑1 (MQ=255)
                                                                                                                                      ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt             >  4:2997/1‑149 (MQ=255)
                                                                                                                                      ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt             <  4:655175/149‑1 (MQ=255)
                                                                                                                                       tCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTcc                                                         <  4:330881/104‑1 (MQ=255)
                                                                                                                                       tCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTcc                                                         >  3:330881/1‑104 (MQ=255)
                                                                                                                                            ccTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTcc       <  1:751818/149‑1 (MQ=255)
                                                                                                                                                cagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc   >  1:822423/1‑149 (MQ=255)
                                                                                                                                                cagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc   >  1:73662/1‑149 (MQ=255)
                                                                                                                                                cagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc   <  1:607456/149‑1 (MQ=255)
                                                                                                                                                 agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCa  <  1:658134/149‑1 (MQ=255)
                                                                                                                                                 agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCa  <  3:83055/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCA  >  NC_000913/1129248‑1129541

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