Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,335,035 T→G N186H (AAC→CAC)  eptA ← lipid A phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,335,0350TG98.5% 247.9 / ‑3.7 68N186H (AAC→CAC) eptAlipid A phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (37/30);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CGTTTTTCGTTCTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTTATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCAATCA  >  NC_000913/4334890‑4335181
                                                                                                                                                 |                                                                                                                                                  
cGTTTTTCGTTCTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGc                                                                                                                                               >  1:896836/1‑151 (MQ=255)
    tttCGTTCTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGa                                                                                                                                            >  2:67998/1‑150 (MQ=255)
    tttCGTTCTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGa                                                                                                                                            >  1:813320/1‑150 (MQ=255)
       cGTTCTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGaaca                                                                                                                                         >  2:685092/1‑150 (MQ=255)
           cTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACg                                                                                                                                      >  2:386164/1‑149 (MQ=255)
           cTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACg                                                                                                                                      >  1:271098/1‑149 (MQ=255)
            tGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAgg                                                                                                                                   >  1:578090/1‑151 (MQ=255)
            tGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAgg                                                                                                                                   >  2:69081/1‑151 (MQ=255)
                  aaCGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAg                                                                                                                              <  1:450467/150‑1 (MQ=255)
                  aaCGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAg                                                                                                                              <  1:708283/150‑1 (MQ=255)
                    cGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCt                                                                                                                           <  2:634721/151‑1 (MQ=255)
                      ggTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCttt                                                                                                                         >  2:893296/1‑151 (MQ=255)
                          gCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTtata                                                                                                                      >  1:762499/1‑150 (MQ=255)
                                cgcgTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGt                                                                                                               <  1:437647/151‑1 (MQ=255)
                                   gTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTgcg                                                                                                            <  1:923317/151‑1 (MQ=255)
                                       tCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTgcggcg                                                                                                         <  2:148604/150‑1 (MQ=255)
                                       tCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTgcggc                                                                                                          <  2:843497/149‑1 (MQ=255)
                                        cACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGAc                                                                                                       >  1:479145/1‑151 (MQ=255)
                                          ccAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGAcc                                                                                                      >  2:317511/1‑150 (MQ=255)
                                           cAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGCGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGAcc                                                                                                      >  1:120052/1‑149 (MQ=255)
                                            aaTTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGc                                                                                                   <  2:134567/151‑1 (MQ=255)
                                            aaTTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGc                                                                                                   <  2:773853/151‑1 (MQ=255)
                                                 gCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAAt                                                                                              >  2:630860/1‑151 (MQ=255)
                                                      aGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCagta                                                                                         <  2:822656/151‑1 (MQ=255)
                                                          gCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTAc                                                                                     <  1:936147/151‑1 (MQ=255)
                                                            aGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAg                                                                                   <  2:355852/151‑1 (MQ=255)
                                                                ttGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGaaa                                                                                <  2:556551/150‑1 (MQ=255)
                                                                   ccAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAAcc                                                                              >  1:449652/1‑149 (MQ=255)
                                                                     aGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGa                                                                           >  2:704390/1‑150 (MQ=255)
                                                                     aGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGa                                                                           >  2:893374/1‑150 (MQ=255)
                                                                       tCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAat                                                                         <  2:712840/150‑1 (MQ=255)
                                                                        cGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAata                                                                        <  1:720426/150‑1 (MQ=255)
                                                                           tGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTgg                                                                    >  2:208068/1‑151 (MQ=255)
                                                                            gATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTg                                                                     <  2:359085/149‑1 (MQ=255)
                                                                            gATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGc                                                                   <  2:778291/151‑1 (MQ=255)
                                                                               gggAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTc                                                                 <  2:505479/150‑1 (MQ=255)
                                                                                     aCCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGa                                                           <  2:834189/150‑1 (MQ=255)
                                                                                      ccATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGaa                                                          >  1:231924/1‑150 (MQ=255)
                                                                                         tGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGaaa                                                         >  2:503341/1‑148 (MQ=255)
                                                                                             cAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAacac                                                  >  2:228627/1‑151 (MQ=255)
                                                                                                 tGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCACCACGGAAAAGAACACTgcg                                              >  1:723040/1‑151 (MQ=255)
                                                                                                  ggCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTgcgc                                             <  1:492304/151‑1 (MQ=255)
                                                                                                   gCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCa                                            <  2:420253/151‑1 (MQ=255)
                                                                                                       caaTGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGAcg                                        >  2:619280/1‑151 (MQ=255)
                                                                                                         aTGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGAcgcg                                      <  1:468976/151‑1 (MQ=255)
                                                                                                           gCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGAcgcg                                      >  2:120154/1‑149 (MQ=255)
                                                                                                              gTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTg                                 <  2:695654/151‑1 (MQ=255)
                                                                                                                 agagGGTCTTAAGGTTTTCACCAGCTGTCTTTGATTTCGGCACAACGTTGAGTAGTCTTTATAACACAGCGCTGCGACCAGCAAAATTAGTACTACAGATACCAGAATATTCGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGc                               >  1:127291/1‑150 (MQ=255)
                                                                                                                   agGGGCTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAgg                            >  2:482941/1‑151 (MQ=255)
                                                                                                                     gggCTTAAGGATTTCACCAGCTCTTTTTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGtt                          <  1:548019/151‑1 (MQ=255)
                                                                                                                        cTTAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTg                        >  2:812139/1‑150 (MQ=255)
                                                                                                                        cTTAAGGATTTCACCAGCTCTTTGTGATTGAGGAACAACGAGGCGAAGTCTTTATACAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTg                        >  1:687294/1‑150 (MQ=255)
                                                                                                                         ttAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGAt                      >  1:769620/1‑151 (MQ=255)
                                                                                                                          tAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGa                       <  2:486968/149‑1 (MQ=255)
                                                                                                                          tAAGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGAtt                     >  2:769922/1‑151 (MQ=255)
                                                                                                                           aaGGATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGAttt                    >  1:770711/1‑151 (MQ=255)
                                                                                                                             ggATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGAttttt                  >  1:410545/1‑151 (MQ=255)
                                                                                                                             ggATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGAttttt                  <  1:416966/151‑1 (MQ=255)
                                                                                                                              gATTTCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGAttttt                  >  2:284836/1‑150 (MQ=255)
                                                                                                                                 ttCACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATcc              >  1:431477/1‑151 (MQ=255)
                                                                                                                                   cACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATcc              >  1:521442/1‑149 (MQ=255)
                                                                                                                                   cACCAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCca             <  2:309537/150‑1 (MQ=255)
                                                                                                                                      cAGCTCTTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCcagcag         <  1:685073/151‑1 (MQ=255)
                                                                                                                                        gCTCTTTGTTATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGc       >  1:889509/1‑151 (MQ=255)
                                                                                                                                           ctTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCa      <  2:435519/149‑1 (MQ=255)
                                                                                                                                           ctTTGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCAAt    >  2:153862/1‑151 (MQ=255)
                                                                                                                                             ttGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCaa     <  2:317862/148‑1 (MQ=255)
                                                                                                                                             ttGTGATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCAATc   >  1:712186/1‑150 (MQ=255)
                                                                                                                                              tgtgATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCAATc   >  2:262810/1‑149 (MQ=255)
                                                                                                                                                 tATGGCGGAACCAGGTGCCGTAGTCTTTATAATACAGTGCGCCTACCAGCAACATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGAGGCAGGTATGATCTTTATTCAGCAGCCAACCa  >  1:176064/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGTTTTTCGTTCTGCATTAACGGGTTGCGGTGCGCGTCTTCACCAATTCGCACCAGCGGCAGATTTGCCAGTCGCTGATGGGAGTACCATGACCAGCTGGCAACAATGCTGTTAGAGGGGCTTAAGGATTTCACCAGCTCTTTGTTATTGCGGAACAACGAGGCGTAGTCTTTATAAAACAGTGCGGCGACCAGCAAAATCAGTAGTACAGAAACCAGAATATTGGCTCCACGGAAAAGAACACTGCGCAGACGCGAGGTGGCAGGTTTGATTTTTATCCAGCAGGCAATCA  >  NC_000913/4334890‑4335181

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