Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,129,394 G→C 5.5% 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 5.5% 585.4 / 19.6 165intergenic (+20/+20)murJ/flgNputative lipid II flippase/export chaperone for FlgK and FlgL
Reads supporting (aligned to +/- strand):  ref base G (74/82);  new base C (4/5);  total (78/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 178 total aligned reads displayed.

TCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCA  >  NC_000913/1129246‑1129541
                                                                                                                                                    |                                                                                                                                                   
  aTTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtg                                                                                                                                                   >  2:650160/1‑149 (MQ=255)
   ttGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgc                                                                                                                                                  >  2:931811/1‑149 (MQ=255)
     gggTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgcag                                                                                                                                                >  2:629240/1‑149 (MQ=255)
       gTAACATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGc                                                                                                                                              >  4:1581092/1‑149 (MQ=255)
         aCCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAAGAATGCATTCCGGCCTGCAGTGCAGGCCg                                                                                                                                            <  2:351405/149‑1 (MQ=255)
           cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                          >  2:154681/1‑149 (MQ=255)
           cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                          <  1:154681/149‑1 (MQ=255)
           cATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                          <  1:675002/149‑1 (MQ=255)
            aTGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCCCTGCAGGCCGgaa                                                                                                                                         <  1:1415202/149‑2 (MQ=255)
             tGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGgaga                                                                                                                                        <  1:1724637/149‑1 (MQ=255)
                  tGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATc                                                                                                                                   >  1:39476/1‑149 (MQ=255)
                      gTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCa                                                                                                                               >  3:1167851/1‑149 (MQ=255)
                         tACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGAt                                                                                                                            <  1:650160/149‑1 (MQ=255)
                            tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGa                                                                                                                         >  4:1510913/1‑149 (MQ=255)
                            tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGa                                                                                                                         >  3:1431262/1‑149 (MQ=255)
                            tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGa                                                                                                                         >  3:496213/1‑149 (MQ=255)
                            tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGa                                                                                                                         >  1:1623941/1‑149 (MQ=255)
                            tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGa                                                                                                                         >  4:201107/1‑149 (MQ=255)
                            tGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGga                                                                                                                                          <  2:1215109/132‑1 (MQ=255)
                             gCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCCGCAGTGCAGGCCGGAGATAATCTTCAGATCGaa                                                                                                                        <  3:1093239/149‑1 (MQ=255)
                               gTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAAt                                                                                                                      <  3:73364/149‑1 (MQ=255)
                               gTTTAATGACGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAAGTTGCCCGCCGGGCGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATTATATTCAGATCGAAAt                                                                                                                      >  4:494645/1‑149 (MQ=255)
                                   aaTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCttt                                                                                                                  <  2:409019/149‑1 (MQ=255)
                                    aTGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                 <  3:1049429/149‑1 (MQ=255)
                                          gTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGCGCAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAccgc                                                                                                           >  3:988030/1‑149 (MQ=255)
                                             gTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACcgccgc                                                                                                        <  2:358812/149‑1 (MQ=255)
                                                        gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                  >  1:969512/1‑144 (MQ=255)
                                                        gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                                                                  <  2:969512/144‑1 (MQ=255)
                                                        gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttgtt                                                                                             <  3:1320799/149‑1 (MQ=255)
                                                          ttGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGa                                                                                           >  2:1632767/1‑149 (MQ=255)
                                                          ttGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAAAAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGa                                                                                           <  4:1402452/149‑1 (MQ=255)
                                                                 gTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                          >  4:487128/1‑95 (MQ=255)
                                                                 gTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGa                                                                                    <  4:1364011/149‑1 (MQ=255)
                                                                 gTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGa                                                                                    <  4:1567960/149‑1 (MQ=255)
                                                                  tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCGTCAGACCGAAATCTCTTTACCGCCGCTAAGGGTTGTAGAGGACTggt                                                                                   <  2:1070231/149‑3 (MQ=255)
                                                                    cTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAACg                                                                                 >  3:1368320/1‑149 (MQ=255)
                                                                     ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  1:854696/1‑149 (MQ=255)
                                                                     ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  4:1605394/1‑149 (MQ=255)
                                                                     ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  4:414386/1‑149 (MQ=255)
                                                                     ttCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGt                                                                                >  4:1137662/1‑149 (MQ=255)
                                                                       cGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTc                                                                              <  2:500722/149‑1 (MQ=255)
                                                                         ctccACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGc                                                                            <  4:697075/146‑1 (MQ=255)
                                                                         cTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                          >  4:702366/1‑135 (MQ=255)
                                                                         cTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAg                                                                                          <  3:702366/135‑1 (MQ=255)
                                                                         cTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGc                                                                            >  1:1367587/1‑149 (MQ=255)
                                                                              cTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCat                                                                       <  3:1581092/149‑1 (MQ=255)
                                                                               tGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCata                                                                      >  1:196733/1‑149 (MQ=255)
                                                                                     tACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGgaat                                                                                                                                        >  1:1453793/1‑75 (MQ=255)
                                                                                      aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc                                                               <  2:1371293/149‑1 (MQ=255)
                                                                                      aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc                                                               <  1:1207790/149‑1 (MQ=255)
                                                                                      aCTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTGTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc                                                               <  4:1443485/149‑1 (MQ=255)
                                                                                          ggCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGtt                                                           >  3:838933/1‑149 (MQ=255)
                                                                                           gCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTc                                                          <  4:1588878/149‑1 (MQ=255)
                                                                                             ttCAAAGTTAAAGATTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                                                                          >  3:218384/1‑67 (MQ=255)
                                                                                             ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCAtata                                                                    >  3:1286291/1‑137 (MQ=255)
                                                                                             ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCAtata                                                                    <  4:1286291/137‑1 (MQ=255)
                                                                                             ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                        <  2:1129677/149‑1 (MQ=255)
                                                                                             ttCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt                                                        >  3:949695/1‑149 (MQ=255)
                                                                                                aaaGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt                                                     <  1:153543/149‑1 (MQ=255)
                                                                                                aaaGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt                                                     >  4:918783/1‑149 (MQ=255)
                                                                                                  aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc                                                               >  2:1683296/1‑137 (MQ=255)
                                                                                                   gTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGTGGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGAAGGTTTGTTGAGGTCTGCCCGTTCCCCCGATATAGCGCCGGTTCCTgttaga                                                  >  1:1236460/1‑144 (MQ=255)
                                                                                                    ttAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCgg                                                 >  2:1570288/1‑149 (MQ=255)
                                                                                                      aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGtt                                               >  3:958907/1‑149 (MQ=255)
                                                                                                       aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGttt                                              >  3:1291109/1‑149 (MQ=255)
                                                                                                         gAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCa                                            >  1:1751674/1‑149 (MQ=255)
                                                                                                               gCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACAttt                                      <  2:152634/149‑1 (MQ=255)
                                                                                                                 ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTcc                                    >  1:907820/1‑149 (MQ=255)
                                                                                                                  cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCa                                   <  2:726447/149‑1 (MQ=255)
                                                                                                                           gtgtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGtt                          <  3:656114/149‑1 (MQ=255)
                                                                                                                            tgtAACAATGCATTCCGTCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                 >  1:136233/1‑61 (MQ=255)
                                                                                                                            tgtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                                                                 <  2:136233/61‑1 (MQ=255)
                                                                                                                             gtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGa                        >  3:123322/1‑149 (MQ=255)
                                                                                                                              tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAt                       <  3:1137662/149‑1 (MQ=255)
                                                                                                                              tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAt                       >  3:1121818/1‑149 (MQ=255)
                                                                                                                               aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                             >  1:624856/1‑126 (MQ=255)
                                                                                                                               aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTc                                             <  2:624856/126‑1 (MQ=255)
                                                                                                                               aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAtt                      <  4:648270/149‑1 (MQ=255)
                                                                                                                               aaCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAtt                      <  2:909561/149‑1 (MQ=255)
                                                                                                                                aCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTg                     >  2:1012564/1‑149 (MQ=255)
                                                                                                                                 cAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGTGCCTGTTGATTgc                    >  3:915531/1‑149 (MQ=255)
                                                                                                                                 cAAGGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGGTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCCGTTGATTgc                    >  4:843057/1‑149 (MQ=255)
                                                                                                                                  aaTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcg                   >  4:1640815/1‑149 (MQ=255)
                                                                                                                                   aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcgc                  <  4:958907/149‑1 (MQ=255)
                                                                                                                                   aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcgc                  <  3:1088305/149‑1 (MQ=255)
                                                                                                                                   aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcg                   <  3:668437/148‑1 (MQ=255)
                                                                                                                                   aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTgcg                   >  4:668437/1‑148 (MQ=255)
                                                                                                                                      atctCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTc                                                          <  4:1361876/103‑1 (MQ=255)
                                                                                                                                        ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt             >  3:561197/1‑149 (MQ=255)
                                                                                                                                        ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt             <  3:414386/149‑1 (MQ=255)
                                                                                                                                        ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAAt             >  3:1687208/1‑149 (MQ=255)
                                                                                                                                         tCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATc            >  1:1111840/1‑149 (MQ=255)
                                                                                                                                          ccGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCt           <  4:595313/149‑1 (MQ=255)
                                                                                                                                          ccGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCt           >  1:1048757/1‑149 (MQ=255)
                                                                                                                                            ggCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGt         <  3:771128/149‑1 (MQ=255)
                                                                                                                                             gCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGAtt                      >  2:1386097/1‑135 (MQ=255)
                                                                                                                                              ccTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTcc       <  1:1644925/149‑1 (MQ=255)
                                                                                                                                                 gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTc    <  4:1431262/149‑1 (MQ=255)
                                                                                                                                                 gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTc    <  2:1751674/149‑1 (MQ=255)
                                                                                                                                                   agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCa  <  2:1515189/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCATTGGGTACCATGCCCTGGCGTTTACTGCGTTTAATGGCGGTCGTGCTGGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGCGGTTTCAACATTTCCAGCGCCTGTTGATTGCGCTCAATCTGTCCTTCCA  >  NC_000913/1129246‑1129541

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