Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 975,700 G→T 5.3% D27Y (GAC→TAC) ‡ mukE → chromosome condensin MukBEF, MukE localization factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913975,7000GT 5.3% 260.9 / 4.4 75D27Y (GAC→TAC) ‡mukEchromosome condensin MukBEF, MukE localization factor
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base T (2/2);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

AAGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTA  >  NC_000913/975553‑975846
                                                                                                                                                   |                                                                                                                                                  
aaGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGa                                                                                                                                                   >  4:64773/1‑149 (MQ=255)
aaGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGa                                                                                                                                                   >  4:184500/1‑149 (MQ=255)
  gCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACa                                                                                                                                                 >  1:452681/1‑149 (MQ=255)
      aTTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAgcgc                                                                                                                                             >  2:538071/1‑149 (MQ=255)
       tttCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGcc                                                                                                                                            >  2:21966/1‑149 (MQ=255)
         tCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGTCCAGCGccgg                                                                                                                                          <  2:17754/149‑3 (MQ=255)
         tCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCtt                                                                                                                                          <  3:50148/149‑1 (MQ=255)
           aCCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                                                                                                                                        >  4:233341/1‑149 (MQ=255)
           aCCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTAc                                                                                                                                        >  3:382692/1‑149 (MQ=255)
                 cTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACATCGCCTTACGTTCAg                                                                                                                                  <  2:696883/149‑1 (MQ=255)
                 cTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAg                                                                                                                                  <  1:124716/149‑1 (MQ=255)
                  tGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAgg                                                                                                                                 >  1:35049/1‑149 (MQ=255)
                     cAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACg                                                                                                                              >  4:751349/1‑149 (MQ=255)
                      aGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                             >  4:378692/1‑149 (MQ=255)
                      aGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                             >  2:423390/1‑149 (MQ=255)
                       gCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGcc                                                                                                                            <  1:774773/149‑1 (MQ=255)
                           aaTGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCatat                                                                                                                        >  1:328057/1‑149 (MQ=255)
                              ggCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTAATCGACAGATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTTAGGACGCCATAATgg                                                                                                                     >  3:15848/1‑149 (MQ=255)
                               gCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGc                                                                                                                    <  1:453523/149‑1 (MQ=255)
                                     gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac                                                                                                              >  1:50817/1‑149 (MQ=255)
                                     gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac                                                                                                              >  1:341300/1‑149 (MQ=255)
                                     gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac                                                                                                              <  1:300933/149‑1 (MQ=255)
                                     gATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTcgac                                                                                                              >  3:565963/1‑149 (MQ=255)
                                           tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                             >  3:284005/1‑128 (MQ=255)
                                           tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGc                                                                                                                             <  4:284005/128‑1 (MQ=255)
                                           tGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTg                                                                                                        <  4:168249/149‑1 (MQ=255)
                                                    aGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAAc                                                                                                          <  1:588949/138‑1 (MQ=255)
                                                    aGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAAc                                                                                                          >  2:588949/1‑138 (MQ=255)
                                                    aGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCAt                                                                                               <  3:728083/149‑1 (MQ=255)
                                                     gCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATg                                                                                              >  4:11731/1‑149 (MQ=255)
                                                         aGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCAtt                                                                                          >  1:137837/1‑149 (MQ=255)
                                                          ggTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTc                                                                                         <  3:184500/149‑1 (MQ=255)
                                                          ggTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTc                                                                                         <  3:442845/149‑1 (MQ=255)
                                                              cAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGa                                                                                     <  1:423164/149‑1 (MQ=255)
                                                                 gcgcATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATgg                                                                                  >  3:91245/1‑149 (MQ=255)
                                                                    cATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGAtt                                                                               <  3:233341/149‑1 (MQ=255)
                                                                          aTCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCa                                                                                                >  4:564461/1‑126 (MQ=255)
                                                                          aTCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCa                                                                                                <  3:564461/126‑1 (MQ=255)
                                                                          aTCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAgg                                                                         >  1:272727/1‑149 (MQ=255)
                                                                             gACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAat                                                                      >  3:717914/1‑149 (MQ=255)
                                                                                 aaTATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCt                                                                  <  1:423390/149‑1 (MQ=255)
                                                                                 aaTATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCt                                                                  <  3:378692/149‑1 (MQ=255)
                                                                                  atatTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATAtt                                                                                                                       >  2:19757/1‑95 (MQ=255)
                                                                                  atatTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATAtt                                                                                                                       <  1:19757/95‑1 (MQ=255)
                                                                                    atTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGa                                                               >  4:173367/1‑149 (MQ=255)
                                                                                    atTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGa                                                               >  4:754102/1‑149 (MQ=255)
                                                                                      tGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAag                                                             >  3:7191/1‑149 (MQ=255)
                                                                                         aCAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGt                                                          >  1:551148/1‑149 (MQ=255)
                                                                                                 tGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcg                                                  >  2:710096/1‑149 (MQ=255)
                                                                                                  gCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTAcgcgc                                                 <  1:310636/149‑1 (MQ=255)
                                                                                                    cGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGt                                               <  3:257422/149‑1 (MQ=255)
                                                                                                             tGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGtt                                      >  3:660673/1‑149 (MQ=255)
                                                                                                             tGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGtt                                      <  2:213524/149‑1 (MQ=255)
                                                                                                               gcgcAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTTATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGa                                    <  2:228898/149‑1 (MQ=255)
                                                                                                                    ggCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTa                               <  4:382692/149‑1 (MQ=255)
                                                                                                                                tCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAg                   >  3:196876/1‑149 (MQ=255)
                                                                                                                                  cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCAt                                                                                               >  3:76191/1‑71 (MQ=255)
                                                                                                                                  cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCAt                                                                                               <  4:76191/71‑1 (MQ=255)
                                                                                                                                  cGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGaa                 >  2:583215/1‑149 (MQ=255)
                                                                                                                                   gTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAg                <  3:64773/149‑1 (MQ=255)
                                                                                                                                    ttATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAgg               >  4:724769/1‑149 (MQ=255)
                                                                                                                                     tATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAggg              >  1:367390/1‑149 (MQ=255)
                                                                                                                                        ttCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttc           <  2:732036/149‑1 (MQ=255)
                                                                                                                                         taaggCGCTGTCCAGCGCCCTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGttct          >  4:523832/4‑149 (MQ=255)
                                                                                                                                          aaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGaa                                                              >  3:65053/3‑96 (MQ=255)
                                                                                                                                          aaggCGCTGTCCAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGaa                                                              <  4:65053/94‑1 (MQ=255)
                                                                                                                                           cGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGtt                                      <  4:76613/119‑1 (MQ=255)
                                                                                                                                           cGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGtt                                      >  3:76613/1‑119 (MQ=255)
                                                                                                                                            ggCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTtcttct       <  3:730691/149‑1 (MQ=255)
                                                                                                                                            ggCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTtcttct       >  1:661212/1‑149 (MQ=255)
                                                                                                                                            ggCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTtcttct       <  3:724769/149‑1 (MQ=255)
                                                                                                                                            ggCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTtcttct       <  2:367390/149‑1 (MQ=255)
                                                                                                                                              cgcTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTAt     <  1:487526/149‑1 (MQ=255)
                                                                                                                                               gcTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTAtt    <  3:492063/149‑1 (MQ=255)
                                                                                                                                                 tGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTa  >  4:406547/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGCAGATTTCACCGGACTGCCAGCGAAATGGCAGCCGATTAATGATTACGGAGCCAAGGTACAGGCGCATGTCATCGACAAATATTGAACAAGTGATGCCGGTTAAGCTGGCGCAGGCGCTGGCGAATCCGTTATTTCCGGCGCTGGACAGCGCCTTACGTTCAGGACGCCATATTGGCCTCGACGAACTGGATAATCATGCATTCCTGATGGATTTTCAGGAATATCTGGAAGAGTTTTACGCGCGTTATAACGTTGAGCTTATTCGCGCACCAGAAGGGTTCTTCTATTTA  >  NC_000913/975553‑975846

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