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

TGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAG  >  NC_000913/1129250‑1129542
                                                                                                                                                |                                                                                                                                                    
tGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgca                                                                                                                                                  <  1:495747/149‑1 (MQ=255)
  ggTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAgg                                                                                                                                                >  4:62913/1‑149 (MQ=255)
     aCCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCg                                                                                                                                             <  2:39306/149‑1 (MQ=255)
       cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGga                                                                                                                                           <  4:87110/149‑1 (MQ=255)
         tGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGgaat                                                                                                                                         <  2:159920/149‑3 (MQ=255)
              tGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATc                                                                                                                                    >  3:334491/1‑149 (MQ=255)
                gCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCtt                                                                                                                                  >  2:407533/1‑149 (MQ=255)
                  gTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCa                                                                                                                                >  3:722206/1‑149 (MQ=255)
                   tttACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAg                                                                                                                               <  3:62913/149‑1 (MQ=255)
                        tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTc                                                                                                                                 <  3:647029/142‑1 (MQ=255)
                        tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTc                                                                                                                                 >  4:647029/1‑142 (MQ=255)
                           gTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAAt                                                                                                                       >  1:25430/1‑149 (MQ=255)
                                aTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                  <  1:228168/149‑1 (MQ=255)
                                aTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                  <  3:375489/149‑1 (MQ=255)
                                 tGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttttt                                                                                                                 >  3:473410/1‑149 (MQ=255)
                                 tGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttttt                                                                                                                 <  3:240118/149‑1 (MQ=255)
                                   gCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAc                                                                                                               >  3:228875/1‑149 (MQ=255)
                                               gCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                   >  1:717452/1‑149 (MQ=255)
                                               gCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATgg                                                                                                    >  4:181558/1‑148 (MQ=255)
                                               gCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATgg                                                                                                    <  3:181558/148‑1 (MQ=255)
                                                 ggggATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgtt                                                                                                 >  4:594616/1‑149 (MQ=255)
                                                    gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttgtt                                                                                              <  3:594616/149‑1 (MQ=255)
                                                     aTTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCtt                                                                                                                                  <  4:584132/112‑1 (MQ=255)
                                                     aTTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCtt                                                                                                                                  >  3:584132/1‑112 (MQ=255)
                                                       tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGGCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                           <  4:3955/149‑1 (MQ=255)
                                                       tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                           <  4:228875/149‑1 (MQ=255)
                                                       tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                           <  1:407533/149‑1 (MQ=255)
                                                       tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTCTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGCTGGGTTGTTGAg                                                                                           <  4:475749/149‑1 (MQ=255)
                                                          cgcgTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTc                                                                                        <  1:559148/149‑1 (MQ=255)
                                                              tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAc                                                                                    <  3:716460/149‑1 (MQ=255)
                                                                 ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGa                                                                                                       >  3:591004/1‑127 (MQ=255)
                                                                 ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGa                                                                                                       <  4:591004/127‑1 (MQ=255)
                                                                 ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                 >  2:196104/1‑149 (MQ=255)
                                                                           tGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCata                                                                       <  1:303688/149‑1 (MQ=255)
                                                                              cGGTACTGGGCTTCAAAGTGAAAGAATTTGCCCGCCGGACGGTGTAACAATGCGGTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAg                                                                    >  4:523911/1‑149 (MQ=255)
                                                                                   cTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCg                                                               <  4:722206/149‑1 (MQ=255)
                                                                                         ttCAAAGTTAAAGAATTTGCGCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                         >  1:616270/1‑149 (MQ=255)
                                                                                         ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                           >  3:249612/1‑115 (MQ=255)
                                                                                         ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                           <  4:249612/115‑1 (MQ=255)
                                                                                         ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                         >  4:17079/1‑149 (MQ=255)
                                                                                         ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                         <  4:155735/149‑1 (MQ=255)
                                                                                            aaaGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt                                                      >  1:435007/1‑149 (MQ=255)
                                                                                                ttAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCgg                                                  >  3:213307/1‑149 (MQ=255)
                                                                                                tgaaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCgg                                                  <  4:138540/147‑1 (MQ=255)
                                                                                                  aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                                >  2:642871/1‑149 (MQ=255)
                                                                                                  aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                                >  1:97361/1‑149 (MQ=255)
                                                                                                  aaaGAATTTGCCCGCCGGACGGAGTAACAATGCATTTCGGCCTGCAGTGCAGGCAGGAGCTAATCTTCAGATCGAAATCTTTTTACCGACGCGATGGGCTGGTGAGGTCTGACCGTTCGCCCCATGTAGCGTCAGTTCGTGATGCGGtt                                                >  2:2415/1‑149 (MQ=255)
                                                                                                    aGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                              <  3:523911/149‑1 (MQ=255)
                                                                                                     gAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCa                                             >  3:171053/1‑149 (MQ=255)
                                                                                                      aaTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCaa                                            <  3:398688/149‑1 (MQ=255)
                                                                                                      aaTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCaa                                            <  4:334491/149‑1 (MQ=255)
                                                                                                       aTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAAc                                           <  4:88598/149‑1 (MQ=255)
                                                                                                             ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTcc                                     >  3:316306/1‑149 (MQ=255)
                                                                                                             ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTcc                                     >  3:238721/1‑149 (MQ=255)
                                                                                                             ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTcc                                     >  2:83032/1‑149 (MQ=255)
                                                                                                              cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCa                                    <  3:227507/149‑1 (MQ=255)
                                                                                                                          tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAt                        <  4:749879/149‑1 (MQ=255)
                                                                                                                          tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAt                        >  1:493815/1‑149 (MQ=255)
                                                                                                                           aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAtt                       <  2:661805/149‑1 (MQ=255)
                                                                                                                              aaTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcg                    >  3:476852/1‑149 (MQ=255)
                                                                                                                               aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcgc                   >  1:158718/1‑149 (MQ=255)
                                                                                                                               aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcgc                   >  1:203282/1‑149 (MQ=255)
                                                                                                                               aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcgc                   >  3:143457/1‑149 (MQ=255)
                                                                                                                                   ataccGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTGaa               >  3:504363/4‑149 (MQ=255)
                                                                                                                                    ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                         <  2:94151/138‑1 (MQ=255)
                                                                                                                                    ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                         >  1:94151/1‑138 (MQ=255)
                                                                                                                                    ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt              <  4:213307/149‑1 (MQ=255)
                                                                                                                                    ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt              >  1:329820/1‑149 (MQ=255)
                                                                                                                                    ctCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                   <  4:603597/63‑1 (MQ=255)
                                                                                                                                    ctCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                   >  3:603597/2‑64 (MQ=255)
                                                                                                                                     tCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATc             <  4:140460/149‑1 (MQ=255)
                                                                                                                                     tCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATc             >  2:131211/1‑149 (MQ=255)
                                                                                                                                        ggCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGt          >  4:17304/1‑149 (MQ=255)
                                                                                                                                          ccTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTcc        <  4:316306/149‑1 (MQ=255)
                                                                                                                                             gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTc     <  3:17079/149‑1 (MQ=255)
                                                                                                                                              cagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc    <  4:473410/149‑1 (MQ=255)
                                                                                                                                              cagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTcc    <  2:126065/149‑1 (MQ=255)
                                                                                                                                               agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCa   <  2:435007/149‑1 (MQ=255)
                                                                                                                                               agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCa   <  3:76234/149‑1 (MQ=255)
                                                                                                                                                gtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAg  <  3:371470/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCAG  >  NC_000913/1129250‑1129542

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